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ABSTRACT 


A  set  of  internal  variable  type  constitutive  equations  which  model  large 
elastic-viscoplastic  deformations  of  metals  at  high  tenyeiatutes  is  formulated. 

For  a  first-order  representation  of  large  deformation  constitutive  behavior  we 
use  a  scalar  as  an  internal  variable  to  represent  the  isotropic  resistance  to 
plastic  flow  offered  by  the  internal  state  of  the  material.  Hot  compression 
tests  on  an  iron-2%  silicon  alloy  have  been  conducted.  Based  on  these  experi¬ 
ments,  specific  viscoplastic  constitutive  functions  are  proposed,  and  the  values 
for  the  material  parameters  appearing  in  these  functions  are  determined.  The 
constitutive  equations  with  these  functions  and  material  parameters  accurately 
reproduce  the  basic  tests  and  also  accurately  predict  the  response  of  the  metal 
to  some  critical  experiments  which  were  not  used  to  determine  the  material  param¬ 
eters  in  the  constitutive  functions.  --Although  the  use  of  a  single  scalar  internal 
variable  limits  the  constitutive  equations  to  describing  the  deformation  behavior 
of  initially  isotropic  metals  up  to  deformation  levels  where  significant  internal 
damage  and  polycrystal  texturing  has  not /developed ,  the  constitutive  equations 
presented  here  should  be  useful  for  obtaining  an  improved  analysis  of  various  hot 
deformation  processing  operations  of  medals. 
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PREFACE 


This  final  report  covers  the  effort  of  the  research  project  titled  “Constitutive  Equa¬ 
tions  For  Hot  Working”  conducted  in  the  Department  of  Mechanical  Engineering  at 
Massachusetts  Institute  of  Technology,  under  U.S.  Army  MTL  Grant  No.  DAAG46-83- 
K-0021.  The  performance  period  for  this  grant  was  5-2-83  through  9-30-86. 

This  report  is  excerpted  from  the  Ph.D.  thesis  of  Dr.  Stuart  Brown,  which  was  finan¬ 
cially  supported  by  this  grant.  The  principal  investigator  on  this  project  was  Professor 
Lallit  Anand  of  the  Department  of  Mechanical  Engineeering  at  MIT,  and  the  program 
manager  for  the  project  at  MTL  was  Dr.  Dennis  Tracey. 
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Chapter  1 

Summary  Of  Results 


We  have  formulated  a  state  variable  constitutive  model  for  large  deformation,  isotropic 
thermo-elasto- viscoplasticity.  The  viscoplastic  material  parameters  appearing  in  the  con¬ 
stitutive  functions  are  determined  from  experimental  data  obtained  from  hot,  isothermal, 
strain  rate  jump,  load-unload-hold-reload,  and  constant  true  strain  rate  experiments  con¬ 
ducted  on  an  Fe-2%  Si  alloy.  The  constitutive  model  is  briefly  described  below.  The 
details  of  our  work  our  reported  in  the  chapters  to  follow. 

Constitutive  Model 

The  state  variables  are  taken  to  be  {T,  0.  s},  where  T  is  the  Cauchy  stress,  0  is  the 
absolute  temperature,  and  s  is  a  scalar  internal  variable  with  dimensions  of  stress,  called 
the  isotropic  deformation  resistance.  The  internal  variable  s  represents  an  aveiaged 
isotropic  resistance  to  macroscopic  plastic  flow  offered  by  the  underlying  “isotropic” 
strengthening  mechanisms  such  as  dislocation  density,  subgrain  size,  grain  size,  solid 
solution  strengthening  effects,  etc. 

The  evolution  equations  for  the  three  state  variables  are: 

•  Evolution  equation  for  the  stress: 

Tv  =  C[D  -  Dp]  -  no, 
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where  =  T  —  WT  +  TW  is  the  Jauinann  derivative  of  Cauchy  stress  T; 
C  =  2 pX  +  {k  —  (2/3)/i}l  ®  1  is  the  fourth  order  isotropic  elasticity  tensor;  p  = 
p(9),K  =  k(0)  are  the  elastic  shear  and  bulk  moduli,  respectively;  II  =  (3kq)1  is 
the  stress-temperature  tensor;  a  =  a(0)  is  the  coefficient  of  thermal  expansion;  D 
is  the  stretching  tensor;  W  is  the  spin  tensor;  X  is  the  fourth  order  identity  tensor; 
and  1  the  second  order  identity  tensor. 

The  plastic  stretching  tensor  Dp  is  given  by  the  flow  rule: 

D”  =  (3/2)f  (T'/S),  with  ?  =  /(( |),«), 

where  T'  denotes  the  stress  deviator,  and  a  =  ^/(3/2)T'  •  T'  denotes  the  equivalent 
tensile  stress. 

•  Evolution  equation  for  s: 

This  is  taken  to  be  given  by 

s  =  h(a ,  0 ,  s)  t  —  r(0,  s ), 

where  h  is  a  hardening  function  which  accounts  for  any  dynamic  recovery,  and  r  is 
a  static  recovery  function  in  the  absence  of  stress. 

•  Evolution  equation  for  0: 

From  an  approximate  balance  of  energy 

0  =  (pc)"1  {div (A grad 0)  +  wT'  Dp} , 

where  p  is  the  mass  density,  c  =  c(0)  is  the  specific  heat.  A  =  A(0)  is  the  thermal 
conductivity,  and  w,  a  scalar  in  the  range  0.85  <  u>  <  1,  is  the  fraction  of  plastic 
work  converted  to  heat. 
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Based  on  our  hot  compression  experiments  on  an  Fe-2%Si  alloy  (performed  in  the 
homologous  temperature  range  0.6  to  0.9  and  the  strain  rate  range  10-3  to  10°  sec-1), 
we  have  proposed  the  following  specific  constitutive  functions  for  f  and  s: 


•~p 

e 


=  A  exp 


-no)  [sinh  Of). 


<r\  1  l/m 


Our  experiments  show  that  for  at  least  the  Fe-2%Si,  the  static  recovery  function  f  has 
negligible  contribution  to  the  evolution  of  s,  and  that  s  may  be  adequately  represented 

by 

i  =  (  h0 


_±)  sign  (l -£)}?, 


with 


The  list  of  material  parameters  in  these  constitutive  equations  are:  A,Q,m,£,ho, 
a,s ,  and  n.  Also  R  is  the  universal  gas  constant.  We  have  formulated  a  systematic 
procedure  to  determine  these  material  parameters  from  data  obtained  from  isothermal, 
strain  rate  jump,  and  constant  true  strain  rate  experiments.  These  material  parameters 
are  listed  below. 
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TABLE  1 


Material  Parameters  for  Fe  -  2%  Si 


Material  Parameter 

Value 

A 

6.346  x  1011  sec-1 

Q 

312.35  kJ/mole 

3.25 

m 

0.1956 

s 

125.1  MPa 

n 

0.06869 

ho 

3093.1  MPa 

a 

1.5 

The  correlation  between  the  model  and  the  experimental  data,  upon  which  the  model 
is  based  and  from  which  the  material  parameters  of  the  model  have  been  determined, 
is  very  good.  Also,  we  have  also  performed  a  variety  of  experiments  that  were  not  used 
to  determine  the  material  parameters  in  our  constitutive  functions.  These  experiments 
were  performed  to  test  the  predictive  capability  of  the  constitutive  model.  We  find  that 
our  constitutive  model  accurately  predicts  the  response  of  the  material  to  these  critical 
experiments. 

It  is  expected  that  specific  forms  of  constitutive  equations  should  be  the  same  for  dif¬ 
ferent  materials,  as  long  as  the  underlying  operative  physical  mechanisms  are  the  same. 
Since  the  mechanisms  are  generally  the  same  for  a  given  class  of  materials  in  a  given  range 
of  temperatures  and  strain  rates,  we  anticipate  that  the  constitutive  equations  proposed 
above  should  be  applicable  for  modeling  the  deformation  behavior  of  other  cubic  metals 
at  temperatures  greater  than  half  the  melting  tempeiature  in  degrees  absolute,  and  in 
the  strain  rate  range  10-3  to  10°  sec-1.  This  range  of  strain  rates  and  temperatures  is 
of  great  practical  interest  for  numerous  low-rate  hot  working  processes  (e.g.,  isothermal 
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forging).  Finally,  although  the  use  of  a  single  scalar  internal  variable  limits  the  con¬ 
stitutive  equations  to  describing  the  deformation  behavior  of  initially  isotropic  mela.s 
upto  deformation  levels  v.here  significant  polycrystal  texturing  has  not  developed,  the 
constitutive  equations  developed  during  the  course  of  our  research  should  be  useful  for 
obtaining  improved  analyses  of  various  hot  deformation  processing  operations  of  metc.ls. 
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Chapter  2 
Introduction 


Hot  working  refers  to  a  collection  of  metal  forming  processes  which  take  advantage  of  the 
lower  flow  stress  and  greater  ductility  of  metals  at  high  temperatures.  It  encompasses 
such  operations  as  hot  rolling,  extrusion,  and  forging,  where  metals  are  heated  to  a  laige 
fraction  of  their  noting  temperatures  and  then  subjected  to  large  deformations.  Hot 
working  is  a  common  operation  in  the  processing  of  metals.  It  has  been  estimated  that 
well  over  80  percent  of  all  metallic  products  undergo  some  form  of  hot  working  during 
their  fabrication  history  [Semiatin  et  al.  1981]. 

The  benefits  of  hot  working  are  extensive.  Large  changes  in  shape  may  be  achieved 
with  each  processing  step.  The  machine  loads  required  to  achieve  these  changes  arc 
much  lower  than  would  be  encountered  at  room  temperature.  There  can  be  substan¬ 
tial  savings  in  mateiial,  with  properly  designed  processes  producing  little  waste  metal. 
There  is  also  the  opportunity  to  control  material  properties  during  the  working  operation. 
Many  hot  working  processes  produce  beneficial  distributions  of  oriented  grains  and  in¬ 
troduce  deformation- induced  microstructures  which  increase  the  subsequent  deformation 
and  fracture  resistance  of  the  material. 

There  are  several  features  common  to  most  hot  working  processes.  Most  hot  work¬ 
ing  produces  large  deformations,  with  equivalent  strains  easily  exceeding  1  or  2.  Most 
processes,  to  achieve  lower  working  loads,  heat  the  worked  metal  to  temperatures  from 
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one-half  to  nine-tenths  its  melting  (homologous)  temperature.  Strain  rates  are  similarly 
high,  ranging  from  10-2  per  second  for  isothermal  forgings  to  103  or  more  per  second 
in  finish  rolling  and  high  speed  hammer  operations.  Many  working  processes  are  in¬ 
terrupted,  incorporating  multiple  deformation  passes  on  the  same  workpiece  with  hold 
periods  between  passes.  Temperatures  may  vary  dramatically,  as  a  function  of  both  time 
and  position  within  the  workpiece.  Temperature  changes  may  occur  due  to  both  heat 
losses  to  tooling  and  the  atmosphere  and  the  conversion  of  plastic  work  into  heat.  Finally, 
the  mode  of  deformation  may  be  extremely  inhomogeneous,  with  very  complex  material 
flow  paths  involving  large  rotations,  strains,  and  changes  in  strain  rate. 

Partially  due  to  the  conservatism  of  the  industry  and  partially  due  to  the  complex¬ 
ity  of  the  hot  working  process,  much  of  the  knowledge  for  the  design  of  hot  working 
processes  has  been  both  experiential  and  proprietary.  The  modelling  of  such  processes 
requires  the  solution  of  very  complex,  strongly  coupled  mechanical  and  thermal  boundary 
value  problems.  The  advent  of  digital  modelling  techniques  such  as  the  finite  element 
method  is  beginning  to  permit  the  simulation  of  working  operations  which  before  were 
frequently  designed  on  a  trial-and-error  basis.  The  problems  associated  with  such  simula¬ 
tions,  however,  are  far  from  resolved.  Several  difficulties  remain,  many  which  encompass 
experimental,  theoretical,  and  computational  disciplines.  These  difficulties  include  the 
proper  accomodation  of  kinematics  of  large  deformation,  the  complexity  of  material  be¬ 
havior  possible  in  hot  working,  coupled  thermo-mechanical  deformation,  the  complex 
friction  and  thermal  boundary  conditions,  and  element  distortion  and  remeshing. 

This  investigation  concentrates  on  the  central  issue  of  appropriate  large  deformation 
constitutive  equations  for  the  high  temperature  deformation  of  metals.  Currently  used 
constitutive  models  (e.g.  in  the  F.E.M.  program  ALPID)  are  usually  simple  three  dimen¬ 
sional  generalizations  of  uniaxial  stress  versus  strain  rate  relationships  of  the  power  law 
form  where  the  stresses  and  strain  rates  are  those  obtained  as  “steady  state”  values  after 
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any  transients  have  decayed  [Thomsen,  et.  al.  1965;  Jonas,  1969,  Sellars  1972,  1978].  Al¬ 
though  these  models  consider  only  steady  state  behavior,  conditions  described  above  for 
hot  working  are  inherently  non-steady.  The  interrupted,  inhomogeneous  deformations 
resulting  from  hot  working  may  prevent  the  numerous  operative  physical  mechanisms 
and  thermally  activated  processes  from  ever  reaching  a  steady  state.  The  final  state  of  a 
hot  worked  piece  of  metal  therefore  is  strongly  affected  by  its  deformation  history.  This 
history  dependence  suggests  a  constitutive  model  of  the  internal  variable  type,  where  the 
current  state  or  condition  of  the  microstructure  is  tracked  via  a  list  of  variables  which 
evolve  with  the  deformation  and  temperature  history.  These  internal  variables  are  not 
necessarily  directly  measureable,  but  they  are  intended  to  represent  some  averaged  ma¬ 
terial  property,  such  as  a  generalized  resistance  to  plastic  deformation.  At  any  point, 
the  internal  variables  are  assumed  to  represent  the  condition  of  the  metal;  there  is  no 
need  to  know  the  prior  thermo- mechanical  history.  A  second  list  of  variables  is  added 
to  the  list  of  internal  variables  to  complete  the  constitutive  model.  These  imposed  or 
“external”  variables  such  as  strain  rate,  stress,  and  temperature  are  assumed,  with  the 
internal  variables,  to  be  sufficient  to  represent  the  current  mechanical  response  of  the 
material  and  to  govern  the  evolution  of  both  the  external  and  internal  variables. 

Several  investigators  have  recently  proposed  internal  variable  constitutive  models  for 
the  high  temperature  deformation  of  metals  [Hart,  1970;  Bodner  and  Partom,  1975; 
Miller,  1976;  Kocks,  1976;  and  Anand  1982].  Virtually  all  models  except  that  of  Anand 
have  been  proposed  for  application  in  the  creep  deformation  regime,  at  temperatures 
and  strain  rates  lower  than  that  encountered  in  hot  working.  The  issues  associated  with 
formulating  an  internal  variable  model  in  either  regime  are  the  same.  However,  due  to  lack 
of  experimental  data  in  the  hot  working  regime,  it  is  not  evident  that  a  model  formulated 
for  creep  may  be  extrapolated  for  hot  working  simulation.  The  model  proposed  by  Anand 
for  hot  working  was  based  on  a  reduced  set  of  moderate  strain,  material  tests  representing 
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a  small  subset  of  the  hotworking  regime.  Anand  by  necessity  assumed  functional  forms 
for  his  material  model  which  he  then  fit  to  the  experimental  data.  No  data  were  available 
fo-  model  validation  independent  of  that  used  for  the  fitting  of  model  parameters. 

The  goals  of  the  research  reported  here  were  to: 

1.  Develop  a  high  temperature  mechanical  test  system  to  characterize  metal  behavior 
in  the  hot  working  regime,  including  high  homologous  temperatures,  moderate 
strain  rates,  and  large  strains. 

2.  Perform  experiments  necessary  to  investigate  the  constitutive  response  of  a  repre¬ 
sentative  metal  in  the  hot  working  regime. 

3.  Formulate  an  internal  variable  constitutive  model  from  the  experimental  data,  with 
the  following  considerations: 

a.  The  model  should  represent  large,  three  dimensional  deformations  in  the  hot 

working  regime. 

b.  The  model  should  be  able  to  represent  the  most  significant  aspects  of  the  under¬ 

lying  microstructural  state  of  the  metal  and  its  evolution  during  deformation. 
Physical  phenomena  to  be  modelled  include  strain  hardening,  strain  rate  and 
temperature  sensitivity  including  history  effects,  and  both  dynamic  and  static 
recovery  processes. 

c.  Material  parameter  determination  should  be  straightforward  and  require  a  min¬ 

imum  of  experimental  testing.  This  should  of  course  be  commensurate  with 
the  degree  of  complexity  that  is  to  be  modelled. 

d.  Finally,  the  model  should  be  formulated  with  due  consideration  of  issues  regard¬ 

ing  its  numerical  implementation  in  finite  element  programs. 
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4.  Evaluate  the  model  through  appropriate  validation  experiments,  including  an  in¬ 
vestigation  of  the  model’s  ability  to  predict  some  aspects  of  the  internal  state  of 
the  material. 

The  following  chapter  describes  a  high  temperature  test  system  and  describes  the 
compression  testing  procedures  associated  with  an  iron  -  2%  silicon  alloy.  Chapter  4 
discusses  the  structure  of  the  constitutive  model,  describes  a  set  of  experiments  used  to 
evaluate  the  material  functions,  and  proposes  specific  forms  for  those  functions.  Chapter 
5  discusses  the  reduction  of  the  data  obtained  from  the  mechanical  tests  and  presents 
a  procedure  for  determining  material  constants  for  the  constitutive  model.  Chapter  6 
presents  a  comparison  of  material  response  predicted  by  the  constitutive  model  with  ex¬ 
periments  involving  boundary  conditions  not  used  for  the  material  parameter  evaluation. 
This  chapter  also  discusses  means  of  correlating  predicted  values  of  the  internal  variable 
with  physically  measureable  quantities.  Finally,  we  conclude  this  report  with  some  final 
remarks  and  discussion  in  Chapter  7. 
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Chapter  3 

Experimental  Apparatus 


This  chapter  begins  with  a  description  of  the  high  temperature  test  system  designed 
and  used  to  evaluate  material  response  under  hot  working  conditions.  Issues  related  to 
compression  testing  are  discussed,  and  the  model  material  used  in  this  investigation,  an 
iron  -  2%  silicon  alloy,  is  described. 

3.1  High  Temperature  Test  System 

Any  system  used  for  basic  high  temperature  mechanical  testing  and  for  the  simulation 
of  hot  working  should  possess  the  following  minimal  characteristics: 

1.  Isothermal  testing  to  temperatures  exceeding  .9  0m  of  the  material  of  interest,  where 
0m  is  the  melting  temperature  in  degrees  absolute. 

2.  True  strain  rates  ranging  from  10“4  to  at  least  10°  sec-1  and  preferably  greater. 

3.  Variable  testing  modes,  including  compression,  torsion,  and  tension. 

4.  Rapid  quenching  of  specimens  from  test  temperatures. 

5.  Complex  loading  histories,  both  to  evaluate  material  response  and  to  simulate  non¬ 
steady  forging  operations. 
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6.  Controllable  environments  to  reduce  the  effect  of  oxidation. 


7.  High  speed  data  collection. 

Figure  3.1  shows  a  photograph  of  the  test  system  constructed  to  meet  the  above 
characteristics.  The  basic  system  consists  of  an  Instron  Corporation,  Model  1322  ten¬ 
sion/torsion  servohydraulic  test  machine.  The  system  actuator  has  a  frequency  response 
of  approximately  10  hertz,  which  effectively  limits  the  maximum  controllable  strain  rate 
for  ordinary,  laboratory  sized  axial  specimens  to  approximately  1  per  second.  The  lower 
strain  rate  limit  is  approximately  10-4  per  second  due  to  the  limits  of  resolution  for  the 
actuator  displacement  and  hydraulic  servovalve.  Approximately  the  same  strain  rates 
are  obtainable  in  torsion. 

The  analog  controllers  for  the  servohydraulic  machine  permit  external  control  signals 
in  the  form  of  voltages.  These  may  be  provided  either  by  digital  or  analog  function 
generators,  permitting  arbitrarily  complex  loading  histories.  Most  of  the  control  signals 
for  this  investigation  were  obtained  from  an  analog  function  generator,  designed  and 
assembled  by  Mehrdad  Haghi  (of  our  laboratory),  following  a  simpler  design  used  by 
Immarigeon  and  Jonas  [Immarigeon  and  Jonas  1975).  The  function  generator  uses  a 
simple  RC  circuit  to  provide  a  decaying  exponential  voltage  which  simulates  a  constant 
true  strain  rate  compression  test.  The  function  generator  is  also  designed  with  set  points 
which  may  either  hold  the  voltage  at  that  point  or  change  the  strain  rate,  which  enables 
strain  rate  jump  and  recovery  tests. 

The  temperature  and  environmental  control  for  the  test  system  is  provided  by  a  high 
temperature,  vacuum  furnace  built  by  Centorr  Associates,  Suncook,  New  Hampshire, 
according  to  our  specifications.  The  furnace,  pictured  in  Figures  3.2  through  3.4  is  a 
stainless  steel,  double-walled  chamber  with  a  center  heat  zone.  The  furnace  is  mounted 
in  the  frame  of  the  servohydraulic  machine,  as  illustrated  in  Figure  3.1.  The  vacuum 
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system  consists  of  a  roughing  pump  and  a  4  inch  diffusion  pump  which  is  attached  to 
the  back  plate  of  the  furnace.  Maximum  vacuum  attainable  was  approximately  10"6 
torr.  The  furnace  is  equipped  with  flexible  bellows  above  and  below  through  which  pass 
loading  rods  composed  of  a  high  temperature  molybdenum  alloy,  TZM.  The  loading  rods 
pass  into  the  furnace  through  the  bellows  and  are  hermetically  sealed  with  a  pair  of 
differentially  pumped  Viton  quad  rings.  The  lower  seal  is  designed  to  permit  the  lower 
loading  rod  to  both  rotate  and  move  in  and  out  of  the  furnace  without  disturbing  the 
vacuum. 

Heating  is  provided  by  a  cylindrical  heating  zone  consisting  of  tungsten  mesh  heat¬ 
ing  elements  surrounded  by  a  multiple  layer  set  of  molybdenum  heat  shields.  Power 
is  provided  by  an  SCR-controlled  A/C  power  supply  which  caused  extensive  problems 
with  electrical  noise.  Electrical  noise  is  a  continual,  unavoidable  problem  with  SCR  con¬ 
trollers.  In  retrospect  the  system  should  have  been  designed  with  a  D/C  power  supply. 
A  Leeds  and  Northrup  Electromax  V  temperature  controller  is  used  with  chromel-alumel 
thermocouples  to  maintain  specimen  temperatures.  A  thermocouple  spot  welded  to  each 
specimen  serves  as  the  feedback  transducer  to  the  controller.  The  heat  zone  is  designed 
to  reach  temperatures  up  to  1800  Celsius,  but  the  loading  rod  materials,  94%  alumina 
and  TZM,  begin  to  creep  and  recrystallize,  respectively,  at  approximately  1200  Celsius. 
All  tests  were  therefore  run  at  or  below  1200  Celsius.  The  furnace  is  water-cooled,  water 
being  circulated  between  the  double  walls  of  the  furnace  while  at  temperature. 

The  furnace  was  also  equipped  with  an  insulated  Conflat  induction  heating  feedthrough. 
Induction  heating  in  atmosphere  or  vacuum  is  therefore  another  possible  means  of  reach¬ 
ing  high  temperatures.  Tests  involving  quenching  require  induction  heating  since  the 
tungsten  mesh  heating  elements  and  molybdenum  heat  shields  can  not  withstand  the 
rapid  contraction  associated  with  querch  cooling.  Induction  heating  possesses  the  same 
disadvantage  as  the  SCR  heating  system  in  that  it  generates  a  significant  amount  of 
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electrical  noise,  obscuring  transducer  signals. 

Other  features  of  the  furnace  include  ports  to  accomodate  gas  quenching,  12  ther¬ 
mocouple  feedthroughs,  and  mounting  brackets  and  feedthroughs  appropriate  for  a  high 
temperature,  MTS  extensometer. 

Data  acquisition  for  most  of  this  project  was  accomplished  with  an  IBM  PC/XT 
equipped  with  a  Metrabyte  Corporation  DASH- 16  analog- to-digital  converter.  The  DASH- 
16  was  provided  with  BASIC  language  callable  subroutines  which  permitted  sampling 
rates  exceeding  1  kilohertz.  An  approximate  limit  of  30,000  data  points  could  be  col¬ 
lected  during  one  sampling  session,  the  maximum  number  controlled  by  the  memory 
remaining  after  subtracting  the  memory  requirements  of  the  BASIC  language  and  the 
BASIC  sampling  program  from  a  64I<  sector  of  IBM/PC  memory.  All  data  was  trans¬ 
ferred  to  a  Data  General  MV4000  computer  for  analysis;  the  IBM/PC  was  used  only  for 
data  acquisition. 

The  high  temperature  system  is  unique  in  the  variety  of  high  temperature  test  ca¬ 
pabilities.  Other  systems  exist  which  permit  one  mode  of  testing,  heating,  or  control, 
but  we  know  of  no  other  system  with  the  same  collected  capabilities.  All  of  the  data 
presented  in  this  report  were  obtained  using  this  test  system. 

3.2  Compression  Testing 

The  role  of  the  compression  test  in  the  study  of  the  large  deformation  of  metals  is  well  es¬ 
tablished  [ASTM  Standard  E209-65;  Chait  and  Papirno,  1983;  Hsu,  1969].  Compression 
testing  achieves  large  strains  due  to  the  absence  of  necking  and  permits  a  fairly  simple 
means  of  measuring  strain  by  determining  the  relative  displacement  of  the  compression 
platens.  Also  it  is  also  easier  to  perform  relative  to  torsion  testing  since  it  requires  sim¬ 
pler  and  more  common  experimental  apparatus  and  requires  less  complicated  alignment 
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of  loading  fixtures.  Compression  testing  is  particularly  appropriate  for  large  deformation 
testing  at  high  temperatures,  where  temperature  uniformity  is  difficult  to  maintain  in 
tension  and  torsion  specimens,  and  strain  measurement  in  tension  and  torsion  requires 
complex  extensometry. 

Unlike  tension  and,  to  a  lesser  extent,  torsion  testing,  compression  testing  frequently 
is  not  selected  for  material  testing  due  to  its  propensity  for  different  modes  of  inhomo¬ 
geneous  deformation.  Temperature  gradients,  lateral  buckling,  and  friction  between  the 
compression  specimen  and  compression  platen  all  may  create  inhomogeneous  deforma¬ 
tion. 

When  a  material  demonstrates  substantial  rate  sensitivity,  temperature  gradients, 
specifically  axial  gradients,  may  produce  conical  specimens,  illustrated  in  Figure  3.5. 
The  gradient  in  this  case  was  caused  by  the  quenching  of  the  upper  end  of  the  specimen 
by  a  cooler  compression  platen.  Extensive  finite  element  heat  transfer  analyses  of  the 
Centorr  furnace  heat  zone  and  the  initially  installed  TZM  alloy  loading  rods  prr  rided  with 
the  furnace  indicated  that  the  greatest  source  of  heat  loss  was  through  the  water-cooled 
loading  rods.  The  conduction  through  the  rods  was  enough  to  p-oduce  temperature 
gradients  in  excess  of  20  degrees  Celsius  along  the  length  of  the  compression  specimens. 
This  gradient  was  reduced  to  ±2  degrees  Celsius  by  installing  one  inch  diameter,  94% 
alumina  ceramic  rods  in  place  of  the  TZM  loading  rods.  The  alumina  rods  fit  within  the 
mounting  sockets  of  the  TZM  rods,  which  were  withdrawn  from  the  furnace  heat  zone. 
The  change  to  alumina  rods  reduced  the  thermal  conductivity  of  the  heat  zone  loading 
rod  material  from  approximately  100  W/m  °K  to  5  W/m  °K.  The  new  loading  rod  system 
is  shown  in  Figure  3.6.  TZM  was  still  used  as  the  compression  platen  material,  since  it 
could  be  easily  polished  to  provide  a  smooth  compression  surface,  and  since  it  acted  to 
distribute  the  specimen  contact  load  evenly  over  the  alumina  rod  crossection.  It  should 
be  noted  that  a  much  easier  method  to  eliminate  temperature  gradients  is  to  employ 
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a  three  zone  furnace.  Such  a  furnace  eliminates  any  sensitivity  to  load  train  materials, 
permitting  the  use  of  high  temperature  materials,  such  as  thoriated  tungsten  or  graphite 
without  concern  for  these  materials’  high  thermal  conductivities. 

Lateral  buckling  was  encountered  in  compression  specimens  whose  height  to  diameter 
ratio  exceeded  1.5.  This  height-to- diameter  limit  appears  rather  inflexible,  for  we  fre¬ 
quently  obtained  shear  buckling  in  specimens  with  a  1.6  ratio.  A  ratio  of  1.5  eliminated 
any  ostensible  buckling. 

A  more  significant  deformation  inhomogeneity  resulted  from  friction  between  the 
specimen  and  the  compression  platens.  Figure  3.7  indicates  the  deformation  pattern 
which  results  from  excessive  friction  between  the  specimen  and  platens.  Figure  3.8  illus¬ 
trates  the  shape  of  a  comparably  barrelled  specimen.  The  effect  of  friction  in  preventing 
homogeneous  compression  was  reduced  through  the  use  of  grooves  in  the  ends  of  the 
compression  specimens  in  conjunction  with  high  temperature  lubricants.  This  technique 
has  been  used  by  many  investigators  [Uvria,  1968;  Sherby,  1980;  Hsu,  1969]  and  appears 
to  work  very  well  in  producing  homogeneous  compression  for  compressive  true  strains 
exceeding  —1.0.  Shallow,  concentric  grooves  are  machined  on  the  end  faces  of  the  spec¬ 
imens  to  hold  high  temperature  lubricant.  Figure  3.9  provides  the  dimensions  of  the 
specimen  and  the  grooves  on  the  specimen  faces.  The  specimen  dimensions  were  selected 
to  provide  sufficient  grains  in  the  cross-section  to  permit  a  continuum  assumption,  the 
number  of  grains  encountered  across  the  diameter  being  approximately  fifteen. 

The  lubricants  used  were  mixtures  of  powdered  glass  and  boron  nitride  powder.  This 
particular  combination  of  glass  and  boron  nitride  appears  to  act  as  a  high  temperature 
analog  to  the  common  room  temperature  lubricants  of  molybdenum  disulfide  or  graphite 
in  a  grease  carrier.  At  high  temperatures  the  glass  melts  to  form  a  viscous,  lubricating 
film  which  also  carries  the  boron  nitride.  Boron  nitride  possesses  a  hexagonal  basal  plane 
crystalline  structure,  providing  lubrication  in  the  same  manner  as  graphite  or  molybde- 
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num  disulfide  [Niedenzu,  et  al.  1965].  The  particular  glass/boron  nitride  mixtures  used 
as  lubricants  at  the  different  test  temperatures  are  listed  in  Table  3.1.  Our  general  ex¬ 
perience  indicates  that  a  good  criterion  for  lubricant  selection  is  that  the  melted  glass 
should  have  a  viscosity  of  approximately  104  poise  at  test  temperature,  and  that  the 
glass-to-boron  nitride  weight  ratio  should  be  greater  than  4-to-l,  and  should  probably 
be  closer  to  8  or  10-to-l. 

It  is  also  possible  to  eliminate  the  boron  nitride  completely  and  still  obtain  uniform 
deformation  using  only  glass.  Boron  nitride  extends  the  effective  usable  lubrication  tem¬ 
perature  range  of  a  particular  glass.  Boron  nitride  may  also  react  with  oxygen  at  high 
temperatures,  so  its  usefulness  may  be  limited  to  vacuum  or  inert  environments. 

Figure  3.10  shows  a  typical  compression  specimen  before  and  after  an  essentially 
homogeneous,  isothermal,  constant  true  strain  rate  compression  test  to  a  true  strain  of 
-100%.  Notice  that  the  deformed  specimen  is  still  cylindrical,  and  that  there  is  no  notice¬ 
able  rollover  of  the  sides  of  the  specimen,  which  would  have  resulted  were  there  friction 
between  the  specimen  end  faces  and  the  compression  platens.  Figure  3.11  shows  the 
microstructure  of  one  such  specimen.  Notice  that  all  of  the  grains  possess  approximately 
the  same  aspect  ratio  as  the  deformed  specimen. 

Additional  factors  influenced  deformation  homogeneity  of  the  compression  specimens. 
It  appears  to  be  very  important  for  the  compression  platens  to  be  well  polished.  The  TZM 
platens  used  for  the  experiments  described  within  were  polished  for  each  experiment  to  an 
approximately  10  micron  surface  finish.  The  ends  of  the  compression  specimens  were  also 
finished  to  a  600  grit  roughness.  It  is  also  important  for  the  opposing  surfaces  of  the  load 
train  to  be  parallel  and  remain  so  during  the  test.  The  end  of  the  compression  specimens 
should  also  be  parallel.  Groove  spacing  and  depth  appeared  to  be  less  sensitive  controlling 
parameters  than  was  initially  anticipated.  Approximate  bounds  for  these  dimensions  are 
a  groove  spacing  of  fifteen  to  thirty  thousandths  of  a  inch  and  a  groove  depth  of  four  to 
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eight  thousandths  of  an  inch. 

Compressive  strains  throughout  this  investigation  were  measured  by  subtracting  the 
effect  of  test  machine  compliance  from  the  displacement  of  the  servohydraulic  actuator. 
The  large  displacements  required  for  large  strain  compression  tests  permitted  this  indirect 
means  of  measuring  deformation,  since  the  errors  associated  with  variations  and  nonlin¬ 
earities  in  compliance  and  the  resolution  of  the  actuator  LVDT  were  relatively  small.  The 
errors  were  more  significant,  however,  at  the  beginning  of  each  test  when  the  tolerances 
were  taken  up,  and  when  we  desired  very  accurate  strain  measurements.  Extensive  digi¬ 
tal  smoothing  techniques  permitted  greater  resolution  in  measuring  displacements  than 
was  possible  using  analog  data  collection  equipment.  Using  smoothing,  we  could  easily 
resolve  strains  less  than  10~4,  although  the  degree  of  resolution  was  sensitive  to  both 
strain  rate  and  rate  of  data  acquisition.  More  exact  resolution  of  strains  would  be  aided 
by  some  means  of  measuring  the  relative  displacement  of  the  compression  platens,  such 
as  a  high  temperature  LVDT,  extensometer,  or  strain  gauge.  Implementation  of  these 
transducers  is  not  trivial,  although  they  are  necessary  for  small  strain  resolution. 

3.3  Iron  -  2%  Silicon  Model  Material 

An  iron  silicon  alloy  was  selected  as  the  model  material  in  this  investigation  for  several 
reasons.  First,  although  an  iron  based  material,  it  retains  a  body-centered  cubic  structure 
up  to  the  solidus  temperature.  This  property  is  useful,  for  we  wished  to  be  able  to 
deform  our  model  material  at  hot  working  temperatures  and  then  rapidly  quench  the  test 
specimens  to  preserve  the  hot  worked  microstructure.  The  martensitic  transformation 
associated  with  many  iron  alloys  would  erase  any  deformation-induced  microstructure. 
The  equilibrium  diagram  for  the  iron/silicon  system  shown  in  Figure  3.12  indicates  the 
persistence  of  the  alpha  phase  at  all  temperatures  for  the  above  composition. 
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The  iron  2%  silicon  was  also  selected  for  the  high  stacking  fault  energy  associated 
with  body-centered  cubic  materials.  Our  intention  was  to  model  the  evolution  of  struc¬ 
ture  including  the  mechanisms  of  hardening,  dynamic  recovery,  and  static  recovery.  We 
wished  to  avoid  dynamic  recrystallization,  which  can  introduce  such  complications  as  an 
oscillating  stress/strain  response  and  deformation  localization  [Jonas,  1969].  Figure  3.13 
illustrates  the  oscillations  which  result  during  the  deformation  of  a  material  which  ex¬ 
hibits  dynamic  recrystallization,  in  this  case  1018  plain  carbon  steel.  High  stacking  fault 
energy  materials  generally  do  not  demonstrate  dynamic  recrystallization  since  the  rate  of 
recovery  is  believed  to  be  sufficiently  high  to  prevent  the  accumulation  of  a  deformation 
structure  which  may  trigger  recrystallization  [McQueen,  1982]. 

The  deformation  response  of  iron  silicon  alloys  has  also  been  studied  extensively 
using  both  transmission  electron  microscopy  [e.g.  Hu,  1964]  and  dislocation  etch  pits 
[Lytton,  et  al.,  1965].  We  felt  that  extensive  experience  in  the  examination  of  iron  silicon 
microstructures  would  facilitate  efforts  at  similar  examinations  should  we  wish  to  do  so. 

The  particular  alloy  used  in  this  investigation  was  provided  by  Armco  Steel  Company, 
Middletown,  Ohio.  It  was  provided  in  rolled  plates,  approximately  1  inch  thick,  12  inches 
wide,  and  3  to  4  feet  long.  The  nominal  composition  is  listed  in  Table  3.3.  All  specimens 
were  machined  such  that  their  axis  of  symmetry  was  oriented  across  the  width  of  the 
plate,  illustrated  in  Figure  3.14.  The  material  possessed  a  large  grain  size.  Figure  3.15 
illustrates  the  polished  and  etched  (Fry’s  reagent)  cross-section  of  a  typical  compression 
specimen. 

All  specimens  were  annealed  for  1  hour  at  700  degrees  Celsius.  This  schedule  is 
identical  to  that  used  by  Young  and  Sherby  [Young  and  Sherby,  1973]  in  an  investigation 
of  a  similar  iron  silicon  alloy.  Annealing  at  higher  temperatures  and  for  a  longer  time 
was  avoided  to  prevent  excessive  grain  growth. 
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TABLE  3.1 

Lubricant  Composition 


Temperature 

(Celsius) 

Constituents 

Composition  Ratio 
by  Weight 

700 

BN:1190:Acetone 

1:4:16 

800 

BN:1190:Acetone 

1:4:16 

900 

BN:0010:Acetone 

1:8:32 

1000 

BN:0010:Acetone 

1:8:32 

1100 

BN:0010:Acetone 

1:8:32 

1200 

BN:0010:Acetone 

1:8:32 

Notes: 


1.  BN  =  boron  nitride  powder,  purchased  from  Union  Carbide  Corporation,  Grade  HCP. 

2.  Four  digit  number  refer  to  glass  classification  codes.  Glass  was  purchased  as  a  powder, 
325  screen  size,  from  Corning  Corporation,  Corning,  New  York. 
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TABLE  3.2 

Composition  of  Iron  -  2%  Silicon 
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gure  3.1  High  temperature  mechanical  test  system. 
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Figure  3.2  Vacuum  furnace  exterior 


Figure  3.3  Vacuum  system  for  high  temperature 
furnace. 
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Figure  3.4  Vacuum  furnace  interior. 


Figure  3.5  Conical  specimen  resulting  from  axial 
temperature  gradient. 


Figure  3.6  Load  train  for  isothermal  compression 

testing. 


Figure  3.7  Pattern  of  grain  deformation  in 
compression  specimen  with  significant  friction  between 
specimen  and  compression  platen. 


Figure  3.8  Barrelled  compression  specimen  resulting 
from  friction  between  specimen  and  compression 
platen. 
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0.750 


Figure  3.9 


Dimensions  of  Fe  -  2%  Si  Compression  Specimen. 
Test  results  insensitive  to  groove  geometry. 

(All  dimensions  in  inches.) 
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Figure  3. 10  Compression  specimen  before  and  after 
homogeneous,  isothermal,  constant  true  strain  rate 
compression  test. 


Figure  3  J 1  Grain  deformation  obtained  from  a 
homogeneous,  isothermal,  constant  true  strain  rate 
compression  test. 
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Fe-Si  Iron-Silicon 


Figure  3.12  Equilibrium  diagram  of  iron/silicon 
system. 
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Strain 

Figure  3.13  Isothermal,  constant  true  strain  rate  compression  data 
for  1018  carbon  steel.  Oscillations  are  due  to 
recrystallization. 
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Rolling  direction 


<■ 


Compression  specimen 
axis 


Figure  3.14  Orientation  of  compression  specimens 
relative  to  rolled  Fe  -  2%  Si  plate. 
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Figure  3. 1 5  Crosssection  of  compression  specimen 
etched  to  indicate  grain  size. 
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Chapter  4 

Experimental  Determination  of  the 
Material  Response  Functions 


The  basic  form  of  the  constitutive  model  follows  the  mathematical  structure  proposed  by 
Anand  for  an  isotropic  metal  [Anand,  1985].  The  model  assumes  three  state  variables: 

{T, «,<>}, 

where  T  is  the  Cauchy  stress,  0  is  the  absolute  temperature,  and  s  is  a  scalar  inter¬ 
nal  variable  which  represents  an  isotropic  resistance  to  plastic  deformation.  Evolution 
equations  for  the  three  variables  are: 

•  Evolution  equation  for  the  stress: 

Tv  =  £[D  -  Dp]  -  HO, 

where 


TV  =  t  _  WT  +  TW 

C  =  2fj,X  -f  {k  —  (2/3), u}l  <8  1 

/i  =  /t(0),  k  =  k(0) 

II  =  (3kqt)1 
or  =  a(0) 


Jaumann  derivative  of  Cauchy  stress  T; 
fourth  order  isotropic  elasticity  tensor; 
shear  and  bulk  moduli; 
stress- temperature  tensor; 
coefficient  of  thermal  expansion; 
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L  =  grad  v 
D  =  sym  (L) 
W  =  skew  (L) 


0 

X 

1 


spatial  gradient  of  velocity; 
stretching  tensor; 
spin  tensor; 
absolute  temperature; 
fourth  order  identity  tensor; 
second  order  identity  tensor. 


The  constitutive  equation  for  Dp  is: 


where 


eP  =  f(*As)  >  0, 
T' 

a=J(3/2)T-T 


a  <  s,  equivalent  plastic  tensile  strain  rate; 
deviator  of  the  Cauchy  stress; 
equivalent  tensile  stress. 


•  Evolution  equation  for  the  internal  variable  s : 

s  =  s{a,0,s), 


•  Evolution  equation  for  the  absolute  temperature: 

Given  by  the  following  approximate  energy  balance  equation  (this  is  not  a  consti¬ 
tutive  equation): 

0  =  (pc)-1  {div  (A  grad#)  +  u  T  •  Dp} , 

Here  p  =  p(0)  is  the  mass  density,  e  —  c(0)  is  the  specific  heat,  A  =  A (0)  is  the 
thermal  conductivity,  and  u>,  a  scalar  in  the  range  0.85  <  u>  <  1.0,  is  the  fraction 
of  plastic  work  converted  to  heat. 
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The  main  task  of  this  investigation  then  is  to  evaluate  the  functional  forms  for  the 
equivalent  plastic  strain  rate  f(a ,  0 ,  s)  and  the  evolution  equation  for  the  internal  variable 
s:  s(a,d,s). 

This  constitutive  model  employs  only  one  scalar  internal  variable,  s,  for  several  rea¬ 
sons.  First,  we  wished  to  determine  whether  one  scalar  internal  variable  would  be  suf¬ 
ficient  to  describe  the  major  features  of  metal  behavior  during  hot  working.  The  well- 
known  correlation  between  the  dislocation  density  and  flow  stress  suggested  that  a  single 
internal  variable  might  be  adequate.  There  is  additionally  sufficient  uncertainty  of  what 
role  a  second  internal  variable  would  play.  Some  investigators  [Kocks,  1966]  suggest 
that  obvious  secondary  parameters  such  as  dislocation  cell  size  may  not  contribute  to 
deformation  resistance  and  are  merely  consequences  of  dislocation  density. 

Second,  the  task  of  determining  the  evolution  of  internal  variables  is  much  more 
straightforward  if  there  is  only  one.  We  wished  to  be  able  to  perform  tests  which  would 
suggest  the  functional  form  for  the  evolution  equation  for  the  internal  variable.  More 
than  one  internal  variable  would  necessitate  tests  which  would  somehow  decouple  the 
evolution  of  the  two  or  more  internal  variables.  We  felt  that  the  extent  of  knowledge 
concerning  the  development  of  state  during  large  deformations  provides  very  little  guid¬ 
ance  for  performing  such  experiments. 

The  internal  variable  was  selected  to  be  scalar  since  the  task  of  accomodating  the 
largest  source  of  anisotropy,  deformation-induced  texture,  is  still  not  well  understood. 
Small  strain  anisotropy  associated  with  Bauschinger  type  effects  were  ignored  since  its 
effect  was  assumed  to  saturate  at  a  small  fraction  of  the  strains  encountered  in  hot 
working. 


39 


Here,  we  focus  our  attention  on  the  two  undetermined  response  functions  of  the 
constitutive  model,  the  rate  or  flow  equation: 

e=f(v,s,9),  4.1 

and  the  evolution  equation: 

s=g(a,s,6).  4.2 

Given  the  above  structure,  we  now  must  deduce  the  functional  forms  of  the  two  equa¬ 
tions.  Ideally,  we  would  like  to  have  an  adequate  understanding  of  the  underlying  mecha¬ 
nisms  governing  high  temperature  behavior  such  that  our  task  would  be  only  to  perform 
the  experiments  to  determine  the  few  specific  material  constants  for  the  metal  of  in¬ 
terest.  Realistically,  our  limited  knowledge  of  these  mechanisms  precludes  this,  both 
because  our  understanding  of  the  different  mechanisms  is  incomplete  and  because  the 
actual  mechanisms  are  associated  with  complex  structures  and  processes  certainly  not 
totally  representable  by  a  single  internal  variable  model.  We  therefore  choose  to  moti¬ 
vate  functional  forms,  where  possible,  from  elementary  assumptions  about  deformation 
mechanisms,  guided  by  appropriate  experimentation  to  modify  these  forms  and  provide 
material  parameters. 

We  make  the  following  assumptions  concerning  the  forms  of  equations  4.1  and  4.2: 

1.  We  define  the  scalar  internal  variable  $  to  be  a  stress  valued  quantity  called  the  de¬ 
formation  resistance.  As  a  structure  parameter,  it  represents  a  generalized  isotropic 
resistance  to  plastic  flow  which  contains  contributions  from  dislocation  density,  sub¬ 
grain  size,  grain  size,  solid  solution  strengthening,  intrinsic  lattice  resistance,  etc. 
For  pure,  single  phase  materials,  we  expect  it  to  be  highly  correlated  with  the 
dislocation  density  and  subgrain  size,  which  are  expected  to  be  the  primary  contri¬ 
butions  to  deformation  resistance  at  these  temperatures. 
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2.  The  state  variable  enters  into  the  rate  equation  only  as  a  ratio  with  the  equivalent 
tensile  stress.  That  is, 

?  =  4.3 

This  form  has  been  suggested  by  Kocks,  Argon,  and  Ashby  [1975]  in  conjunction 
with  a  model  for  mechanically-activated  dislocation  glide.  Rice  [1970]  and  others 
have  suggested  this  structure  as  well,  although  there  seems  to  be  no  fundamental 
requirement  for  a  scalar  internal  variable  to  be  incorporated  in  this  way.  We  shall 
consider  functional  forms  for  which  equation  4.3  may  be  inverted  such  that  we 
may  write  a  =  cs,  where  c  is  a  function  of  strain  rate  and  temperature.  This 
structure  parallels  an  expanding  or  shrinking  isotropic  yield  surface,  where  the 
plastic  resistance  acts  as  the  flow  strength.  Of  course,  the  model  presented  here 
does  not  include  a  yield  surface;  plastic  deformation  is  assumed  to  occur  at  any 
non- zero  value  of  equivalent  stress. 

3.  The  evolution  equation  for  the  internal  variable  is  assumed  to  be  of  the  form: 

s  =  h(cr,s,0)t  —  r(s,0).  4.4 

Assuming  three  primary  mechanisms,  hardening,  dynamic  recovery,  and  static  re¬ 
covery,  we  associate  h(cr,s,0)  with  dynamic  processes,  i.e.,  hardening  and  dynamic 
recovery.  Static  recovery  is  accommodated  through  r(s,  0).  This  form  is  reminis¬ 
cent  of  the  Bailey-Orowan  equation  [Bailey,  1926,  Orowan,  1945],  and  is  almost 
universally  employed  in  internal  variable  evolution  equations  for  high  temperature 
deformation.  The  assumption  that  the  static  recovery  function  is  independent  of 
stress  reflects  our  desire  for  this  function  to  represent  evolution  processes  which 
occur  in  unloaded  material. 

Given  the  above  assumptions  about  the  forms  of  the  rate  and  evolution  equations, 
we  now  ask  what  experiments  may  be  performed  to  complete  their  specification.  Given 
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that  the  two  differential  equations  4.3  and  4.4  are  coupled,  it  seems  reasonable  that  we 
should  seek  some  way  of  decoupling  the  two  equations  so  that  the  functional  forms  may 
be  investigated  separately.  Given  two  equations  of  the  form  4.3  and  4.4,  then  following 
a  similar  development  presented  by  Ruina  [1983],  we  may  make  the  statements: 

•  The  value  of  the  internal  variable  is  uniquely  determined  if  all  of  the  external 
variables  (a,6,eP)  are  known.  This  is  portrayed  in  Figure  4.1  where  dotted  lines  of 
constant  internal  variable  s  are  indicated  in  stress/strain  rate/temperature  space. 
The  lines  are  straight  only  for  illustration,  there  is  no  physical  reason  why  they 
should  be  straight.  The  dashed  line  exists  for  those  materials  which  reach  some 
saturation  condition  where  the  value  of  the  internal  variable  attains  some  constant 
value,  here  denoted  as  s*. 

•  Experiments  where  we  apply  fast  changes  in  the  external  variables  (so  that  we 
may  assume  constant  structure)  should  produce  the  relationship  between  them 
expressed  by  equation  4.3.  If  we  are  confident  that  the  structure  has  not  changed 
in  the  time  necessary  to  apply  the  change,  then  the  changes  in  the  external  variables 
are  governed  only  by  the  rate  equation.  This  then  suggests  a  series  of  experiments 
where  we  rapidly  change  either  the  strain  rate,  temperature,  or  stress  to  different 
values,  with  the  values  of  the  other  external  variables  being  the  same  just  before  the 
change.  This  guarantees  that  the  initial  state  is  unchanged,  and  the  rate  equation 
is  decoupled  from  the  evolution  equation. 

Using  the  assumptions  listed  above  and  exploiting  the  characteristics  of  the  functional 
forms  outlined,  the  following  sections  propose  a  series  of  tests  for  evaluating  the  consti¬ 
tutive  functions.  At  least  three  modes  of  testing  are  required  to  evaluate  the  functional 
form  of  both  the  rate  equation  and  evolution  equation  for  the  internal  variable.  The 
following  sections  detail  the  three  series  of  tests  used  in  this  investigation:  (a)  isothermal 
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strain  rate  jump  tests  from  the  same  state  to  evaluate  g  in  4.3,  (b)  isothermal  load- 
unload-hold-reload  tests  to  evaluate  r  in  4.4,  and  (c)  isothermal,  constant  true  strain 
rate  tests  to  evaluate  h  in  4.4. 
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4.1  Evaluation  of  the  Rate  Equation:  /(j,0) 


As  just  argued,  if  one  wishes  to  experimentally  probe  the  behavior  of  the  rate  equation 
expressed  by  4.3,  one  must  design  experiments  such  that  material  state  is  held  constant. 
One  can  not  otherwise  separate  the  effect  of  an  evolving  structure  from  the  instantaneous 
response  at  constant  structure.  Very  few  experiments  have  been  performed,  particularly 
within  the  hot  working  regime,  which  have  accomplished  this  separation.  This  is  probably 
due  to  the  fact  that  internal  variable  models  have  not  been  applied  to  hot  working; 
experimentalists  therefore  have  not  designed  their  investigations  with  an  internal  variable 
formulation  in  mind.  It  is  striking,  however,  that  there  are  not  more  constant  structure 
data  a  'ociated  with  the  creep  regime,  for  internal  variables  have  been  applied  to  creep 
constitutive  models  for  many  years. 

Investigators  have  used  different  experimental  techniques  to  characterize  the  rate  de¬ 
pendent  deformation  of  metals  at  assumed  constant  internal  state.  Most  techniques 
involve  an  abrupt  change  in  either  strain  rate  or  stress  and  measurement  of  the  instan¬ 
taneous  or  “short  time”  response  in  the  associated  stress  or  strain  rate.  None  of  the 
techniques,  whether  stress-drop,  relaxation,  or  jump  tests,  is  without  experimental  un¬ 
certainties.  These  uncertainties  range  from  difficulty  in  extracting  the  effect  due  to  a 
particular  mechanism  when  many  possible  mechanisms  may  be  operating,  to  operational 
problems  associated  with  test  machine  compliance,  resolution  of  small  strains,  acquisition 
of  data  from  rapid  transients,  and  computational  uncertainties  in  determining  rate  mea¬ 
sures  from  data  obtained  as  a  function  of  time.  As  a  result  no  one  technique  is  universally 
accepted,  and  the  type  of  test  performed  by  an  investigator  is  generally  a  function  of  the 
type  of  test  equipment  which  is  available  instead  of  the  most  desirable  test  technique. 
All  of  the  tests  mentioned  herein  assume  that  internal  state  does  not  change  significantly 
during  the  period  of  the  transient,  an  assumption  which  is  only  approximately  satisfied 
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in  all  tests. 

The  following  paragraphs  describe  briefly  the  different  tests  used  to  correlate  stress 
and  strain  rate  at  constant  structure,  list  some  of  the  investigators  associated  with  each 
test  type,  and  also  give  some  of  the  uncertainties  associated  with  each  type.  The  sec¬ 
tion  finally  presents  some  arguments  proposing  that  the  strain  rate  jump  test  may  have 
some  clear  advantages  over  stress-drop  and  relaxation  tests  in  determining  isostructural 
properties. 

Mitra  and  McLean  [1967]  used  stress-drop  tests  to  evalute  the  stress  dependence  of  the 
creep  strain  rate  in  aluminum  and  nickel  by  correlating  the  strain  rate  resulting  after  the 
stress  drop  with  the  magnitude  of  the  stress  after  the  drop.  Mitra  and  McLean  assumed 
an  initial  incubation  period  after  which  they  measured  the  strain  rate.  Assuming  that 
the  state  of  the  material  has  not  changed  significantly  during  the  incubation  period,  the 
resulting  stress/strain  rate  curve  may  be  taken  to  be  the  stress  dependence  for  the  given 
initial  internal  state.  Investigators  have  used  the  same  procedure  but  have  measured  the 
“instantaneous  creep  rate”  at  the  new  stress  without  any  assumed  incubation  period,  e.g. 
Sherby  [1957]  and  Robinson  [1969,74,75],  summarized  by  Takeuchi  and  Argon  [1976]  and 
Sherby  [1977].  The  above  mentioned  incubation  period  is  problematic  in  the  interpre¬ 
tation  of  stress  drop  data.  The  incubation  period  is  assumed  necessary  to  let  anelastic 
effects  subside,  but  there  is  no  clear,  and  certainly  no  consensus,  criterion  on  what  that 
period  should  be.  Frequently  investigators  do  not  indicate  what  delay,  if  any,  they  have 
employed  before  determining  the  post-stress-drop,  strain  rate.  Complicating  the  uncer¬ 
tainty  about  incubation  is  the  obvious  possibility  of  structure  evolution  occurring  during 
the  incubation  period,  thereby  voiding  the  constant  structure  assumption  of  the  test. 

Of  the  investigators  who  have  used  stress  drop  tests  to  evaluate  rate  dependence, 
most  of  the  results  mentioned  in  the  papers  listed  above  indicate  a  constant  structure 
power  law  stress  dependence  of  the  strain  rate  of  6  to  8.  Figures  4.2a  and  4.2b  illustrate 
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the  power  law  dependence  of  constant  structure,  stress  drop  data.  This  dependence  in  all 
cases  was  greater  than  the  steady  state  stress  exponent  obtained  for  each  material.  The 
power  law  strain  rate/stress  relationship  for  constant  structure  has  not  been  universally 
confirmed.  Gibeling  and  Nix  [1982],  recently  performed  a  set  of  stress  drop  tests  on  pure 
aluminum  and  found  that  the  isostructural  stress  dependence  of  the  strain  rate  followed 
an  exponential  relationship.  [Figure  4.2c] 

Stress  drop  tests  are  difficult  to  perform  for  hot  working  investigations  due  to  the 
magnitude  of  the  associated  strain  rates.  Evaluation  of  strain  rates  following  a  change 
in  stress  requires  numerical  differentiation  of  a  strain  versus  time  record,  which  in  turn 
requires  sampling  of  strain  over  very  small  time  intervals. 

Hart  [1979]  and  others  have  popularized  the  relaxation  test  as  another  test  of  stress 
dependence  at  constant  structure.  In  this  procedure,  a  specimen  is  deformed  to  a  given 
strain  and  the  test  machine  is  halted.  The  specimen  strain  rate  may  be  determined 
by  subtracting  a  machine  compliance  effect  from  the  load  versus  time  curve  obtained 
during  relaxation.  Hart  estimates  that  the  inelastic  strain  accumulated  during  a  typical 
relaxation  test  is  about  10-4.  From  this  Hart  states  that  such  a  small  strain  increment 
produces  negligible  strain  hardening,  so  the  test  can  be  considered  to  be  at  a  constant 
state.  This  assumption  is  valid  only  as  long  as  static  recovery  effects  can  be  neglected, 
which  is  certainly  not  possible  in  the  range  of  homologous  temperatures  associated  with 
hot  working.  Other  investigators  using  this  technique  include  Alexopoulos  [1982]  and 
C.  Li  [1981].  The  relaxation  test  is  problematic  in  that  some  investigators  have  pro¬ 
posed  using  it  to  characterize  recovery  mechanisms,  an  application  directly  contrary  to 
Hart’s  assumption  of  zero  state  change.  Other  investigators  have  suggested  that  the 
small  strains  associated  with  relaxation  tests  do  not  accurately  probe  the  large  strain 
constitutive  response  of  a  material  since  little  dislocation  motion  is  required  to  unload 
the  test  specimen  [Rhode  et  al.  1973]. 
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Another  procedure  used  to  measure  the  constant  structure,  strain  rate  or  stress  de¬ 
pendence  is  to  suddenly  change  the  strain  rate.  The  instantaneous  change  in  stress 
accompanying  a  change  in  strain  rate  may  be  a  more  faithful  representation  of  the  stress 
dependence  than  that  which  results  from  the  stress  drop  test,  since  there  is  no  associated 
incubation  time.  Klepaczko  and  Duffy  [1977,  1982,  1986]  have  probably  performed  the 
most  comprehensive  set  of  strain  rate  jump  tests,  having  documented  the  behavior  of  both 
FCC  and  BCC  materials  over  a  wide  range  of  temperatures  and  strain  rates,  although 
most  of  their  test  temperatures  were  below  one  half  the  homologous  temperature. 

It  should  be  noted  that  very  few  of  the  investigators  mentioned  above  have  employed 
any  of  the  “jump”  tests  described  above  in  a  manner  appropriate  to  decoupling  state 
from  instantaneous  rate  dependence.  Such  decoupling  requires  jumps  from  the  same 
initial  state.  In  the  context  of  the  single,  isotropic  internal  variable  mode)  under  con¬ 
sideration  here,  this  means  from  the  same  value  of  temperature,  strain  rate,  and  stress. 
Given  that  the  real  material  may  (and  probably  does)  exhibit  history  dependence  beyond 
that  representable  by  a  single  internal  variable,  the  jumps  should  occur  after  identical 
deformation  histories  to  assure  as  identical  a  state  as  possible.  In  a  case  where  a  steady 
state  is  reached,  it  is  useful  to  impose  the  sudden  change  from  steady  state  conditions 
where  the  jump  state  is  more  reliably  reproduced.  A  set  of  such  jump  tests  is  indicated 
in  Figure  4.3.  Figure  4.4  illustrates  a  set  of  jump  tests  on  the  schematic  introduced  in 
Figure  4.1,  where  the  test  begins  with  the  value  of  the  internal  variable  so  and  eventually 
saturates  at  s*. 

Interpretation  of  jump  test  data  is  complicated  by  different  back-extrapolation  meth¬ 
ods  used  to  determine  the  instantaneous  stress  dependence,  especially  when  there  is  an 
overshoot  and  subsequent  drop  in  stress.  Kocks,  Argon,  and  Ashby  [1975]  review  the 
different  methods  and  conclude  provisionally  that  the  method  of  extrapolation  is  not  of 
critical  importance.  The  jump  test  data  obtained  in  this  investigation  does  not  demon- 
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strate  any  overshoot.  The  data  show  a  linear  segment  which  then  increases  monotonically 
with  increasing  strain.  One  may  therefore  select  a  simple  and  straightforward  technique 
of  defining  the  instantaneous  stress  response  as  that  resulting  from  a  0.2  percent  strain 
offset,  illustrated  in  a  schematic  strain  rate  jump  test  in  Figure  4.9,  where  a/  is  used 
as  the  constant  structure  stress  response.  This  necessarily  assumes  that  one,  there  is 
no  change  in  state  during  this  strain  increment,  and  two,  that  the  transition  to  the  new 
strain  rate  has  been  accomplished  over  the  small  strain  increment.  The  0.2%  strain  offset 
criteria  was  chosen  over  some  back  extrapolation  procedure  since  the  presence  of  exten¬ 
sive  hardening  made  extrapolation  from  some  point  further  from  the  jump  point  difficult 
to  justify. 

Very  few  investigators  have  attempted  to  examine  the  constant  structure  temperature 
dependence  of  the  flow  behavior  of  metals  within  the  regime  of  temperatures  and  strain 
rates  associated  with  creep  or  hot  working.  Direct  measurement  of  this  dependence, 
and  thus  the  activation  energies  associated  with  the  dislocation  glide  mechanisms,  re¬ 
quires  temperature  jump  tests.  The  activation  energies  obtained  from  steady  state  data 
represent  the  combined  contribution  of  both  glide  kinetics  and  the  thermal  processes 
controlling  the  evolution  of  microstructure.  We  did  not  perform  anj  temperature  jump 
testing  in  this  investigation  because  the  time  constant  associated  with  changing  temper¬ 
atures  in  the  test  equipment  was  significantly  larger  than  the  recovery  rates  associated 
with  the  model  material. 

4.1.1  Strain  Rate  Jump  Tests 

As  discussed  in  the  previous  section,  jump  tests  measure  material  response  at  a  given  con¬ 
stant  structure,  thereby  characterizing  the  rate  equation  associated  with  our  constitutive 
model.  A  series  of  strain  rate  jump  tests  were  performed  on  the  test  system  described  in 
Chapter  3.  All  strain  rate  jumps  were  programmed  to  occur  at  the  same  value  of  strain 
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for  a  given  temperature,  after  the  stress  variation  with  strain  had  reached  an  apparent 
constant  value  at  the  initial  strain  rate.  Jumping  from  such  apparently  steady  state 
conditions  minimized  the  variation  in  initial  state  among  tests  at  the  same  temperature 
and  initial  strain  rate.  Figures  4.5  to  4.8  show  the  strain  rate  jump  test  results  at  700, 
800,  900  and  1000  degrees  Celsius.  All  jumps  in  strain  rate  were  from  the  same  initial 
strain  rate  of  0.2  x  10-4  per  second.  The  maximum  final  strain  rate  was  selected  to  be 
1  per  second  due  to  the  limit  of  the  servohydraulic  dynamic  response.  This  strain  rate 
also  corresponded  to  the  rate  at  which  noticeable  adiabatic  heating  occurred,  producing 
an  effective  limit  to  the  isothermal,  constant  true  strain  rate  test  conditions. 

The  displacement  command  signal  for  the  strain  rate  jump  tests  was  provided  by 
the  analog  function  generator  described  in  Chapter  3.  Data  acquisition  was  performed 
on  a  IBM  PC/XT  equipped  with  a  Metrabyte  Dash-16  data  acquisition  board.  Two 
programs  were  used  on  the  PC/XT,  both  of  which  are  listed  and  described  in  Appendix 
A  .  Both  the  load  cell  and  LVDT  output  signals  were  also  conditioned  by  a  60  Hertz 
cut-off,  low  pass  filter  to  reduce  the  electrical  noise  generated  by  the  vacuum  furnace 
heating  elements.  The  filter  frequency  was  selected  to  permit  as  much  noise  reduction  os 
possible  without  distorting  transducer  output  signals. 

Once  stored  on  the  PC/XT,  the  load/displacement  data  were  transferred  to  a  Data 
General  MV4000  where  they  were  converted  to  true  stress/strain  data.  The  stress  and 
strain  values  were  then  independently  smoothed,  using  sample  number  as  the  independent 
variable,  and  then  the  plastic  strain  and  strain  rate  were  calculated  at  each  sample  point. 
The  programs  and  procedures  used  for  this  analysis  is  described  in  Appendix  A. 

The  stress  before  and  immediately  after  the  jump  was  then  measured.  The  stress 
after  the  jump  was  selected  to  be  the  stress  corresponding  to  a  0.002  offset  strain  from 
the  jump  strain,  as  illustrated  in  Figure  4.9.  It  was  frequently  difficult  to  determine  this 
value  precisely  due  to  the  absence  of  a  clearly  linear,  elastic,  initial  portion  of  the  jump 
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in  stress.  The  initial  portion  of  the  curve,  although  linear,  decreased  in  slope  between 
tests  as  the  final  strain  rate  was  increased.  In  any  event,  the  linear  section  of  the  curve 
was  assumed  to  be  elastic,  and  the  offset  was  taken  from  this  section.  The  uncertainty  in 
this  linear  region  did  not  significantly  change  the  value  of  this  offset  stress.  Examination 
of  the  strain  rate  data  indicated  that  the  final  strain  rate  was  achieved  at  approximately 
the  same  time  as  the  0.002  strain  offset  stress  was  reached.  Given  the  assumption  of 
no  change  in  state  up  to  that  point,  the  offset  stress  could  then  be  assumed  to  indicate 
the  stress  sensitivity  of  the  material  at  a  given  internal  state.  Figure  4.10  shows  the 
constant  state  lines  constructed  from  this  stress/strain  rate  data.  Table  4.1  lists  the  data 
associated  with  these  tests  which  were  used  to  produce  Figure  4.10. 

Three  rate  equations  were  considered  for  correlation  with  the  constant  structure  data. 
One,  a  simple  model  representing  a  stress  and  structure  dependent  thermal  activation: 

u  Ae*P  Hw  (' -  f)]  •  4-Ll 

Here,  A  and  Q  are  material  parameters.  Two,  a  modification  of  4.1.1  to  reflect  the 
phenomena  of  power  law  breakdown: 

j=/l(f),exp[--|(l-f)],  4.1.2 

where  an  additional  parameter  q  has  been  added.  Three,  an  alternative  to  4.1.1  which 
separates  temperature  and  stress  dependence  into  a  simple  Arrhenius  term  involving 
a  constant  activation  energy  and  a  function  depending  only  on  stress  and  the  internal 
variable: 

Four  material  parameters  are  required:  A,  Q,£,  and  m.  R  is  the  gas  constant  in  these 
equations  (8.314  x  10-3  k Joules/mole).  Equation  4.1.1  reflects  a  simplification  of  a  form 
proposed  by  Kocks,  Argon,  and  Ashby  [1975]  for  the  jerky  glide  of  dislocations.  Here, 
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s  represents  a  generalized  obstacle  resistance  to  dislocation  motion.  Equation  4.1.2, 
proposed  by  Lee  and  Zaverl  [1978],  follows  the  same  motivation,  but  assumes  a  priori  that 
the  rate  equation  should  reflect  the  transition  from  power  law  to  exponential  behavior, 
and  therefore  includes  a  stress-dependent  pre-exponential  to  accommodate  the  power 
law.  Equation  4.1.3  makes  the  same  a  priori  assumption  but  uses  a  modification  of  the 
following  hyperbolic  sine  form  first  proposed  by  Garofalo  [1963]  to  model  steady  state 
behavior  into  the  power  law  breakdown  regime: 

=  Ass  exp  (-■ [sinh(c>f<7ss)]1/m  ,  4.1.4 

where  the  “ss”  subscript  denotes  quantities  relevant  only  to  steady  state  conditions.  This 
form  accomodates  power  law  breakdown  since  for  low  values  of  (acr5J)  it  approximates: 

6SS  =  Ass  exp  (~^)  4.1.5 

and  for  large  values  of  the  same  argument  it  approaches: 


There  appears  to  be  no  widely  accepted  mechanistic  derivation  of  the  hyperbolic  sine 
form  of  the  rate  equation  represent  by  4.1.4.  Gittus  [1976]  has  proposed  an  explanation 
of  power  law  breakdown  founded  on  the  dominance  of  forward  and  backward  vacancy 
formation  rates.  A  hyperbolic  sine  form  results  naturally  from  this  assumption  of  forward 
and  backward  mechanical  activation.  Many  other  investigators,  [Sherby,  1968;  Frost  and 
Ashby,  1982]  associate  power  law  breakdown  with  a  transition  from  climb- dominated 
deformation  to  dislocation  glide  domination.  In  any  event,  it  is  well-known  that  the 
hyperbolic  sine  form  provides  a  convenient  phenomenological  representation  of  power 
law  breakdown. 

We  settled  upon  the  hyperbolic  sine  form,  4.1.3,  for  the  rate  equation  for  the  following 
reasons: 
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1.  The  steady  state  stress/strain  rate  relationship  for  pure  aluminum  and  Fe  -  2% 
Si  suggests  a  constant  activation  energy.  Figure  4.11  illustrates  the  steady  state 
stress/strain  rate  relationship  for  aluminum  over  a  wide  range  of  strain  rates  and 
temperatures,  where  the  steady  state  strain  rate  is  normalized  by  a  constant  ac¬ 
tivation  energy  [Jonas,  1969].  In  light  of  the  success  of  this  normalization,  we 
found  it  reasonable  to  extend  this  relationship  to  the  rate  equation  and  assume  a 
relationship  of  the  form  of  4.1.4. 

The  jump  test  data  for  the  Fe  -  2%  Si  further  suggests  a  functional  dependence  of 
the  form  of  4.1.3.  Figure  4.12  plots  the  constant  structure  jump  test  data  where 
the  strain  rate  has  been  normalized  by  a  constant  activation  energy,  taken  in  this 
case  to  be  that  for  self  diffusion  in  alpha  iron.  It  is  obvious  that  a  power  law 
does  not  represent  the  dependence  between  stress  and  the  normalized  strain  rate. 
With  the  exception  of  the  700  degree  Celsius  data,  each  constant  structure  curve 
is  concave  upward  in  a  manner  commensurate  with  power  law  breakdown.  The 
700  degree  data  is  suspect  for  reasons  associated  with  the  dynamic  response  of  the 
servohydraulic;  it  is  included  however  for  both  for  completeness  and  because  we  are 
not  certain  that  the  data  is  actually  incorrect.  It  would  be  interesting  to  perform 
constant  structure  tests  at  higher  strain  rates  to  investigate  this  behavior. 

Although  the  constant  structure  data  for  the  Fe  -  2%  Si  may  be  represented  by  an 
exponential,  constant  structure  stress  drop  tests  performed  on  other  metals  suggest 
a  power  law  asymptote  in  the  low  strain  rate  creep  regime  [Sherby  et  al.  1977].  For 
this  reason  we  selected  equation  4.1.3  which  may  accomodate  this  transition  from 
an  exponential  to  a  power  law. 

2.  Regardless  of  the  form  of  the  rate  equation,  we  insist  that  the  steady  state  value 
of  the  internal  variable  (sJ3)  increase  monotonically  with  steady  state  stress  (crS3). 
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This  not  only  makes  sense  physically  but  also  corresponds  to  measurements  of  room 
temperature  yield  stress  as  a  function  of  steady  state  stress  [Young  and  Sherby, 
1973].  One  may  check  whether  the  three  proposed  rate  equations  accomodate  this 
requirement  by  assuming  steady  state  conditions  and  equating  each  rate  equation 
to  known  empirical  relations  for  steady  state  behavior.  By  eliminating  strain  rate 
as  the  common  variable  we  may  determine  whether  the  resulting  functional  de¬ 
pendence  between  the  steady  state  values  of  stress  and  the  internal  variable  is 
reasonable. 

Following  this  procedure,  equation  4.1.1  requires  an  unreasonable  relationship  be¬ 
tween  steady  state  stress  and  steady  state  value  of  the  internal  variable  in  order 
to  reproduce  power  law  breakdown.  The  steady  state  strain  rate/stress  relation¬ 
ship  for  many  materials  including  most  metals  in  the  hot  working  regime  may  be 
represented  by  equation  4.1.4: 

=  Ass  exp  [sinh(a<7ss)]1/m  ,  4.1.4 

where  the  “ss”  subscript  denotes  quantities  relevant  only  to  steady  state  conditions. 
If  we  assume  steady  state  conditions  and  then  apply  4.1.4  for  the  steady  state  strain 
rate  in  4.1.1,  we  obtain: 

s“  =  m17"  [ln {t  exp (“iw22) •  4.1.7 

Wong  and  Jonas  [1968]  have  fit  4.1.4  to  pure  aluminum  within  a  wide  range  of 
temperatures,  strain  rates,  and  stresses  to  find  the  following  constants: 

Q,s  =  156.0  kJ/mole, 
a  =  .0446  MPa-1, 

1/m  =4.70, 
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Aas  =2.34  x  lO10  sec"1. 


If  we  assume  additionally  the  following  values  for  equation  4.1.1: 

A  =  1  x  lO10  sec-1,  and 
Q  =  142.0  kJ/mole,  (self-diffusion), 

then  substitution  of  these  values  into  4.1.7  yields,  for  a  representative  temperature 
of  200  degrees  Celsius: 

s4J  =  36.0<7j4  [in  |0.067[sinh(.0446cr„)]4’7||  .  4.1.8 

The  steady  state  stress/internal  variable  relationc.iip  expressed  by  4.1.8  is  shown  in 
Figure  4.13.  The  negative  values  of  state  and  the  singularity  indicates  that  equation 
4.1.1  is  inadequate  to  simulate  hot  working. 

3.  Generally  accepted  notions  of  the  stress  dependence  of  the  pre-exponential  create 
a  situation  for  equation  4.1.2  similar  to  that  described  in  item  2  above  for  equation 
4.1.1.  Stress  dependence  of  the  pre-exponential  is  generally  derived  to  be  propor¬ 
tional  to  stress  raised  to  a  low  power,  either  one  or  two  [Argon,  1975].  If  we  assume 
two  to  be  an  upper  bound  and  then  follow  a  procedure  identical  to  that  followed  in 
reason  two  above,  we  obtain  the  following  expression  correlating  steady  state  stress 
and  steady  state  deformation  resistance: 

'4"exp  (-%)  lsi"h(aa“)),/m  =  A  ©”eXp  [~W  ('  -  Z)\  ’  4-U 

where  the  variables  have  the  same  interpretations  as  in  item  1.  Rearranging  and 
using  the  same  values  for  aluminum  used  before  yields: 

0.067[sinh(.0446a„)]4’7  =  (— )  exp  36.0—  .  4.1.10 

\Sjj/  sJS . 

The  values  of  $„  corresponding  to  variation  in  crss  may  be  determined  numerically 
from  4.1.10.  Figure  4.14  plots  the  relationship  resulting  from  such  a  determination. 
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Once  again  the  variation  in  the  value  of  the  internal  variable  at  steady  state  with 
respect  to  steady  state  stress  is  not  what  we  should  reasonably  expect. 


4.  The  fact  that  the  ratio  of  equivalent  stress/internal  variable  appears  twice  in  equa¬ 
tion  4.1.2  complicates  the  determination  of  material  parameters  significantly  when 
fitting  the  rate  equation  to  actual  material  data.  As  will  be  seen,  much  of  the 
following  analysis  exploits  a  desirable  feature  of  the  rate  equation  to  be  inverted  to 
produce  an  analytic  relationship  of  the  form: 


a  =  cs, 

where  c  depends  only  on  temperature  and  strain  rate.  Equation  4.1.2  does  not 
permit  such  an  inversion. 

The  material  parameters  associated  with  equation  4.1.3  may  be  evaluated  directly  from 
the  jump  test  data.  The  procedures  and  numerical  routines  used  to  obtain  these  pa¬ 
rameters  are  described  in  Chapter  5.  Figure  4.15  illustrates  the  correspondence  between 
the  experimental  data  and  the  rate  equation  with  the  following  values  for  the  material 
parameters: 

Q  =  247.5  kJ/mole 
A  =  1.26  x  108  sec-1 
m  =  .22793 

The  value  of  Q  is  within  the  range  of  activation  energies  obtained  by  other  investigators 
for  self-diffusion  in  alpha  iron:  239  to  251  kj/mole  [Frost  and  Ashby,  1982).  It  is  lower 
than  the  value  of  333.6  kj/mole  obtained  by  Uviia  and  Jonas  [1968]  by  fitting  equation 
4.1.4  to  a  set  of  steady  state  stress  and  stress  rate  data  on  a  Fe  -  3%  Si  alloy.  There  is 
no  reason  to  expect  the  activation  energies  obtained  in  this  investigation  and  by  Uvira 
and  Jonas  to  be  similar.  The  value  r’  ^ined  through  the  jump  tests  is  a  consequence  of 
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constant  structure  rate  dependence  while  that  obtained  from  steady  state  data  provides 
a  weighted  average  of  both  this  dependence  and  the  thermally  activated  processes  con¬ 
trolling  the  evolution  of  state.  The  correlation  between  the  data  and  the  fit  rate  equation 
is  excellent.  We  therefore  adopted  equation  4.1.3  as  the  rate  equation.  The  parameter  £ 
in  the  rate  equation  is  indeterminate.  We  chose  to  include  £  in  order  to  restrict  values  of 
the  internal  variable  s  to  be  at  all  times  greater  than  a.  This  corresponds  to  requiring 
the  proportionality  constant  c  to  be  less  than  one,  where  a  =  cs  for  isothermal,  constant 
true  strain  rate  conditions.  For  4.1.3,  c  is: 


c  = 


?Sinh_1 


4.1.12 


For  the  rest  of  this  investigation  we  assume  a  value  of  f  =  5.0. 

To  recapitulate,  the  rate  equation  which  will  be  used  throughout  the  rest  of  this 
investigation  is: 


e  =  A  exp  (--|)  [sinh  («f) 


1/m 


4.1.3 


where  A,  Q ,  and  m  are  material  parameters  which  should  be  determined  from  constant 
structure,  jump  test  data. 


4.2  Evaluation  of  the  Static  Recovery  Function:  r(0,  s) 

Several  means  are  available  for  decoupling  the  dynamic  and  static  terms  of  the  evolution 
equation  for  s  represented  by  equation  4.4,  restated  here: 

s  =  h(a,s,e)‘f  -r(s,0).  4.4 

Certainly,  one  may  assume  that  one  term  predominates  within  a.  particular  tempera- 
turc/strain  rate  regime.  We  expect  the  effect  of  the  static  recovery  function  r  to  diminish 
as  the  strain  rate  increases  or  as  the  temperature  decreases.  One  may  therefore  assign 
a  regime  where  static  recovery  may  be  neglected  and  the  strain  hardening  data  used  to 
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characterize  only  the  dynamic  term  of  the  evolution  equation.  We  may  similarly  perform 
tests  where  the  hardening  function  h  may  be  considered  neligible.  As  will  be  shown  be¬ 
low,  this  may  be  accomplished  either  via  some  extrapolation  technique  or  by  static  hold 
tests  which  enforce  a  zero  plastic  strain  rate. 

Experimental  determination  of  a  static  recovery  function  of  the  form 

r  =  r(6,s )  4.2.1 

is  particularly  arduous  since  recovery  testing  typically  yields  only  one  data  point,  i.e.  a 
“recovered”  state,  per  test.  Characterization  of  a  recovery  function  therefore  requires  nu¬ 
merous  tests  to  completely  capture  the  effect  of  each  of  the  relevant  parameters.  Different 
investigators  have  proposed  different  experiments  to  examine  static  recovery  phenomena. 
Most  assume  a  Bailey-Orowan  form  for  the  combined  effect  of  hardening  and  recovery, 
and  then  propose  a  procedure  whereby  the  hardening  rate  may  be  considered  to  be  negli¬ 
gible  relative  to  the  recovery  term.  We  have  not  found  any  investigator  who  has  evaluated 
recovery  data  directly  in  order  to  determine  an  internal  variable  based,  static  recovery 
function. 

Perhaps  the  first  comprehensive  phenomenological  investigation  of  hardening  and 
recovery  was  performed  by  Mitra  and  McLean  [1966],  who,  following  a  suggestion  of 
Cottrell  and  Aytekin  [1950],  used  stress  drop  experiments  to  evaluate  a  recovery  function 
based  on  stress.  They  assumed  that  the  recovery  function  could  be  represented  as: 


They  then  performed  a  series  of  stress  drop  tests  from  the  same  initial  stress,  tempera¬ 
ture,  and  strain,  and  then  associated  the  drop  in  stress  (A a)  with  the  time  increment 
(A t)  required  to  reach  a  new  steady  state  strain  rate.  Figure  4.16  illustrates  one  such 
stress  drop  test.  They  then  plotted  the  stress  reductions  versus  the  associated  recovery 
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time.  The  instantaneous  recovery  rate  was  assumed  to  be  slope  of  the  resulting  curve 
extrapolated  back  to  a  zero  stress  reduction: 

r  =  lim  (Aa/ At),  4.2.2 

which  is  illustrated  in  Figure  4.17. 

Kocks  [1975]  proposed  two  methods  of  evaluating  r.  Method  one  consisted  of  imposing 
strain  rate  jumps  from  identical  steady  state  conditions  and  measuring  the  associated 
instantaneous  rate  of  change  in  stress.  Back  extrapolation  of  the  resulting  data  to  zero 
strain  rate  as  illustrated  in  Figure  4.18  would  then  provide  the  recovery  rate.  Method  two 
consisted  of  measuring  the  rate  of  change  of  stress  immediately  following  a  large  drop  in 
strain  rate  from  steady  state  conditions  (where  the  rate  of  recovery  supposedly  equals  that 
of  hardening).  If  the  drop  in  strain  rate  is  at  least  one  order  of  magnitude,  then  Kocks 
hypothesized  that  since  the  rate  of  hardening  is  assumed  in  a  Bailey-Orowan  formulation 
to  be  proportional  to  the  strain  rate,  then  the  hardening  term  would  be  negligible  relative 
to  the  static  softening  term,  which  would  be  unaffected  by  the  reduction  in  strain  rate. 
The  rate  of  change  of  stress  would  therefore  be  due  purely  to  static  recovery  in  the 
absence  of  significant  hardening  or  dynamic  recovery.  Both  of  these  methods  possess  the 
operational  disadvantage  of  attempting  to  measure  the  rate  of  change  of  stress  associated 
with  a  very  short  transient.  No  investigations  have  been  found  which  used  either  of  these 
two  methods  to  evaluate  the  recovery  function  r. 

An  additional  test  of  static  softening  mechanisms  which  has  been  employed  exten¬ 
sively  in  investigations  of  metal  working  behavior  is  the  load-unload-hold-reload  test 
[McQueen.  1985;  Luton,  et  al.,  1980;  Petkovic,  1979].  Employed  primarily  to  simulate 
multistage  hot  working,  specimens  are  deformed,  usually  isothermally,  either  in  compres¬ 
sion  or  torsion  to  a  given  strain,  unloaded,  held  for  varying  time  periods,  then  reloaded. 
The  load-unload-reload  cycle  may  be  repeated  many  times  sequentially  on  the  same 
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specimen  to  simulate  complicated,  multistage  hot  working  processes.  The  change  in 
state  resulting  from  a  particular  hold  period  may  be  characterized  by  the  yield  stress  ob¬ 
tained  on  reloading  the  specimen.  The  test  possesses  the  benefit  of  avoiding  the  multiple 
uncertainties  (and  controversies)  associated  with  both  measuring  and  interpreting  the 
transients  associated  with  changing  the  stress  or  strain  rate.  The  test,  however,  suffers 
the  same  shortcoming  as  the  jump  test,  since  some  semi- arbitrary  criteria  must  be  used 
to  measure  the  yield  stress  upon  reloading.  The  same  arguments  can  be  made  as  in  the 
case  of  the  jump  test  that  the  variation  in  yield  stress  is  relatively  small  given  sufficient 
resolution  on  the  reloading  transient. 

Rather  than  fixing  a  final  form  for  the  recovery  function  r  and  then  using  recovery 
data  to  determine  the  parameters  of  that  function,  the  following  technique  is  proposed  to 
characterize  the  static  recovery  function  assuming  the  following  intermediate  structure 
for  the  evolution  equation  for  the  internal  variable: 

s  =  h{a,0,s)f  -r(s)exp  (-■ .  4.2.4 

Here,  we  assume  that  the  dependence  of  static  recovery  on  the  temperature  and  the 
internal  variable  may  be  expressed  through  the  product  of  an  Arrhenius  term  with  a 
constant  activation  energy,  Qr,  and  a  function  depending  only  on  s.  The  activation 
energy  Qr  is  expected  to  be  different  from  that  associated  with  the  rate  equation,  Q. 
Since  static  recovery  is  generally  assumed  to  be  the  result  of  dislocation  climb  [Hirth, 
1982]  Qr  should  be  close  to  the  activation  energy  for  the  self-diffusion  of  vacancies.  We 
therefore  assume  Qr  =  300  kJ/mole,  which  is  an  average  value  for  self-diffusion  in  alpha 
iron  [Frost  and  Ashby.  1982].  Chapter  5  describes  a  series  of  tests  which  validates  this 
assumption.  The  expectation  is  that  we  may  design  an  appropriate  technique  so  that  we 
only  assume  a  priori  equations  4.3  and  4.2.4,  and  then  let  the  experimental  data  suggest 
a  form  for  r(s). 
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The  proposed  procedure  exploits  the  consequence  of  the  above  equations  that  the 
internal  variable  5  is  at  all  times  proportional  to  the  stress  a  during  an  isothermal, 
constant  true  strain  rate  test.  That  is, 


a  —  cs,  4.2.5 

where  c  is  a  function  of  the  rate  equation  which  includes  strain  rate  and  temperature 
dependence  (equation  4.1.12).  This  proportionality  is  central  to  the  evaluation  of  both  the 
static  recovery  function  r  and  the  hardening  function  h,  for  it  allows  us  to  use  a  measure 
of  stress  to  determine  the  internal  variable  s  through  the  constant  c.  We  therefore  do 
not  have  to  measure  s  directly;  we  need  only  measure  stress  and  assure  that  isothermal, 
constant  true  strain  rate  conditions  exist.  This  assumption  may  be  made  as  long  as  the 
functional  dependence  in  the  rate  equation  between  the  equivalent  stress  and  the  internal 
variable  is  one-to-one  for  a  given  temperature  and  strain  rate. 

If  we  now  perform  load-unload-hold-reload  experiments  in  isothermal,  constant  true 
strain  rate  conditions,  the  difference  between  the  stress  just  before  unloading  and  the  yield 
stress  upon  reloading  represents  the  change  in  state  due  only  to  static  recovery.  We  can 
then  vary  temperature  (0),  hold  time  (At),  and  stress  before  unloading  (cr,),  maintaining 
a  constant  value  of  c,  and  thus  determine  the  temperature  and  state  dependence  of  the 
recovery  function  by  measuring  the  value  of  stress  upon  reloading  (cr/). 

The  procedure  proposed  assumes  that  we  have  already  determined  the  functional 
dependence  on  temperature,  which  as  stated,  will  be  represented  by  an  Arrhenius  ex¬ 
pression  with  an  activation  energy  equal  to  that  of  self-diffusion.  The  data  required  is  a 
set  of  load-unload-hold-reload  tests,  where  the  stress  before  unloading,  the  temperature, 
and  strain  rate  are  held  constant.  Figure  4.19  illustrates  a  series  of  such  load-unload- 
hold-reload  tests.  The  only  variable  changed  between  tests  is  the  hold  time,  At.  The 
dependent  variable  is  the  stress  obtained  upon  reloading  the  specimen,  Of.  The  proce- 
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dure  effectively  measures  the  integrated  effect  of  recovery,  then  differentiates  it  to  obtain 
the  recovery  function.  This  is  in  contrast  to  the  methods  described  at  the  beginning  of 
this  section  which  attempt  to  measure  recovery  function  directly. 

In  the  absence  of  hardening  or  dynamic  recovery  (eP  =  0),  the  evolution  equation 
becomes: 


s  = 


4.2.6 


Integration  yields 


exp 


Qr\  ds_ 
'(s) 


9l)  r 

ROJ  4 


=  -  At. 


4.2.7 


However,  since  we  only  know  a  and  not  s,  we  may  express  this  as: 


(Qr\  [°^c  d(cr / c) 

exp  [w  L  WR) 


-At. 


4.2.8 


From  load-unload-hold-reload  tests  we  may  determine  and  At.  Differentiating  the 

integral  on  the  left  hand  side  of  4.2.8  with  respect  to  07/c  yields: 

d  n'c  dja/c)  =  1  A  2  9 

d(<rj/c)  Joi/c  r(cr/c)  f(<r//c)’ 


or, 


6XP  (ro)  r(<7f/c)  d(of/c)  At 


4.2.10 


'i<7j/c)  d(<Jf/c) 

If  we  now  know  (or  assume)  Qr ,  we  may  use  4.2.8  to  plot  the  hold  period  —At  versus  cr^/c, 
the  derivative  of  which,  according  to  4.2.10,  is  the  reciprocal  of  the  recovery  function: 


r{af/c)  =  exp  4.2.11 

'  d(oj/c) 

Figure  4.20  shows  a  schematic  plot  of  hold  period  (At)  versus  the  stress  measured  upon 
reloading  (^"y),  which  may  be  obtained  from  a  series  of  tests  like  those  illustrated  m 
Figure  4.19.  Figure  4.21  illustrates  the  result  of  treating  the  data  in  Figure  4.20  in  the 
manner  proposed  by  equation  4.2.11.  The  data  plotted  in  Figure  4.21  then  represents 
the  recovery  function  r(s,0). 
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This  procedure  may  be  applied  to  the  static  recovery  function  proposed  by  both 
Prinz  and  Argon  [1984]  and  Nix  and  Gibeling  [1985].  Slightly  generalized,  the  form  of 
the  function  is: 


r  = 


4.2.12 


Here,  p  is  the  shear  modulus  and  B  and  p  and  material  constants.  There  is  an  additional 
temperature  dependence  of  the  pre-exponential  in  the  Prinz  and  Argon  and  the  Nix  and 
Gibeling  models  which  has  been  neglected  in  4.2.12.  Given  load-unload-hold-reload  data 
as  described  above,  4.2.12  applied  to  4.2.10  yields: 


-dAt  /<3r\  ^  (a*\  P 
d{<r,/c)  ~  GXP  \R0J  fiB  V  c  ) 


4.2.13 


The  above  procedure  requires  enough  tests  at  a  given  temperature  to  be  able  to 
determine  the  derivative  on  the  right  hand  side  of  4.2.10.  It  also  presupposes  knowledge 
of  the  value  of  c  and  Qr.  Once  again,  c  is  the  proportionality  constant  defined  by  equation 
4.1.12  and  obtained  from  the  strain  rate  equation.  QT  may  be  estimated  as  the  activation 
energy  for  self- diffusion  for  the  material  of  interest.  It  may  also  be  determined  by  a  series 
of  tests  where  through  trial  and  error  the  same  values  of  a,  and  oy  are  obtained  for  two 
tests  at  different  temperatures  but  the  same  value  of  c.  This  can  be  accomplished  only 
by  allowing  a  longer  recovery  time  for  the  lower  temperature  test.  If  these  conditions  are 
met,  then  for  the  two  tests  at  two  temperatures: 


exp  (■  r£>  A(i = exp  (~m)  a<2’ 


4.2.14 


and  Qr  may  be  determined. 


4.2.1  Recovery  Tests 

A  set  of  load-unload-hold-reload  tests  were  performed  on  the  silicon  iron  according  to 
the  paradigm  described  above.  Before  performing  these  tests,  however,  we  performed  a 


series  of  tests  to  determine  the  temperature  and  strain  at  which  static  recrystallization 
initiated.  We  then  selected  a  temperature  and  strain  below  these  values  to  assure  that 
the  softening  measured  due  to  the  hold  period  would  be  due  solely  to  static  recovery. 

The  initiation  of  static  recrystallization  was  determined  phenomenologically  according 
to  a  procedure  described  by  McQueen  [McQueen,  1982;  Ryan  et  ah,  1983].  A  series  of 
tests  were  performed  at  increments  of  100  degrees  Celsius  where  a  compression  specimen 
was  deformed  at  a  constant  strain  rate  and  temperature  for  a  given  strain  increment, 
unloaded,  and  held  for  a  specified  time.  This  was  repeated  for  the  same  strain  increment 
until  a  true  strain  of  approximately  —1.2.  The  fractional  softening  for  each  hold  period 
was  then  plotted  as  a  function  of  strain  for  each  test.  Fractional  softening  is  defined  as: 

FSn  =  ~ <Tr(n+1)  4.2.15 

O’un  CTyO 

where 

crun  =  stress  before  unloading  nth  strain  increment, 

0>(n+i)  =  yield  stress  upon  reloading  from  nth  hold  period, 

CTyo  =  original  yield  stress. 

McQueen  found  that  during  a  given  test,  fractional  softening  remained  relatively  con¬ 
stant  if  only  static  recovery  was  present.  Static  recrystallization,  however,  past  a  certain 
initiation  strain  caused  the  fractional  softening  to  increase. 

The  results  of  our  multiple  hold  tests  are  included  in  Figures  4.22  to  4.25.  The 
fractional  softening  versus  strain  is  plotted  in  Figure  4.26.  Notice  that  the  fractional 
softening  begins  to  increase  at  1000  degrees  Celsius  at  a  true  strain  of  0.6.  We  therefore 
chose  our  recovery  test  conditions  to  be  at  900  degrees  Celsius,  at  a  strain  rate  less  than 
that  used  .  r  the  multiple  hold  tests,  (—0.02  per  second  instead  of  —0.1  per  second), 
and  at  a  strain  of  -0.3.  Examination  of  the  etched  grains  of  specimens  tested  at  this 
conditions  did  not  show  any  regions  of  recrystallized  grains. 
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The  results  of  the  recovery  tests  for  these  conditions  are  presented  in  Figure  4.27.  One 
obvious  result  of  these  tests  is  that  the  stress  measured  upon  reloading  is  not  significantly 
different  from  that  measured  just  before  unloading,  even  for  relatively  long  hold  periods. 
We  should  therefore  expect  the  rate  of  static  recovery  to  be  commensurately  low. 

Another  result  from  these  recovery  tests  is  that  when  the  “yield  stress”  is  measured 
upon  reloading,  the  greatest  incremental  decrease  in  stress  occurs  during  the  shortest 
hold  period.  In  other  words,  the  rate  of  static  recovery  decreases  dramatically  as  the 
hold  period  increases.  This  is  not  necessarily  unexpected,  for  we  may  hypothesize  short 
range  rearrangements  due  to  the  polarization  of  a  deformation-induced  microstructure 
or  anelastic  effects  which  produce  a  rapid  apparent  initial  rate  of  recovery.  This  result 
was  not  obtained  by  Petkovic  [1979]  who  performed  similar  load-unload-hold-reload  tests 
on  polycrystalline  copper.  Petkovic  found  that  some  incubation  time  was  necessary  after 
unloading  before  any  noticeable  recovery  was  obtained  upon  reloading.  This  incuba¬ 
tion  period  extended  to  as  long  as  100  seconds  at  500  degrees  Celsius  (0.6  homologous 
temperature). 

The  question  arises  however  whether  we  should  model  this  behavior.  We  decided  to 
ignore  this  initial  rapid  unloading  to  the  following  reasons: 

1.  Insufficient  recovery  testing  was  performed  to  characterize  the  extent  of  this  effect. 

2.  If  this  initial  high  recovery  rate  is  due  to  microstructure  polarization  or  anelastic 
effects,  we  do  not  expect  these  eifects  to  significantly  alter  the  large  deformation 
material  response. 

3.  If  the  procedure  described  in  Section  4.2  for  evaluation  of  the  static  recovery  func¬ 
tion  is  followed,  the  high  apparent  initial  recovery  rate  produces  an  unreasonably 
strong  dependence  of  recovery  on  the  internal  variable.  This  dependence  corre¬ 
sponds  to  a  power  law  dependence  of  recovery  rate  on  the  internal  variable  s  raised 
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to  a  power  exceeding  100. 


We  therefore  chose  to  consider  data  for  recovery  hold  periods  exceeding  20  seconds  to  fit 
a  material  function  for  static  recovery. 

Following  the  procedure  described  in  Section  4.2,  a  0.002  strain  offset  stress  was 
measured  from  the  reloading  segment  of  each  test.  The  variation  in  this  “yield”  stress 
was  then  used  as  a  measure  of  the  variation  in  state  due  to  different  hold  periods.  This 
data  is  plotted  in  Figure  4.28.  Continuing  with  the  development  presented  in  Section 
4.2,  we  assume  that  the  recovery  function  may  be  represented  by  equation  4.2.6: 

i  =  — r(s)  exp  (-^)  •  4.2.6 

Equation  4.2.11  then  expresses  the  recovery  function  in  terms  of  the  derivative  of  the 
data  in  Figure  4.28: 

ffVcJex  p(-|l)  =  -Jg-.  4.2.11 

d(of/c) 

Figure  4.29  contains  the  data  treated  in  the  manner  suggested  by  4.2.11,  which  then 
represents  the  recovery  function.  Due  to  the  small  number  of  tests  the  derivative  was 
taken  to  be  the  slope  of  the  line  drawn  between  sequential  data  points. 

The  static  recovery  function  which  we  considered  is  a  generalized  form  of  a  function 
proposed  both  by  Nix  and  Gibeling  [1985]  and  Prinz  and  Argon  [1984]  based  on  models 
of  dislocation  climb  causing  static  recovery  in  dislocation  cell  walls: 

r  =  -„Bexp(^)(i)  ,  4.2.12 

where  p  is  the  shear  modulus  and  B  and  p  are  material  parameters.  Both  of  the  above 
pairs  of  investigators  hypothesize  a  value  of  p  between  four  and  six. 

Fitting  equation  4.2.13  to  the  recovery  data  illustrated  in  Figure  4.29  produced  the 
following  values  for  the  material  parameters: 
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B  =  1.08  x  lO60  MPA/sec 

Qr  =  300  kJ/mole, 

p  =  20.6. 

The  procedures  and  numerical  routines  used  to  obtain  these  parameters  are  described 
in  Chapter  5.  Figure  4.30  indicates  the  fit  of  equation  4.2.12  to  the  recovery  data  using 
the  above  parameters.  The  lack  of  correlation  with  the  high  recovery  rate  portion  of  the 
data  reflects  the  previously  described  decision  to  ignore  the  initial  rapid  static  recovery 
data. 

Notice  that  the  exponent  p  is  significantly  greater  than  that  predicted  by  Argon  and 
Nix.  It  is  unclear  whether  this  is  due  to  a  stronger  dependence  of  recovery  on  state  than 
was  anticipated  or  whether  some  characteristic  of  the  experimental  procedure  produces 
this  dependence.  Since  these  experiments  are  rather  arduous  when  performed  with  Fe  - 
2%  Si,  further  investigation  of  both  the  phenomena  and  the  experimental  method  would 
be  easier  if  the  model  material  were  a  lower  melting  temperature  metal  such  as  aluminum. 

Nevertheless,  equation  4.2.12  with  the  above  material  parameters  provides  a  measure 
of  static  recovery,  we  believe  the  first  measure  of  static  recovery  employing  an  internal 
variable  formulation.  To  summarize,  the  static  recovery  function  for  the  internal  variable 
s  is: 


where  B,Qr ,  and  p  are  material  parameters. 

4.3  Evaluation  of  Dynamic  Hardening  and  Recov¬ 
ery  Function:  hU f,  s ,  6) 

If  we  now  consider  an  experiment  where  the  temperature  and  true  strain  rate  are  held 
constant,  then  according  to  section  4.1  for  fully  developed  flows  where  e  «  f  there  is  a 
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unique  proportional  relationship  between  the  equivalent  stress  and  the  internal  variable 
throughout  the  test: 


s  = 


4.3.1 


where  c  is  constant  for  a  given  temperature/strain  rate  pair.  We  can  therefore  express 
the  evolution  equation  for  the  internal  variable  in  terms  of  stress: 


a  =  ch  (a,  — ,  9^jeP  —  cr  ,  o'j  ,  4.3.2 

0" 


The  slope  of  the  stress  versus  plastic  strain  data  (a,  ep),  or  hardening  data,  from  an 
isothermal,  constant  true  strain  rate  test  therefore  reflects  the  internal  variable  evolution 
equation  for  the  above  model.  The  hardening  data  however  incorporates  contributions 
from  both  functions  h  and  r  representing  the  evolution  of  structure.  Additional  testing  or 
assumptions  must  be  made  to  further  separate  these  equations  for  evaluation.  Figure  4.31 
illustrates  an  isothermal,  constant  true  strain  rate  test  path  on  the  schematic  introduced 
in  Figure  4.1. 

A  series  of  isothermal,  constant  true  strain  rate  tests  have  been  performed  on  the  iron 
-  2%  silicon.  These  tests  were  performed  over  a  temperature  range  of  700  to  1200  degrees 
Celsius  (.55  to  .83  Tm)  and  a  strain  rate  range  of  10~3  to  1  sec-5 .  700  Celsius  represented 
the  lowest  temperature  at  which  tests  could  be  performed  due  to  the  load  limit  of  the 
ceramic  load  train.  We  did  not  perform  tests  above  1200  Celsius  due  to  excessive  creep 
of  the  load  train.  The  lower  limit  on  the  strain  rate  of  10"3  sec-1  was  lowest  true  strain 
rate  obtainable  with  the  analog  function  generator.  The  highest  strain  rate  represented 
both  the  limit  of  the  frequency  response  of  the  servohydraulic  and  the  strain  rate  where 
adiabatic  heating  invalidates  isothermal  test  conditions. 
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Figures  4.32  through  4.36  contain  the  stress/strain  data  for  the  different  temperatures 
at  constant  strain  rate,  while  Figures  4.37  through  4.40  illustrate  the  same  data,  except 
for  the  various  strain  rates  at  constant  temperature.  Notice  that  the  stress  eventually 
reaches  an  apparent  steady  state  value  at  strains  approaching  0.5.  Notice  also  that  we 
have  chosen  to  neglect  data  below  an  initial  strain  of  0.01.  Since  the  specimens  are 
grooved  at  the  ends,  we  anticipate  that  at  the  test  start,  some  “settling”  of  the  ends 
will  occur,  rendering  accurate  interpretation  of  the  small  initial  strains  doubtful  if  not 
impossible.  Table  4.2  lists  the  isothermal,  constant  true  strain  rate  test  parameters  and 
values  of  saturation  stress  (<r*)  reached  in  each  test. 

Figures  4.41  through  4.44  represent  the  hardening  data  for  the  constant  true  strain 
rate  tests.  Elastic  strains  were  subtracted  before  the  hardening  (da/dep)  was  determined. 
The  program  used  to  numerically  determine  the  hardening  is  included  in  Appendix  A. 
Elastic  constants  for  the  Fe  -  2%  Si  alloy  are  listed  as  a  function  of  temperature  in  Table 
4.3. 


The  hardening  data  was  used  to  determine  a  functional  form  for  the  evolution  equation 
for  the  internal  variable.  Section  4.0  indicated  that  given  our  chosen  form  for  the  rate 
equation,  we  may  represent  the  evolution  equation  for  the  internal  variable  s  in  terms  of 
hardening  and  stress: 


4.3.3 


Since  we  have  performed  experiments  to  evaluate  the  rate  equation  and  static  recovery 
function,  we  may  adjust  the  hardening  data  for  the  contribution  of  static  recovery  and 
evaluate  the  form  of  h(a,s,0). 

One  important  result  of  the  recovery  tests  is  that  the  static  recovery  function  does 
not  appreciably  affect  the  hardening  data.  The  integrated  contribution  of  the  recovery 
function  to  reducing  the  internal  variable  over  a  typical  isothermal,  constant  true  strain 
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rate  test  is  less  than  1  percent  of  the  instantaneous  value  of  the  internal  variable.  The 
static  recovery  function  therefore  may  be  ignored  rather  than  added  back  to  the  hardening 
data,  and  the  hardening  data,  in  the  absence  of  any  contribution  of  static  recovery,  then 
represents: 


lkKch(d'-c'0)  4-3-4 

The  form  of  h{a,s,9)  which  we  first  investigated  was  the  simple  linear  softening  function 
proposed  by  Anand  [1982]: 


ds 

d& 


4.3.5 


Here,  h0  represents  a  constant  rate  of  athermal  hardening  which  via  mechanistic  argu¬ 
ments  is  generally  assumed  to  be  of  the  order  of  one  hundredth  the  shear  modulus  (/i/100) 
[Kocks,  1966].  The  ratio  {s/s')  represents  deformation  dependent  softening  phenomena. 
The  denominator  s'  is  a  saturation  value  of  s  associated  with  a  given  temperature,  and 
strain  rate  (or  stress);  The  combined  term  ho(s/s')  represents  the  rate  of  dynamic  soft¬ 
ening.  For  isothermal,  constant  true  strain  rate  tests  where  a  =  cs,  equation  4.3.5  may 
be  expressed  as . 


where  a'  =  cs'  is  the  saturation  or  steady  state  stress.  Equation  4.3.5  is  attractive 
on  practical  grounds  since  it  supposes  two  separate  phenomena,  hardening  and  dynamic 
softening,  which  interact  linearly  and  for  which  models  may  be  independently  formulated 
and  then  summed.  The  constant  c  represents  the  contribution  of  constant  structure  defor¬ 
mation  kinetics,  manifested  here  as  a  rate  sensitivity  of  strain  hardening.  Equation  4.3.6 
may  be  integrated  directly  to  yield  the  Voce  [1955]  equation  for  stress/strain  behavior: 


a  =  a'  —  {o'  —  <7q)  exp 


4.3.7 


where  a0  =  cs0  and  s0  is  the  initial  value  of  s. 
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Examination  of  the  hardening  data,  Figures  4.41  to  4.44,  however,  indicates  that  the 
softening  is  not  linear;  there  is  a  substantial  change  in  the  softening  rate  as  the  steady 
state  stress  is  approached.  We  chose  to  accomodate  this  change  in  the  softening  behavior 
via  a  modification  to  equation  4.3.5  through  the  introduction  of  an  exponent  a  to  yield: 

ds 


4.3.8 


or,  for  an  isothermal,  constant  true  strain  rate  test: 

da 


deP  \  a* ) 


4.3.9 


Figure  4.45  indicates  the  behavior  of  the  hardening  version  of  equation  4.3.9  as  tne  value 
of  a  is  varied.  Equation  4.3.9  also  retains  the  feature  t.f  a  saturation  stress,  which  is 
supported  by  the  Fe  -  2%  Si  data,  where  the  stress  appears  to  reach  a  constant  value 
given  isothermal,  constant  strain  rate  conditions.  This  feature  does  not  necessarily  extend 
to  low  («  0.3Tm)  temperatures  where  many  materials  do  not  appear  to  demonstrate  a 
saturation  stress,  even  for  very  large  strains  [Hecker,  1981].  Figure  4.46  indicates  the 
correspondence  between  equation  4.3.9  and  hardening  data  for  a  given  strain  rate  and 
a  range  of  temperatures.  A  value  of  a  =  1.5  was  found  to  best  represent  the  hardening 
curves  for  the  Fe  -  2%  Si. 

For  low  values  of  the  internal  variable  s  (s  <  s*),  the  value  of  h0  may  be  interpreted 
as  the  rate  of  athermal  hardening.  The  value  of  ho  found  for  the  Fe  -  2%  Si  was  3498.0 
MPa,  which  corresponds  to  a  value  of  approximately  /i/20.  This  is  greater  than  the 
hardening  rate  obtained  and  derived  for  State  II  hardening  in  single  crystals  which  as 
stated  previously  is  expected  to  be  in  the  range  of  /i/100  to  /i/200. 


The  integrated  form  of  the  equation  4.3.9  yields: 

5-  =  a*  -  [(<7*  -  cr0)(1-a)  +  (a  -  1)  {(c/io)(d*)-0}  ?] 1_a , 
where  ctq  =  cso  and  So  is  the  initial  value  of  s. 


4.3.10 
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This  formulation  also  requires  an  expression  for  s*,  which  may  be  interpreted  as  the 


saturation  state.  We  propose  a  phenomenological  form  for  sm(0,a,  s,e)  where 


s  =  s 


[ap  ^\ln 

[iexpb?,)J  • 


4.3.11 


This  form  is  motivated  by  requiring  that  the  saturation  state  increases  with  increasing 
strain  rate  and  decreasing  temperature,  and  is  influenced  by  thermally  activated  processes 
represented  by  Qj},  the  activation  energy  associated  with  dynamic  softening. 

Once  again,  since  we  may  represent  s*  by  5^/c,  we  may  plot  s*  versus  strain  rate  in 
Figure  4.47  by  dividing  the  saturation  stress  obtained  during  isothermal,  constant  true 
strain  rate  tests  by  the  constant  c  determined  from  the  rate  equation.  The  figure  indicates 
that  a  power  law  relationship  between  strain  rate  and  saturation  state  is  justified.  Figure 
4.48  illustrates  the  correlation  between  the  steady  state  data  and  a  nonlinear  least  squares 
fit  of  equation  4.3.11  to  the  same  data.  Chapter  5  describes  the  routines  and  procedures 
used  to  fit  the  data  to  the  material  functions.  The  material  constants  obtained  from  the 
nonlinear  least  squares  fit  to  4.3.11  are: 


s  =  36.3  MPa 

n  =  .06272 

Qds  =  439.4  kJ/mole 

The  value  of  the  activation  energy  associated  with  the  steady  state  relation  (439.4 
kJ/mole)  is  almost  twice  that  obtained  for  the  rate  equation  (247.5  kJ/mole).  The 
large  value  of  493.4  kJ/mole  is  an  artifact  of  the  form  of  equation  4.3.11.  If  we  cast  this 
equation  in  the  form: 

~p 

s*  =5  1 

5  5  [a 

then  we  obtain  a  value  of  Qj,  =  27.6  kJ/mole  with  no  change  to  the  values  of  the  other 
two  material  parameters. 


(Qds\ 

exp  (W  ’ 
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Figure  4.49  demonstrates  the  relationship  between  steady  state  stress  and  the  steady 
state  value  of  the  internal  variable  for  different  values  of  strain  rate  and  temperature. 
An  interesting  feature  of  this  plot  is  the  linear  relation  between  state  and  stress  which 
reaches  from  the  power  law  into  the  power  law  breakdown  regime.  One  interpretation  of 
this  result  is  that  power  law  breakdown  is  a  consequence  of  the  kinetics  of  the  deformation 
process  (incorporated  through  the  parameter  c)  and  is  not  due  to  any  change  in  the  basic 
processes  governing  the  evolution  of  state.  As  a  result  power  law  breakdown  should  be 
represented  by  the  rate  equation  in  the  manner  provided  by  equation  4.1.3. 

The  form  of  the  dynamic  term  of  the  evolution  equation  was  finally  taken  to  be: 


h{B,s,o)  =  h0 


4.3.13 


with  sm  given  by  4.3.11.  In  4.3.13  |x|  denotes  the  absolute  value  and  sgn(x)  is  the  sign 
function  defined  as: 

,  v  f  1  x  >  0 

=  {_1  x~0 

The  absolute  value  and  sign  function  are  necessary  to  accomodate  strain  softening  when 
5*  is  less  than  s.  Such  a  situation  is  certainly  possible  in  hot  working  processes  where 
the  strain  rate  decreases  or  the  temperature  increases  at  a  material  point.  Such  a  change 
may  cause  the  current  value  of  s  to  be  greater  than  the  asymptotic  value  represented  by 
s*  calculated  under  the  new  conditions.  Chapter  6  contains  the  results  from  strain  rate 
decrement  tests  which  indicate  the  necessity  of  this  form. 
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4.4  Summary  of  Constitutive  Model 


Once  the  different  material  functions  were  determined,  they  were  combined  into  the  the 
constitutive  equations  summarized  below. 

The  rate  equation: 

f  =  A  exp  (-rs)  [sinh  ((l) 

The  evolution  equation  for  the  internal  variable: 


where 


4.3.11 


The  following  chapter  details  the  procedures  used  to  determine  the  model  parame¬ 
ters  from  a  set  of  appropriate  experimental  data.  The  list  of  model  parameters  are: 
A,  Q,  m ,  £,  h0 ,  a,  s,  Qjs,  n,  B ,  Qr,  and  p.  Also,  p(0)  is  the  elastic  shear  modulus  used  purely 
for  purposes  of  normalizing  the  value  of  s.  Its  temperature  dependence  is  listed  in  Table 


4.3. 
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TABLE  4.1 

Iron  -  2%  Silicon  Jump  Test  Data 


Test 

Number 

Temperature 

(Celsius) 

Initial 
Strain  Rate 
(sec-1) 

Final 

Strain  Rate 
(sec-1) 

Stress 

before  jump 
(MPa) 

Stress 
after  jump 
(MPa) 

1 

700. 

2.00E-4 

2.00E-4 

74.0 

74.0 

2 

2.00E-4 

2.02E-3 

74.0 

90.5 

3 

2.00E-4 

9.42E-3 

74.0 

102.3 

4 

2.00E-4 

8.66E-2 

74.0 

131.0 

5 

2.00E-4 

8.99E-1 

74.0 

182.0 

6 

800. 

2.20E-4 

2.20E-4 

30.0 

30.0 

7 

800. 

2.20E-4 

2.11E-3 

30.0 

41.6 

8 

800. 

2.20E-4 

1.02E-2 

30.0 

49.5 

9 

800. 

2.20E-4 

9.35E-2 

30.0 

65.4 

10 

800. 

2.20E-4 

.953E-0 

30.0 

78.3 

11 

900. 

1.90E-4 

1.90E-4 

12.0 

12.0 

12 

900. 

1.90E-4 

2.07E-3 

12.0 

20.0 

13 

900. 

1.90E-4 

1.97E-2 

12.0 

28.0 

14 

900. 

1.90E-4 

1.01E-1 

12.0 

35.0 

15 

900. 

1.90E-4 

1.11E-0 

12.0 

47.0 

16 

1000. 

2.35E-4 

2.35E-4 

7.0 

7.0 

17 

1000. 

2.35E-4 

2.00E-3 

7.0 

10.6 

1  18 

1000. 

2.35E-4 

1.00E-2 

7.0 

15.5 

19 

1000. 

2.35E-4 

8.79E-2 

7.0 

21.6 

20 

1000. 

2.35E-4 

.899E-0 

7.0 

28.4 
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TABLE  4.2 
Iron  -  2%  Silicon 

Isothermal,  Constant  True  Strain  Rate  Tests 


Test  ID 

Temperature 

(Celsius) 

Strain 

Rate 

(sec-1) 

a* 

(MPa) 

8508132 

700 

1.0 

229.97 

8511291 

750 

1.0 

180.85 

8506191 

800 

1.0 

131.38 

8511281 

800 

1.0 

138.25 

8506201 

800 

.1 

93.73 

8504132 

800 

.1 

91.36 

8054011 

800 

.01 

68.85 

8503311 

798 

.001 

46.16 

8511292 

850 

1.0 

108.44 

8504012 

850 

.01 

49.63 

8506141 

900 

1.0 

84.38 

8502071 

900 

.1 

55.29 

8501231 

900 

.01 

34.77 

8504131 

899 

.01 

36.60 

8501221 

900 

.001 

21.11 

8412151 

950 

.01 

25.18 

8511261 

1000 

1.0 

57.57 

8511301 

1000 

1.0 

53.96 

8503011 

1000 

.316 

44.71 

8501081 

1000 

.1 

32.57 

8502051 

1000 

.0316 

26.78 

8501241 

1000 

.01 

19.77 

8502061 

1000 

.00316 

15.76 

8501121 

1000 

.001 

10.75 

8511271 

1100 

.502 

34.11 

8504041 

1100 

.316 

27.99 

8504061 

1100 

.316 

28.12 

8502081 

1100 

.1 

21.33 

8504021 

1100 

.01 

12.24 

8504121 

1098 

.001 

6.63 

8508131 

1200 

1.0 

25.66 

8508091 

1200 

.1 

14.72 
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TABLE  4.3 

Elastic  Constants  for  Fe  -  2%  Si 


Temperature 

(Celsius) 

(GPa) 

E 

(GPa) 

V 

700 

57 

150 

MU  - 

50 

135 

H: 

42 

115 

.38 

1000 

37 

105 

.41 

1100 

32 

93 

.44 

1200 

30 

87 

.47 

Source:  ASM  Handbook,  Volume  1,  pg  641. 
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Constant  Temperature  Plane 


Figure  4.1  Lines  of  constant  internal  variable  resulting 
from  specification  of  values  of  stress,  strain  rale,  and 
temperature  (Lines  are  straight  only  for  illustration.) 


c  /fi,  10"’ 


Stress  dependence  of  ihe  initial  creep  rate  after 
various  amounts  of  stress  reduction  at  steady-state  creep 
of  Al.  Solid  marks  indicate  initial  states.  (>asfied  lines  are 
the  stress  dependence  of  tt  for  each  temperature  (from 
11541). 


Figure  4.2a  Constant  structure  data  for  aluminum. 
Reference  [Mitra  and  McLean,  1967J 
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-Dependence  of  temperature-compensated  creep  rate 
on  applied  stress  for  specimens  vith  different  structures 
and  internal  stresses  Dashed  and  dotted  curxcs— tempera¬ 
ture  compensated  creep  rate  vs  effective  stress  determined 
by  two  methods. 


Figure  4.2b  Constant  structure  data  for  304  stainless. 
Reference  [Cuddy,  1970) 


Figure  4.2  Constant  substructure  data  for  aluminum. 
Reference  (Gibeling  and  Nix,  1982) 
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Stress  (MPa) 


Figure  4.3  Schematic  of  series  of  strain  rate  jump  tests  where 
jumps  occur  at  the  same  initial  state. 


Constant  Temperature  Plane 


Figure  4.4  Constant  stale,  strain  rate  jump  tests 
illustrated  in  strain  rale,  stress  and  temperature 
space 
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Strain 


Figure  4.5  Strain  rate  jump  tests  with  Fe  -  2%  Si  to  determine  constant 
structure  rate  dependence. 
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Figure  4.6 


Strain  rate  jump  tests  with  Fe  -  27.  Si  to  determine  constant 
structure  rate  dependence. 
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Stress  (MPa) 
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Figure  4.7  Strain  rate  jump  tests  with  Fe  -  2%  Si  to  determine  constant 
structure  strain  rale  dependence. 


Figure  4.8  Strain  rate  jump  tests  with  Fe  -  2%  Si  to  determine  constant 
structure  strain  rale  dependence 
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Log10  Strain  Rate  (seconds-1 
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Logi0  Stress  (MPa) 


Figure  4.12  Stress/strain  rate  dependence  of  Fe  -  2%  Si  at  different 
constant  states  (s,).  The  strain  rate  has  been  normalized  by  an 
Arrhenius  term  with  constant  activation  energy  of  251  kJ/mole, 


Strain  e  Log,0  Strain  Rate  (seconds 


Figure  4.15  Correspondence  between  constant  state  data  and  fit  of  equation 
4.1.3  to  the  same  data. 


Figure  4.16  Schematic  stress  drop  test  illustrating  parameters  used 
by  Mitra  and  McLean  [1966]  to  evaluate  recovery  rale. 


86 


Stress  Rate  (d<r/dt)  Recovery  Time  At 


Stress  Reduction  Act 


Figure  4.17  Illustration  of  back  extrapolation  of  stress  drop  test 
to  determine  static  recovery  rate.  Points  indicated  are  for 
illustration  only;  they  do  not  represent  actual  data. 


0.0 


Strain  Rate  (dy/dt) 

Figure  4.10  Illustration  of  static  recovery  rate  measurement  via  back 
extrapolation  of  strain  rate  jump  data.  Proposed  by  Kocks  [1977], 
Points  are  for  illustration  only  and  do  not  represent  lest  data. 
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Period 


Recovery  Test  Example 

Load-Unload-Hold-Reload  Tests 


Figure  4.19  Schematic  of  load-unload-hold-reload  test  series  used  to 
evaluate  rate  of  static  recovery.  Curves  do  not  represent  actual 
test  data. 


Recovery  Test  Example 

Load-Hold-Unload-Reload  Tests 


Yield  Stress 


Figure  4.20  Hold  period  versus  stress  obtained  immediately  after 

reloading  from  hold  period.  Each  point  represents  one  hold  period. 
Points  are  for  illustration  only  and  do  not  represent  real  data. 


Recovery  Function  (Stress/second) 


Stress  (MPa) 


Stress  (MPa)  Stress  (MPa) 


Figure  4.24  Multiple  hold  test  data  for  Fe  -  2%  Si.  Test  performed  at  1000 
Celsius  and  0.1  sec-1  constant  true  strain  rate.  Duration  of  each  hold 
period  was  20  seconds. 


Strain 


Figure  4.25  Multiple  hold  test  data  for  Fe  -  27,  Si.  Test  performed  at  1100 
Celsius  and  0.1  sec"'  constant  true  strain  rale.  Duration  of  each  hold 
period  was  20  seconds. 
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Fractional  Softening 


Recovery  Function  (MPa/sec) 


Fe  2%Si  Recovery  Tests 

e  =  —  0.02  sec-1 


§  400.0 


Yield  Stress  a,  (MPa) 


Figure  4.28  Hold  period  versus  stress  upon  reloading  for  Fe  -  2%  Si. 
Tests  performed  at  900  degrees  Celsius  and  at  0.02  sec"1  constant 
true  strain  rate. 


Fe  2%Si  Recovery  Tests 

e  =  -  0.02  sec'1 


State  Variable  (s,  =  o>/c)  (MPa) 


Figure  4.29  Recovery  function  dependence  on  slate  for  Fe  -  Z% 
Data  obtained  from  Figure  4.28  according  to  equation  4.2.11 


Recovery  Function 


Strain 


Stress  (MPa) 


Strain 


Stress  (MPa)  Stress  (MPa) 


Strain 


Figure  4.38  Isothermal,  constant  true  strain  rate  tests  on  Fe  -  2%  Si  with 
strain  rate  of  -0.1  sec"1  and  at  different  temperatures. 


Strain 

Figure  4.39  Isothermal,  constant  true  strain  rate  tests  on  Fe  -  2%  Si  with 
strain  rate  of  -0.01  sec"1  and  at  different  temperatures. 
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Strain 


Figure  4.40  Isothermal,  constant  true  strain  rate  tests  on  Fe  -  2%  Si  with 
strain  rate  of  -0.001  sec  _1  and  at  different  temperatures. 


Stress  (MPa) 


dcr/dep  (MPa)  dcr/dep  (MPa) 


Figure  4.43  Plastic  hardening  data  (da/dep)  derived  from  isothermal, 
constant  true  strain  rale  tests  on  Fe  -  2%  Si.  Strain  rate  of  -0.01  sec'1. 


Figure  4.44  Plastic  hardening  data  (do/de!>)  derived  from  isothermal, 
constant  true  strain  rate  tests  on  Fe  -  2%  Si.  Strain  rale  of  -0.001  sec 


(da/dep)/ch 


1.0 


Figure  4.45  Variation  in  equation  4.3.9  as  the  value  of  exponent  a 
is  varied. 


Stress  (MPa) 


Figure  4.46  Correspondence  between  hardening  data  and  equation  4.3.9 
for  a  value  of  a=1.5. 


LogI0  (s  =  or  / c) 


-4.0  -3.0  -2.0  -1.0  0.0  1.0 

Log10  Strain  Rate:  e„ 

Figure  4.47  Saturation  value  of  internal  variable  derived  via 
s'  =  a  / c  versus  strain  rate  for  Fe  -  2%  Si.  Data  obtained 
from  isothermal,  constant  true  strain  rate  tests. 


-4.0  -3.0  -2.0  -1.0  0.0  1.0 

Logi0  Strain  Rate:  e„ 


Figure  4.40  Fit  of  equation  4.3.11  for  saturation  of  internal  variable 
data  for  Fe  -  2%  Si.  Data  obtained  from  isothermal,  constant  true 
strain  rate  tests. 
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Chapter  5 

Determination  of  Material 
Parameters 

Fitting  of  the  parameters  associated  with  our  constitutive  model  follows  the  progression  of 
testing  described  in  the  previous  two  chapters.  Parameters  associated  with  each  material 
response  function  are  determined  from  the  tests  used  to  characterize  that  particular 
aspect  of  material  behavior.  The  rate  equation  parameters  are  therefore  determined  from 
the  jump  test  data,  static  recovery  function  parameters  are  determined  from  the  load- 
unload- hold- reload  data,  and  the  dynamic  hardening  and  recovery  function  is  determined 
from  the  isothermal,  constant  true  strain  rate  tests. 

The  fitting  procedure  exploits  extensively  the  proportionality  between  stress  and  the 
internal  variable  resulting  from  the  consequences  of  an  isothermal,  constant  true  strain 
rate  test.  At  a  fixed  temperature  and  a  fixed  compressive  strain  rate,  equation  4.1.3 
implies: 

s  =  |cr|/c,  5.1 

where  c  is  a  constant  defined  by 


This  relationship  permits  the  substitution  of  (|<rj/c)  for  the  internal  variable  s  while 
determining  model  parameters.  Since  the  internal  variable  s  represents  some  generalized 


resistance  to  plastic  deformation  it  is  difficult,  if  not  practically  impossible,  to  measure  a 
value  of  s  directly.  Ideally  we  would  like  to  have  some  macroscopic  or  microscopic  quantity 
which  we  may  correlate  with  s.  Some  possibilities  for  such  measures  are  discussed  in 
Chapter  6.  In  the  absence  of  a  directly  measureable  quantity,  equation  5.1  permits  the 
determination  of  model  parameters  without  directly  knowing  the  value  of  the  internal 
variable. 

The  following  procedure  provides  a  deterministic  method  for  evaluating  the  material 
parameters  associated  with  the  constitutive  model  presented  in  Chapter  4.  The  stress 
after  a  change  in  strain  rate,  whether  it  be  associated  with  a  strain  rate  jump  or  a 
recovery  test  reloading,  is  determined  in  our  data  as  the  stress  corresponding  to  a  0.2% 
strain  offset  from  the  strain  at  the  change  in  strain  rate. 


1.  Determination  of  rate  equation  parameters  A,Q  and  m : 

The  strain  rate  jump  tests  provide  constant  state  data  for  different  states  at  different 
temperatures.  If  we  normalize  the  stresses  after  each  jump  by  the  stress  before  each 
jump,  then  the  rate  equation  eliminates  the  contribution  of  the  internal  variable 
since  the  test  before  and  after  the  jump  is  both  isothermal  and  at  a  constant  true 
strain  rate: 


Of  Cj  s 
0{  Ci  s  ’ 


5.3 


or, 


where 


sinh-1 

{^exp 

(w) 

FI 

sinh-1 

,{iex  pi 

n 

5.4 


o j  =  stress  immediately  after  jump, 

Oi  =  stress  immediately  before  jump, 
e/  =  strain  rate  after  jump, 
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q  =  strain  rate  before  jump,  and 

0  =  test  temperature. 


Equation  5.4  is  fit  to  the  jump  test  data  (<7/,cr,-,e/,et-,0)  for  each  temperature 
through  a  nonlinear  least  squares  fit  to  provide  the  parameters  m  and  the  combined 


as  a  function  of  temperature.  The  pre-exponential  A  and  the  activation  energy  Q 
are  then  determined  via  a  second  nonlinear  least  squares  fit  of  the  temperature  to 
the  above  term.  The  constant  m  is  taken  as  the  average  value  across  all  temper¬ 
atures.  Nonlinear  fitting  throughout  this  procedure  was  performed  with  an  IMSL 
numerical  subroutine  ZXSSQ  which  provides  for  minimization  of  a  nonlinear  ob¬ 
jective  function  on  the  basis  of  a  least  squares  residue.  The  programs  used  in  the 
fitting,  ISOCON  and  its  associated  subroutine  FITISO,  are  listed  in  Appendix  B. 

The  three  constants  ( A,Q,m )  could  be  determined  simultaneously  via  a  three  pa¬ 
rameter  least  squares  fit,  but  it  was  found  that  the  two  step  procedure  produced  a 
better  correlation  with  the  data. 

Figure  5.1  indicates  the  correlation  between  the  jump  test  data  and  the  rate  equa¬ 
tion  prediction  using  parameters  determined  via  the  above  process. 

It  should  be  stated  that  the  procedure  described  in  this  step  does  not  produce  an 
unambiguous  value  for  the  activation  energy,  Q.  Strain  rate  jump  tests  directly 
measure  the  functional  relation  between  stress  and  strain  rate.  Temperature  jump 
tests  are  necessary  if  one  wishes  to  measure  the  activation  energy  directly.  We  are 
able  to  determine  an  activation  energy  here  since  the  rate  equation  cast  in  the  form 
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of  equation  5.4  retains  a  variation  with  temperature.  A  different  rate  equation 
may  very  well  not  permit  the  determination  of  an  activation  energy  from  strain 
rate  jump  test  data.  Practically,  very  little  error  is  introduced  to  the  model  if  one 
chooses  to  use  the  same  value  of  activation  energy  for  both  Q  and  Qds. 


2.  Determination  of  state/stress  scaling  factor  £: 

As  was  indicated  in  Section  4.1  the  parameter  (  is  a  dimensionless  scaling  factor 
included  in  the  rate  equation  to  enforce  our  requirement  that  the  value  of  the 
internal  variable  $  always  be  greater  than  the  value  of  the  equivalent  tensile  stress 
a.  Motivated  by  an  interpretation  of  s  to  represent  some  athermal  measure  of 
deformation  resistance,  this  corresponds  to  requiring  the  proportionality  constant 
c  to  be  less  than  one,  where  a  —  cs  for  isothermal,  constant  true  strain  rate 
conditions.  The  parameter  does  not  improve  the  predictive  capability  of  the  model; 
it  may  be  set  equal  to  one  with  no  adverse  effect.  A  reasonable  value  of  f  may 
therefore  be  calculated  by  requiring: 


£  >  sinh  1 


where  the  strain  rate  and  temperature  are  the  highest  and  lowest  values  of  their 
applicable  ranges,  respectively. 


3.  Determination  of  static  recovery  term  parameters  B,  Qr  and  p : 

The  static  recovery  function  parameters  were  determined  via  a  nonlinear  least 
squares  fit  of  equation  4.2.13  to  the  recovery  data  analyzed  in  the  manner  described 
in  Section  4.2,  where  4.2.13  is: 


-d/St 

d{aj/c) 


_  v  (Q r\  /*p 
6XP  UJ  fiB  \c  )  ' 


4.2.13 


Here,  c  is  the  stress/internal  variable  constant  defined  by  5.2.  We  assumed  here 
that  Qr  equals  the  activation  energy  for  rdf-diffusion  of  alpha  iron,  300  kJ/mole. 
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[Frost  &  Ashby,  1982]  We  provide  two  justifications  for  this  assumption.  One, 
static  recovery  has  been  observed  to  be  associated  with  non-conservative  climb 
of  dislocations,  [Takeuchi  and  Argon,  1976]  the  relevant  activation  energy  then  is 
that  of  self-diffusion.  Two,  we  checked  this  assumption  via  the  following  test.  We 
performed  two  load-unload-hold-reload  tests  at  the  same  value  of  c  but  at  different 
temperatures  (and  consequently  different  strain  rates).  The  strain  rates  needed 
to  produce  the  same  value  of  c  may  be  determined  from  the  rate  equation  with 
the  parameters  determined  in  Step  1.  If  we  design  the  tests  to  obtain  the  same 
values  of  stress  just  before  unloading  (a;)  and  immediately  after  reloading  (07) 
then,  according  to  our  model,  the  value  of  the  internal  variable  is  the  same  both  at 
the  beginning  and  at  the  end  of  the  hold  period  for  both  tests.  If  these  conditions 
are  met,  then  for  the  two  tests  at  the  two  temperatures  according  to  4.3.17: 

exp  (-;§;)  A‘1=exp  HI) A^  4A14 

Tests  were  performed  according  to  the  above  procedure,  varying  the  hold  times 
until  the  same  value  of  07  was  obtained.  The  two  tests  used  for  the  measurement 
of  Qr  possessed  the  following  experimental  parameters: 


Test 

Temperature 

(Celsius) 

Strain  rate 
(1/sec) 

Hold  strain 

o'; 

(MPa) 

(MPa) 

131 

1 

800. 

■ 

46.7 

352. 

2 

.1 

46.5 

40.8 

20. 

Equation  4.2.14  using  the  values  from  the  above  table  produces  a  value  of  Qr  = 
300.2  kj/mole,  in  excellent  agreement  with  our  assumed  valur  for  self-diffusion. 

The  recovery  data  used  to  fit  4.2.11  is  plotted  in  Figure  4.29,  and  listed  below  in 
Table  5.1.  Fitting  was  accomplished  via  a  nonlinear  least  squares  fit  similar  to  that 
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used  in  Step  1.  The  programs  used  for  the  fitting,  STATCON  and  FITREC,  are 
similarly  listed  in  Appendix  B. 


Table  5.1 
Recovery  Test  Data 
900  Celsius  -  0.02  1/sec 


At 

(seconds) 

(MPa) 

(—dAt)/d(af/c) 

(seconds/MPa) 

20 

33.6 

2.67 

70 

33.1 

31.25 

150 

33.5 

42.19 

300 

32.9 

78.74 

500 

32.2 

86.96 

The  correlation  between  the  recovery  data  and  the  recovery  function  prediction  is 
plotted  in  Figure  5.2.  Notice  that  the  exponent  in  this  case  exceeds  20,  significantly 
greater  than  the  exponent  of  4  proposed  by  Argon  [Prinz  and  Argon,  1984]  and 
Nix  [Nix  and  Gibeling,  1985].  The  point  corresponding  to  the  shortest  hold  period 
was  ignored  in  the  curve  fitting.  We  associated  the  rapid  initial  stress  decrease 
accompanying  short  hold  times  with  very  small  scale  rearrangement  or  anelastic 
effects.  One  consequence  of  including  this  effect  would  be  to  increase  the  value 
of  the  exponent  p  to  magnitudes  of  50  to  100.  We  therefore  chose  to  ignore  this 
effect  with  the  consequence  of  representing  longer  time  recovery  effects  with  greater 
accuracy  than  short  time  effects. 

It  should  be  noted  that  one  may  evaluate  the  static  recovery  function  using  offset 
stresses  measured  at  a  larger  strain  offset  than  the  0.002  strain  offset  used  above. 
This  was  done  with  the  Fe  -  2%  Si  data  using  a  0.02  strain  offset.  The  value  of 
the  exponent  p  resulting  from  the  analysis  exceeded  thirty.  Selection  of  a  different 
strain  offset  therefore  did  not  reduce  the  relatively  high  magnitude  of  the  state 
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variable  exponent. 


4.  Determination  of  the  Saturation  Stress  a*: 

The  value  of  the  saturation  stress  cr*  is  determined  for  each  test  from  the  steady 
state  value  reached  in  the  (cr,  ep)  data1. 

5.  Determination  of  a, cho,  and  <r0: 

The  combined  constant  cho  and  the  parameter  a  for  each  test  are  determined  by 
a  least  squares  fit  of  the  hardening  (<r,da/dep)  data  according  to  equation  4.3.9, 
restated  here: 


After  the  values  of  a  and  ch0  are  obtained  the  value  of  a0  is  obtained  by  a  least 
square  fit  of  the  constant  strain  rate  (cr,  ep)  data  to  equation  4.3.10  for  fixed  values 
of  a  and  ch0.  Alternatively,  cho,  a,  and  cr0  may  be  determined  by  a  similar  least 
squares  fit  of  the  (cr,  ep)  data  to  the  integrated  stress-strain  relation  of  equation 
4.3.10.  We  believe  that  the  first  approach  is  more  physically  motivated,  but  the 
second  should  work  well,  especially  when  it  is  difficult  because  of  lack  of  sufficient 
data  points  or  noisy  data  to  differentiate  numerically  the  (cr,  ep)  data.  Normally 
there  is  variation  in  the  values  of  a  determined  from  different  tests.  The  value 
of  a  is  therefore  determined  as  the  average  of  the  values  determined  from  each 
individual  test.  After  the  value  of  a  is  fixed,  the  least  squares  fit  of  equation  4.3.9 
is  performed  to  determine  the  final  value  of  cho.  For  the  2%  silicon  iron  the  value  of 
a  =  1.5  describes  our  experimental  data  well.  Accordingly,  for  this  material  we  fix 
a  =  1.5.  The  range  of  data  used  to  determine  ch0  and  <r0  corresponds  to  ep  >  1% 
because  there  are  various  uncertainties  associated  with  measuring  small  strains 

'If  a  steady  state  is  not  reached  in  the  experiments,  then  <r*  is  estimated  by  extrapolating  the  a 
versus  dcr/dtp  data  to  da/dcp  =  0. 
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in  compression  by  our  test  techniques.  Figure  4.46  illustrates  the  correspondence 
between  the  hardening  data  and  equation  4.3.9  using  a  =  1.5.  The  program  used 
in  this  step,  QNON,  is  listed  in  Appendix  B. 


6.  Determination  of  h0: 

The  value  of  ho  is  then  obtained  as  the  simple  average: 

,  a.,  «<*>).•/«) 

ho  — - , 


(5.1) 


where  [c/r0],  is  the  value  of  cho  determined  in  step  5  for  test  i,  c,  is  the  correspond¬ 
ing  value  of  c  defined  in  5.2,  and  q  is  the  number  of  tests.  The  value  of  ho  was 
determined  to  be  3498.0  MPa. 

7.  Determination  of  s,n,and  Qd, : 

The  values  of  s,  n  and  the  activation  energy  Qja  are  then  determined  by  a  non-linear 
least  squares  fit  to  equation  4.3.11  cast  in  terms  of  stress: 


cs 


&  ( Qds 

A  eXp  \  R0 


exP  l  -£ir)  • 


The  programs  used  for  this  fit  are  named  SATCONS  and  its  associated  subroutine 
FITSAT.  They  are  included  in  Appendix  B. 

8.  Determination  of  So: 

The  value  of  so  for  each  test,  the  initial  value  of  the  deformation  resistance  s,  is  then 
determined  by  dividing  the  value  of  oo  determined  in  step  5  by  the  corresponding 
value  of  c  defined  in  equation  5.2.  Table  5.2  provides  the  average  so  values  for  each 
test  temperature,  determined  via  the  procedure  described  above. 

The  material  parameters  for  the  2%  silicon  iron  obtained  by  the  method  outlined  above 
are  listed  in  Table  5.3.  The  value  of  the  activation  energy  Q  =  247.5  kJ/molc  is  within 
the  range  237  to  251  kJ/mole  for  the  activation  energy  for  lattice  self-diffusion  in  ferrite 
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(see,  for  example  Ashby  and  Frost  [1982]).  We  note  that  for  a  similar  silicon  steel,  Uvira 
and  Jonas  [1968]  have  previously  obtained  a  value  of  Q  =  333.6  kJ/mole.  As  discussed 
in  Chapter  3,  however,  there  is  no  reason  to  expect  the  activation  energy  associated  with 
the  rate  equation  to  be  equal  to  either  that  of  self-diffusion  or  that  obtained  from  a  fit 
to  steady  state  data. 

Although  there  is  no  physically  justified  reason  for  doing  so,  one  may  set  Qj,  equal  to 
the  value  of  Q  determined  for  the  rate  equation  without  introducing  a  significant  amount 
of  error  in  the  predictive  capability  of  the  model.  The  different  values  of  activation  en¬ 
ergy  in  the  two  equations  reflect  a  difference  in  the  physical  mechanisms  governing  that 
particular  aspect  of  the  material  response.  The  model,  however,  does  not  appear  to  be 
extremely  sensitive  to  averaging  these  two  values  of  activation  energy.  We  shall  adopt 
this  simplification.  With  this  simplification  the  material  parameters  are  those  given  in 
Table  1  of  Chapter  1,  instead  of  those  in  Table  5.3.  Note  that  it  is  not  only  the  activation 
energies  that  have  changed.  A  change  in  the  activation  energies  causes  a  commensurate 
change  in  the  other  material  parameters  as  well. 

5.1  Minimum  number  of  tests 

One  reasonable  question  associated  with  the  determination  of  material  parameters  is: 
What  is  the  minimum  number  of  tests  that  need  be  performed  to  provide  sufficient  data 
to  follow  the  above  procedure?  Certainly  one  need  not  perform  as  many  tests  as  are 
presented  in  this  report. 

If  we  ignore  the  static  recovery  function,  it  is  possible  to  determine  all  of  the  constants 
from  a  series  of  strain  rate  jump  tests  like  those  presented  in  Section  4.1.  The  jump  seg¬ 
ments  provide  constant  structure  data  for  determining  the  parameters  associated  with 
the  strain  rate  equation.  We  may  determine  the  parameters  associated  with  the  satura- 
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tion  state  equation  if  we  perform  the  tests  such  that  the  stress  reaches  a  constant  value 
before  the  jump  and  then  allow  the  stress  to  reach  a  new  steady  state  after  the  jump. 
Finally,  the  hardening  segments  of  the  curves  before  and  after  each  jump  provides  the 
parameters  for  the  evolution  equation  for  the  internal  variable. 

It  is  possible  (although  certainly  not  recommended)  to  determine  all  of  the  material 
parameters  from  six  (6)  strain  rate  jump  tests,  divided  into  two  sets  of  three  jump  tests 
with  each  set  performed  at  a  different  temperature.  This  set  of  tests  represents  the 
minimum  number  required  to  uniquely  determine  each  material  parameter,  once  again 
assuming  that  we  have  neglected  static  recovery.  Obviously,  more  tests  are  required  to 
provide  confidence  in  the  quality  of  the  correlation  between  material  and  model,  with 
appropriate  tests  spanning  the  range  of  strain  rates,  stresses,  and  temperatures  to  be 
represented  by  the  model. 
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TABLE  5.2 

Initial  Values  of  Internal  Variable  s 


Temperature 

Celsius 

so 

MPa 

700 

134.2 

750 

121.1 

800 

102.7 

850 

89.6 

900 

74.6 

950 

79.1 

1000 

66.1 

1100 

58.4 

1200 

60.1 

TABLE  5.3 

Material  Parameters  for  2%  Silicon  Iron 


Material  Parameter 

Value 

A 

1.26  x  108  sec-1 

Q 

247.5  kJ/mole 

e 

5.00 

m 

0.22793 

8 

36.3  MPa 

Qdt 

439.4  kJ/mole 

n 

0.06272 

ho 

3498.0  MPa 

a 

1.5 

B 

1.08  x  1060  (MPa/sec) 

Qr 

300.0  kJ/moie 

P 

20.6 
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Recovery  Function  (MPa/sec)  Logl0  Slrain  Rale  (seconds-t; 


0.5  1.0  1.5  2.0  2.5 


Log,,,  Stress  (MPa) 

Figure  5.1  Correspondence  between  constant  stale  data  and  fit  of  equation 
4.1.3  to  the  same  data. 


Figure  5.2  Fit  of  power  law  static  recovery  function  to  Fe  ~  2%  Si 
recovery  test  data. 


Chapter  6 

Evaluation  and  Validation  of 
Constitutive  Model 


No  consistent  criteria  have  been  established  for  the  evaluation  of  different  constitutive 
models;  most  comparisons  to  date  consist  of  checking  whether  models  represent  basic 
qualitative  aspects  of  material  behavior.  Nevertheless,  the  evaluation  and  validation  of 
a  model  involves  several  considerations  which  include: 

1.  At  the  very  least,  the  model  should  represent  the  physical  data  upon  which  the 
model  is  based  and  from  which  the  material  parameters  of  the  model  are  deter¬ 
mined. 

2.  Comparison  between  the  predictions  from  the  model  and  the  physical  system  should 
encompass  the  variable  space  to  be  modelled.  Given  a  set  of  external  and  internal 
variables,  the  model  should  represent  the  physical  system  thoughout  all  permuta¬ 
tions  of  those  variables.  If  the  model  parameters  have  been  evaluated  using  data 
from  a  subset  of  the  variable  space,  validation  should  include  not  only  that  subset 
but  as  much  of  the  variable  space  as  possible.  Validation  should  therefore  examine 
the  predictive  capability  of  the  model  at  reasonable  points,  including  the  limits,  of 
the  variable  space. 
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3.  Item  one  above  supposes  that  all  of  the  relevant  variables  are  measurable.  Unfor¬ 
tunately  in  our  case  there  is  no  direct  correlation  of  the  internal  variable  s  to  a 
physical  quantity.  One  consideration  therefore  in  our  investigation  is  whether  we 
may  find  some  some  means  of  correlating  s  to  some  physical  measurement  which 
we  may  then  take  to  represent  s. 

The  tests  which  follow  provide  initial  validation  of  the  model  with  consideration  of  the 
above  issues. 

Five  sets  of  tests  follow.  First,  the  model  is  correlated  with  the  original  material  tests, 
and  the  correlation  is  found  to  be  excellent.  Second,  a  load  control  test  is  compared  with 
a  prediction  from  the  constitutive  model.  Data  obtained  from  the  load  controlled  exper¬ 
iments  was  not  used  in  the  evalutation  of  the  material  parameters  for  the  model.  Third, 
strain  rates  jump  tests  where  the  strain  rate  is  increased  twice  during  one  test  is  com¬ 
pared  with  the  predicted  model  response.  Fourth,  a  strain  rate  decrement  test  is  used 
to  confirm  the  presence  of  strain  softening.  Fifth,  a  specimen  is  proposed  and  tested 
to  investigate  both  the  ability  to  predict  a  variation  in  state  subject  to  an  inhomoge- 
nous  deformation  and  to  investigate  different  macroscopic  measurements  which  may  be 
correlated  with  the  internal  variable. 

6.1  Numerical  Integration  of  Constitutive  Equations 

The  simulations  presented  in  this  chapter  were  produced  via  one  of  two  integration  pro¬ 
cedures,  both  of  which  are  represented  by  programs  listed  in  Appendix  C.  The  first 
program,  ISOINT,  uses  an  IMSL  numerical  subroutine  to  evaluate  all  tests  where  the 
true  strain  rate  is  constant  or  is  instantaneously  changed.  Given  a  constant  true  strain 
rate  and  isothermal  conditions,  integration  of  the  constitutive  model  reduces  to  integra¬ 
tion  of  the  evolution  equation  for  the  internal  variable.  Although  the  coupled  system  of 
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equations  is  relatively  “stiff”,  the  evolution  equation  for  s  is  not.  Therefore  the  IMSL 
subroutine  DVERK  which  employs  a  fifth  order  Runge-Kutta  algorithm  for  non-stiff  sys¬ 
tems  of  differential  equations  was  employed.  A  typical  simulation  required  less  than  5 
seconds  of  CPU  time  on  a  Data  General  MV4000. 

Simulations  of  experiments  involving  other  boundary  conditions  were  performed  us¬ 
ing  the  finite  element  program  ABAQUS.  The  constitutive  equations  obtained  in  this 
investigation  have  been  incorporated  by  Anand  et  al.  [1985]  as  a  material  model  via  a 
user-material  interface  in  ABAQUS  developed  by  Hibbitt,  Karlsson,  and  Sorensen,  Inc. 
of  Providence,  Rhode  Island.  The  coupled  differential  constitutive  equations  are  nu¬ 
merically  very  stiff  and  require  special  computational  schemes  for  efficient  integration  in 
time.  Anand  et  al.  use  a  semi-implicit  time-integration  procedure  with  automatic  time 
stepping.  Using  this  time-integration  procedure  in  ABAQUS  we  may  evaluate  material 
response  in  arbitrary  three  dimensional  loading  conditions.  A  simulation  of  an  isother¬ 
mal,  constant  true  strain  rate  test  using  ABAQUS  requires  approximately  10  minutes  of 
CPU  time  as  opposed  to  the  5  seconds  required  by  ISOINT. 

6.2  Simulation  of  Original  Experiments 

Figures  6.1  through  6.14  present  isothermal,  constant  true  strain  rate  test,  recovery  test, 
and  jump  test  simulations  with  the  original  experimental  data.  All  simulations  were 
produced  using  ISOINT.  The  initial  values  of  the  internal  variable,  s0,  used  in  these 
simulations  varied  with  temperature  and  are  included  in  Table  5.2.  The  correlation 
between  the  model  and  the  experimental  data  is  excellent.  The  only  large  discrepancy 
between  the  model  and  the  experiments  occurs  with  the  700  degree  Celsius,  strain  rate 
jump  tests  (Figure  6.10).  This  is  not  unanticipated  since  the  jump  test  data  was  not  well 
represented  by  the  rate  equation  at  that  temperature. 


119 


6.3  Load  Boundary  Condition  Test 

Since  the  independent  variables  in  the  major  portion  of  the  experiments  were  strain  rate 
and  temperature,  it  is  important  to  investigate  the  predictive  capability  of  the  constitu¬ 
tive  model  when  they  are  dependent  variables.  Experiments  which  permit  such  compar¬ 
ison  include  load  or  stress  controlled  tests  and  temperature  variation  tests. 

A  simple  isothermal,  load  controlled  compression  experiment  was  selected  where  the 
load  was  increased  linearly  at  a  rate  of  2000  Newtons  per  second  up  to  20,000  Newtons, 
held  for  10  seconds,  and  then  decreased  linearly  to  zero  load,  illustrated  in  Figure  6.15. 
Displacement  was  recorded  versus  time,  the  true  strain  and  stress  calculated,  and  the 
strain  rate  determined  from  the  strain  versus  time  data.  Total  accumulated  compressive 
strain  in  this  test  exceeded  -0.9.  Figure  6.16  illustrates  the  agreement  between  the  ex¬ 
perimental  strain  rate  and  predicted  strain  rate  versus  time.  The  ABAQUS  input  data 
for  the  simulation  is  included  in  Appendix  D. 

Although  the  initial  loading  segment  of  the  test  produces  the  worst  agreement  between 
the  model  and  the  data,  it  is  also  the  most  severe  test  of  the  model,  since  the  model  is 
being  asked  to  follow  a  rate  of  change  in  one  of  the  state  variables.  The  constant  load 
segment,  where  the  applied  stress  varies  more  slowly,  is  duplicated  very  well  by  the  model. 

6.4  Strain  Rate  Decrement  and  Double  Jump  Tests 

Figures  6.17  to  6.19  show  variations  in  the  strain  rate  jump  tests,  where  the  strain 
rate  is  in  the  first  case  increased  twice  and  in  the  second  case  decreased  once.  The 
model  represents  the  double  jump  test  data  very  well,  duplicating  both  hardening  and 
instantaneous  strain  rate  dependence. 

Although  the  constitutive  model  predicts  the  final  steady  state  stress  in  the  strain 
rate  decrement  test,  the  model  predicts  a  transition  to  that  stress  after  the  strain  rate 
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decrement  which  is  much  different  from  the  actual  data.  The  experimental  data  changes 
abruptly  to  the  new  value  of  steady  state  stress.  This  behavior  was  unexpected.  It  is 
not  clear  whether  this  represents  true  material  behavior  or  is  some  experimental  arti¬ 
fact.  Similar  data  obtained  with  identical  tests  on  1100  aluminum  [Kim  et  al.,  1986] 
demonstrated  a  material  response  like  that  predicted  by  the  constitutive  model.  It  is 
possible  that  the  abrupt  transition  is  a  characteristic  of  the  analog  function  generator  or 
the  servobydraulic,  although  preliminary  examinations  of  the  hardware  failed  to  detect 
any  abnormalities.  It  would  be  useful  to  perform  the  same  tests  on  some  other  metals 
using  different  hardware  to  determine  whether  this  result  is  reproducible. 

The  steady  state  stress  reached  after  the  strain  rate  decrement  corresponds  to  the 
steady  state  value  for  a  monotonic  test  at  the  final  strain  rate  and  temperature.  The 
sgn(-)  function  adopted  in  the  evolution  equation  for  the  internal  variable  simulates 
strain  softening,  permitting  the  model  to  predict  the  new  saturation  state. 

6.5  State  Gradient  Specimen  and  Testing 

One  assumption  of  an  internal  variable  constitutive  model  is  that  the  value  of  the  internal 
variable  provides  some  measure  of  the  mechanical  state  of  the  material.  The  aim  of  this 
particular  test  was  to  determine  the  correlation  between  the  value  of  the  internal  variable 
predicted  by  the  constitutive  model  and  the  variation  in  deformation  resistance  within 
an  inhomogeneously  deformed  test  specimen.  A  simple  conical  compression  specimen 
geometry  was  selected  to  produce  a  gradient  of  deformation  along  the  vertical  axis.  Figure 
6.20  shows  a  typical  gradient  specimen  before  and  after  testing. 

A  typical  test  sequence  consisted  of  compressing  the  specimen  at  elevated  tempera¬ 
ture,  immediately  quenching,  and  then  machining  a  series  of  small  cylindrical  compression 
specimens  from  the  center  axis  of  the  deformed  specimen.  The  yield  stress  obtained  from 
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a  compression  test  on  each  compression  specimen  provides  a  measure  of  the  deformation 
resistance  in  the  region  of  the  gradient  specimen  from  which  the  compression  specimen 
was  machined.  Microhardness  tests  were  also  performed  along  the  center  section  of  the 
compressed  specimen  to  evaluate  whether  microhardness  could  be  correlated  with  the 
internal  variable. 

Since  we  have  evaluated  the  internal  variable  s  by  measuring  a  0.2%  yield  stress,  and 
since  the  yield  stress  is  temperature  dependent,  ideally  we  would  perform  these  compres¬ 
sion  experiments  at  the  original  gradient  test  temperature.  This  is  operationally  difficult 
to  do  since  the  specimen  requires  a  finite  amount  of  time  to  heat  to  test  temperature, 
and  recovery  processes  may  erase  the  variation  in  deformation-induced  microstructure 
between  specimens.  We  may  perform  these  yield  tests  at  room  temperature,  but  then 
we  are  faced  with  the  dilemma  of  relating  the  room  temperature  yield  stress,  where  the 
high  temperature  value  of  c  may  not  be  applicable,  to  the  deformation  resistance. 

A  set  of  initial  tests  were  performed  using  Fe  -  2%  Si  gradient  specimens.  The  spec¬ 
imens  were  heated  using  an  induction  coil  and  then  cooled  with  a  helium  gas  quench 
immediately  after  compression.  Room  temperature  yield  tests  using  compression  speci¬ 
mens  taken  from  the  axis  of  the  quenched  specimens  produced  no  discernable  variation 
in  yield  stress  as  a  function  of  axial  position  in  the  gradient  specimen.  It  was  uncertain 
whether  the  lack  of  variation  was  due  to  inadequate  quench  iate&,  insufficient  grains  per 
compression  specimen  cross-section,  or  the  true  absence  of  a  variation  in  state. 

We  then  decided  to  perform  an  identical  set  of  tests  using  1100  aluminum  in  place  of 
the  Fe  -  2%  Si.  The  use  of  aluminum  permitted  large  gradient  specimens,  faster  quench 
rates,  and  easier  machining  and  testing  of  the  yield  stress  compression  specimens.  A 
set  of  parallel  experiments  following  the  methods  applied  to  the  Fe  -  2%  Si  have  been 
performed  by  Kwon  Hee  Kim  [Kim  et  al.  1986]  to  determine  the  material  parameters 
for  1100  aluminum  for  the  model  developed  in  this  investigation.  The  aluminum  could 
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be  tested  in  an  air  atmosphere  using  a  radiant  lamp,  clamshell  furnace  which  permitted 
heating  rates  more  than  two  orders  of  magnitude  greater  than  was  possible  using  the  vac¬ 
uum  furnace.  The  radiant  furnace  also  permitted  a  water  quench  which  also  increased 
quenching  rates.  Aluminum  thus  permitted  both  the  possibility  of  a  more  reliable  reten¬ 
tion  of  high  temperature  microstructure  and  the  potential  of  performing  the  yield  tests  at 
elevated  temperatures.  Figure  6.21  shows  the  aluminum  gradient  specimen  dimensions. 

The  gradient  specimen  compression  test  was  simulated  using  ABAQUS  and  the  user 
material  subroutine  briefly  described  in  Section  6.1.  Appendix  D  contains  the  ABAQUS 
input  file  describing  the  finite  element  model.  Figure  6.22  illustrates  the  displacements 
resulting  from  a  simulation  of  a  test  on  the  1100  aluminum  at  300  degrees  Celsius  and  at 
a  compression  rate  of  1  millimeter  per  second.  Figure  6.23  provides  a  contours  of  of  the 
internal  variable  for  this  same  test.  Values  provides  in  the  contour  table  are  in  terms  of 
equivalent  shear  stress  and  must  be  multiplied  by  \/3  to  obtain  equivalent  tensile  values. 

Tests  were  performed  with  the  aluminum  gradient  specimen  at  both  room  temperature 
and  at  300  degrees  Celsius.  The  room  temperature  test  was  performed  to  provide  a  large 
variation  in  state  to  compare  with  the  higher  temperature  test  results.  The  elevated 
temperature  specimen  reached  300  degrees  from  room  temperature  in  approximately  500 
seconds  and  was  quenched  from  300  degrees  Celsius  to  room  temperature  in  less  than  10 
seconds.  After  compression,  both  specimens  were  sectioned  along  their  axis  of  symmetry. 
Figures  6.24  and  6.25  show  the  variation  in  microhardness  along  the  longitudinal  axis  of 
both  specimens.  Each  data  point  in  the  figures  represents  the  average  of  8  to  10  readings 
at  each  position.  Although  some  variation  is  evident  in  the  room  temperature  specimen, 
there  was  no  detectable  variation  in  microhardness  in  the  quenched  specimen. 

Compression  specimens  (approximately  7.0  mm  in  diameter  by  7.0  mm  in  length)  were 
machined  coaxial  with  the  axis  of  the  compressed  gradient  specimens  as  indicated  in  Fig¬ 
ures  6.26  and  6.27.  These  figures  indicate  the  variation  in  yield  stress  measured  at  room 
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temperature  from  the  compression  specimens  obtained  from  both  the  room  temperature 
and  300  degree  Celsius  gradient  specimens.  Figure  6.28  shows  the  actual  stress-strain 
data  for  the  five  compression  specimens  machined  from  the  300  degree  gradient  specimen. 
Figure  6.27  also  shows  the  variation  in  the  value  of  the  internal  variable  s  predicted  by 
the  ABAQUS  simulation,  where  the  values  of  s  have  been  multiplied  by  a  constant  (.830) 
such  that  the  model  prediction  and  the  measured  yield  stress  coincide  for  the  rightmost 
data  point.  The  model  prediction  appears  to  represent  qualitatively  the  variation  in  yield 
stress  along  the  length  of  the  specimen  very  well,  duplicating  not  only  the  magnitude  of 
variation  in  stress  but  also  the  axial  position  of  that  variation. 

It  should  be  noted  that  the  constitutive  model  cannot  be  extrapolated  to  room  tem¬ 
perature  given  material  constants  which  were  determined  from  tests  at  higher  homologous 
temperatures.  The  constitutive  model  does  not  predict  the  room  temperature  yield  stress 
variation  from  the  quenched  gradient  specimen,  a  calculation  which  may  be  done  by  mul¬ 
tiplying  the  predicted  value  of  s  by  the  proportionality  constant  c  obtained  by  inserting 
room  temperature  parameters  into  equation  5.2.  Such  a  calculation  overestimates  the 
actual  yield  stress  by  approximately  a  factor  of  2. 

Nevertheless,  the  predicted  variation  in  the  internal  variable  appears  to  indicate  the 
relative  variation  in  plastic  deformation  resistance  very  well. 
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Figure  6.1  Correlation  of  constitutive  model  with  isothermal,  constant  true 
strain  rate  tests  on  Fe  -  2%  Si  at  800  degrees  Celsius. 


Figure  6.2  Correlation  of  constitutive  model  with  isothermal,  constant  tru: 
strain  rate  tests  on  Fe  -  2%  Si  at  900  degrees  Celsius. 
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Figure  6.9  Correlation  of  constitutive  mode)  with  isothermal,  constant  true 
strain  rate  tests  on  Fe  -  2%  Si  with  strain  rate  of  -0.001  sec-1. 


Figure  6.10  Correlation  of  constitutive  model  with  strain  rale  jump  tests 
on  Fe  -  2%  Si  at  700  degrees  Celsius. 
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Figure  6.11  Correlation  of  constitutive  model  with  strain  rate  jump  tests  on 
Fe  -  2%  Si  at  800  degrees  Celsius. 


Figure  6.12  Correlation  of  constitutive  model  with  strain  rate  jump  tests 
on  Fe  -  2%  Si  at  900  degrees  Celsius. 
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Figure  6.13  Correlation  of  eonetitutive  model  with  strain  rate  jump  test 
n  Fe  2%  Si  at  1000  degrees  Celsius. 
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Figure  6.14a  Load-unload-hold-reload  tests  on  Fe  -  2%  Si 


Strain 


Figure  6.14b  Simulation  of  load-unload-liold-reload  tests  on  Fe  -  2 7.  Si. 
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Figure  6.15  Load  hislory  used  on  load  boundary  condition  test 
of  constitutive  model. 


Time  (seconds) 


Figure  6.16  Predicted  and  actual  strain  rates  measured  during  load  boundary 
condition  test  on  Fe  -  Z%  Si  compression  specimen. 
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Figure  6.17  Predicted  and  actual  stress  response  for  double  strain  rate  jump 
lest  on  Fe  -  2%  Si. 


Figure  6.18  Predicted  and  actual  stress  response  for  double  strain  rate  jump 
lest  on  Fe  -  2/5  Si 
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Figure  6.19  Strain  rate  decrement  test  on  Fe  -  2%  Si.  Notice  that  data 
does  not  demonstrate  strain  softening  after  decrement  in 
strain  rate. 
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Figure  6.20  State  gradient  specimen  before  and  after 

testing. 
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Groove  spacing:  0.030 
Groove  depth  :  0.005 


Figure  6.21  Dimensions  of  1100  A1  state  gradient  specimen. 
(All  dimensions  in  inches.) 
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Figure  6.22  ABAQUS  finite  element  mesh  of 
aluminum  state  gradient  specimen.  Figure  includes 
original  mesh  geometry  and  shape  at  end  of  test 
simulation. 
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Figure  6.23  Contours  of  constant  internal  variable 
values  from  ABAQUS  simulation  of  1100  aluminum 
slate  gradient  test.  Values  indicated  in  contour  table 
are  in  terms  of  equivalent  shear  stress. 
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Figure  6.26  Variation  in  room  temperature  yield  stress  along  axis  of 

1100  aluminum  gradient  specimen.  Specimen  compressed  at 
room  temperature. 


Figure  6.27  Variation  in  room  temperature  yield  stress  along  axis  of 

1100  aluminum  gradient  specimen.  Specimen  compressed  at 
300  degrees  Celsius  and  then  quenched. 
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Figure  6.28  Stress/strain  data  from  compression 
specimens  machined  from  quenched  1100  aluminum 
gradient  specimen. 


Chapter  7 

Concluding  Remarks 


The  success  of  this  constitutive  model  in  representing  large  rate-dependent  metal  de¬ 
formations  at  high  temperatures,  indicates  that  a  single  scalar  internal  variable  model 
may  be  adequate  for  analyzing  hot  working  operations.  Such  a  result  is  of  significant 
practical  importance.  Depending  on  the  forms  of  the  material  functions,  a  single  scalar 
model  is  most  likely  easier  to  fit  to  test  data  than  multiple  internal  variable  models.  A 
single  scalar  model  also  provides  a  potential  reduction  in  computational  requirements 
both  by  simplifying  the  updating  of  material  states  from  a  tightly  coupled,  stiff  set  of 

i 

nonlinear  differential  equations  and  by  reducing  the  number  of  variable  values  stored  at 
each  integration  point  in  a  large  finite  element  model. 

The  constitutive  model  developed  in  this  report  has  been  recently  used  both  to  model 
other  metals  and  to  simulate  hot  working  processes.  The  model  has  been  applied  to  a 
face-centered  cubic  metal,  1100  aluminum  [Kim,  et  al.,  1986].  The  correlation  between  the 
model  and  the  1100  aluminum  data  appears  to  be  better  than  the  correlations  presented 
in  this  report  for  the  iron  -  2%  silicon  alloy.  This  model  employing  1100  aluminum 
material  parameters  has  been  used  in  ABAQUS  finite  element  simulations  of  closed  die 
forging  operations  Lush,  Weber  and  Anand  [1987].  The  finite  element  predictions  of  load 
histories  and  material  flow  agrees  well  with  actual  forging  tests. 

A  limitation  of  this  model  derives  from  the  fact  that  anisotropic  material  behavior 
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has  been  ignored,  with  the  consequence  that  errors  may  develop  with  large  deformations. 
Small  scale  anisotropy  associated  with  reversed  loading,  Bauchinger-type  phenomena  was 
reasonably  assumed  to  saturate  after  a  few  percent  strain,  and  large  scale  effects  due  to 
the  development  of  texture  were  ignored  primarily  because  of  the  current  lack  of  simple 
techniques  for  modelling  the  evolution  of  crystalline  orientation.  The  fact  that  the  as¬ 
sumption  of  isotropy  was  enforced  in  this  investigation  through  the  use  of  monotonic  tests 
on  identically-oriented  specimens  does  not  mean  however  that  the  texture  development 
in  the  Fe  -  2%  Si  alloy  is  insignificant.  It  would  be  reasonable  in  further  investigations  to 
investigate  the  accuracy  of  the  constitutive  model  as  a  function  of  variation  of  material 
orientation  after  large  deformations. 

The  model  presented  here  has  been  developed  for  single  phase  metals  under  large 
deformations  where  the  dominant  deformation  processes  are  assumed  to  be  thermally 
activated.  Materials  which  may  possess  a  strong  athermal  component  of  deformation 
resistance,  such  as  precipitation  or  dispersion  strengthened  metals  therefore  may  not  be 
well  represented. 

Several  conclusions  may  be  drawn  about  constitutive  models  for  hot  working  given 
the  research  presented  herein: 

1.  The  dependence  of  strain  rate  upon  stress  at  a  given  constant  structure  is  not  well 
represented  by  a  power  law,  much  less  by  that  power  law  exponent  which  correlates 
the  steady  state  stress  to  the  steady  state  strain  rate.  Furthermore,  given  a  single 
internal  variable  model  for  creep,  where  the  dependence  of  strain  rate  upon  stress 
and  state  is  expressed  in  the  form: 

4-'<0’-  7-‘ 
the  exponent  n  should  not  be  the  steady  state  stress  exponent.  If  we  expect  the 
steady  state  value  of  the  internal  variable  s  to  increase  with  stress,  n  must  be 
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greater  than  the  values  of  3  to  5  normally  associated  with  creep. 


2.  If  power  law  breakdown  is  the  result  of  a  change  in  deformation  mechanism,  that 
mechanism  is  associated  with  the  constant  structure  kinetics  and  not  with  the  evo¬ 
lution  of  the  microstructure.  Power  law  breakdown  therefore  should  be  represented 
in  the  deformation  rate  equation. 

3.  It  may  be  possible  to  neglect  static  recovery,  i.e.,  recovery  in  the  absence  of  load, 
in  simulating  many  hot  working  processes.  The  static  recovery  rate  measured  in 
this  investigation  did  not  affect  the  material  response  over  the  range  of  strain 
rates  and  temperatures  tested.  This  restriction  applies  only  to  recovery  processes; 
recrystallization  if  present  will  most  likely  effect  the  material  response. 


Internal  variable  models  which  postulate  separate  rate  equations  for  deformation  and 
microstructure  evolution  requires  material  data  which  is  generally  not  available.  The 
fundamental  consideration  is  that  one  should  perform  experiments  to  characterize  each 
material  function  unambiguously.  Constant  structure  tests  such  as  the  jump  tests  de¬ 
scribed  in  this  report  are  essential  if  one  wishes  to  evaluate  or  fit  the  deformation  equation 
to  data  which  represent  only  deformation  kinetics.  Similarly,  tests  designed  such  that 
deformation  kinetics  may  either  be  neglected,  as  in  the  load-unload-hold-reload  tests,  or 
compensated  for  in  some  phenomenological  way  (e.g.  using  the  constant  c  in  this  model) 
provide  data  representative  of  the  evolution  of  structure. 

Such  testing  addresses  important  questions  concerning  appropriate  material  func¬ 
tions.  One  obvious  case  in  point  is  the  characterization  of  deformation  rate  equations  like 


those  discussed  m  Section  4.1.1.  Strain  rate,  stress,  and  temperature  jump  tests  provide 
the  capability  of  directly  evaluating  deformation  kinetics.  Such  techniques  are  certainly 
applicable  to  fundamental  investigations  on  single  crystals  as  well  as  polycrystalline  en¬ 
gineering  alloys.  It  is  likely  that  such  experimentation  will  become  more  common  as 
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internal  variable  formulations  gain  acceptance.  Such  experimentation  in  conjunction 
with  increased  understanding  of  the  processes  controlling  the  evolution  of  structure  may 
provide,  in  time,  more  consensus  on  the  form  of  constitutive  models  for  high  temperature 
deformation. 
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APPENDIX  A 

Data  Acquisition  and  Reduction  Programs 
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10  *******  ******  ***********  ***************  ************************ 

20  '  *  * 
30  '*  SAMPL2.  BAS  -  DATA  ACQUISITION  US  I  NS  CHANNELS  0  AND  1  * 

40  ’*  Stuart  Brown  Rev.  1.00  2-23-86  * 

50  ’  ********************************************************************** 
60  ’ 

70  ' - STEP  1  - 

SO  'Load  DASH16.BIN  by  loading  at  92K  (Hex  1700)  outside  basic  workspace 
90  CLEAR,  49152! 

100  DEF  SEG  =  0 

110  SG  =  256  *  PEEK (&H51 1 )  +  PEEK(&H510) 

120  SG  =  SG  +  491521/16 
130  DEF  SEG  =  SG 

140  BLOAD  "DASH16.BIN",  0  'Load  it 


150  DIM  DI07.(8)  ,A7.  (20001)  ’Initialise  data  arrays 

160  DIO’/. (0)  =  &H300  'DASH-16  board  base  address 

170  DI0‘/.(1>  =  2  'Selected  interrupt  level  -for  DASH-16 

180  DIO’/.  (2)  *  3  ’Selected  DMA  level  -for  DASH-16 

190  DASH16  =  0  'Declare  initialise  other  CALL  parameters 

200  FLAG’/.  =  0 

210  MD7.  =  0  'Select  Mode  0  -  initialise  driver 

220  CALL  DASH  16  (MD7.,  DIO’/.  (0)  ,  FLAG*/.)  'do  it 

230  IF  FLAG’/.OO  THEN  PRINT  "INITIALIZATION  ERROR":  STOP  'any  error? 
240  ' 

250  ’ - STEP  2 - 

260  ' 

270  'Set  up  multiplexer  scanning  limits  -for  channels  0  to  1 
280  ’ 

290  MD‘/.=  1  'Mode  1  -  set  scan  limits 


300  DIO’/.  (0>  =  1 
310  DIO*/.  <  1 )  =2 

320  CALL  DASH16  <MD’/.,  D 1 0’/.  ( 0 )  ,  FLAG’/. > 

330  IF  FLAG’/.OO  THEN  PRINT  "Error  in  scan  limits  #  " ;  FLAG’/.  :  STOP 
340  ' 

350  ’ - STEP  3 - 

360  'Prompt  for  sampling  rate 
370  ’ 

380  INPUT  "Enter  sampling  rate  in  Herts  (up  to  1000,  even  number)  :  " ,  SAMRTV. 
390  INPUT  "Enter  number  of  desired  samples  (maximum  10,000)  :  ",NSAM7. 

400  SAMDUR  =  NSAM’/./SAMRT’/. 

410  PRINT 

420  PRINT  "Duration  of  sampling  period  is:  ”, SAMDUR,"  seconds." 

430  PRINT 
440  MD’/.  =  17 

450  DIO’/.  (0)  =  2  :  DI0’/.(1)  =  250000 ! /SAMRT’/. 

460  CALL  DASH16  (MD’/.,  DI0’/.(0)  .FLAG’/.) 

470  IF  FLAG’/.OO  THEN  PRINT"Error  in  scan  rate  set  up  #  " :  FLAG’/.  :  STOP 
480  ’ 

490  ’ - STEP  4  - 

500  PRINT  "Press  any  key  to  start  data  collection" 

510  B*  =  INKEY* 

520  IF  B*=" "  THEN  510 

530  'Do  NSAM’/.  conversions  into  array  A’/.(*)  using  mode  4 

540  DIO’/.  (0)  =  NSAM*/.*2  'Number  of  conversions 

550  DI0’/.(1)  =  VARPTR  (A*/.  (0)  )  'Array  locator 

560  DIO’/.  (2)  =  1  'Trigger  source  -  programmable  timer 

570  MD’/.=4 

580  CALL  DASH  16  < MD*/. ,  DIO’/.  (0)  ,  FuAG’/.) 

590  IF  FLAGV.c. >0  THEN  PRINT  "Error  in  mode  4  #  " ;  FLAG*/.  :  STOP 
600  PRINT  "Data  acquisition  done" 

610  PRINT  "Hit  any  key  to  continue" 
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620 

630 

640 

650 

660 

670 

680 

690 

700 

710 

720 

730 

740 

750 

760 

770 

730 

790 

800 

810 

820 


B*  =  INKEY* 

IF  B*=" "  THEN  620 


' - STEP  5 - 

'Review  even  100th  data  point  voltage  and  store  in  integer  -form 


FOR  17.  =0  TO  NSAM7.*2  STEP  200 
A 1  =  A7.  (150  /204. 7  :  A2  =  A7.  ( I7.+  1 ) /204 . 7 

jy.=r/./2 

PRINT  J7.,A1,A2 
NEXT  17. 

INPUT  "Input  -first  data  point  to  be  saved 
INPUT  "Input  last  data  point  to  be  saved 
INPUT  "Input  name  o-f  -file  to  hold  data 
OPEN  INFILE*  FOR  OUTPUT  AS  #2 
FOR  J’/.  =  FSTPT7.  TO  LSTPT7. 

17.  =  J7.  *  2 

PRINT  #  2 ,  J  7. :  A'/.  <  17. )  :A7.ny.+l) 

NEXT  J7. 

CLOSE  #2 
END 


" .  FSTPT7. 

" ,  LSTPT7. 

" , INFILE* 
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Data  Reduction  of  Raw  Load/Deflection  Data 


The  following  procedure  used  to  obtain  smooth  true  stress  versus  true  strain  data, 
and  hardening  data  from  load/deflection  data  resulting  from  uniaxial  compression  or 
tension  tests.  It  uses  four  programs,  TRUECON,  JUMPPT,  POLYS,  and  HARDEN,  to 
generate  a  generic  data  file  which  may  then  serve  as  a  relatively  “clean”  set  of  uniaxial 
test  data.  The  purpose  of  this  procedure  is  to  provide  data  cleansed  of  any  test  machine 
or  data  acquisition  idiosyncrasies,  providing  a  standard  format  available  for  subsequent 
analysis. 

The  procedure  described  herein  assumes  the  following  format  for  the  raw  data  ob¬ 
tained  from  constant  true  strain  rate  tests: 


First  row:  Number  of  data  points(integer) 

Test  ID  (integer) 

Temperature  in  Celsius 
Strain  rate 

Subsequent  rows:  Data  arranged  in  columns 
Column  Contents 

1  Sample  number 

2  Displacement  in  A/D  integer  units 

3  Load  in  A/D  integer  units 


Entries  on  the  same  row  are  separated  by  commas  or  spaces.  The  programs  treat  the  data 
as  free-formatted.  The  programs  are  structured  such  that  the  output  of  a  particularly 
program  serves  as  the  input  of  the  next  program  in  sequence.  Each  program  produces  an 
output  file  which  consists  of  the  input  file  plus  additional  columns  of  data  which  result 
from  the  analysis  of  that  program.  A  copy  of  each  relevant  program  is  attached  to  the 
end  of  this  appendix. 
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A.l  Converting  to  True  Stress/True  Strain  Data:  TRUECON 


The  program  TRUECON  uses  the  raw  data  file  to  produce  a  file  with  the  follow 
structure  in  free-formatted  form. 


First  row: 
Subsequent  rows: 

Column 

1 

2 

3 

4 

5 


Same  as  raw  data  file 
Data  arranged  in  columns 
Contents 
Sample  number 
Load  in  Newtons 
Displacement  in  mm 
True  stress  in  MPa 
True  strain 


TRUECON  requests  the  following  input: 

1.  Specimen  area  in  square  millimeters. 

2.  Temperature  in  degrees  Celsius. 

3.  Coefficient  of  thermal  expansion. 

4.  Compliance  of  test  machine  in  kiloNewtons/mm. 

5.  Starting  actuator  displacement  (zero  load  point). 

6.  File  name  containing  input  data. 

7.  File  name  to  contain  output  data. 


A. 2  Smoothing  of  True  Stress/Strain  Data:  (JUMPPT  and  POLYS) 

Programs  JUMPPT  and  POLYS  serve  to  smooth  the  true  stress  data  produced  by 
program  TRUECON.  JUMPPT  smoothes  the  strain  data  and  calculates  the  strain  rate. 
POLYS  smoothes  the  stress  data.  Both  strain  and  stress  are  smoothed  versus  sample 
number  (which  may  be  considered  equivalent  to  time).  The  programs  smooth  by  fitting 
a  cubic  polynomial  to  segments  of  the  stress/strain  data,  and  uses  the  value  predicted 
by  thv.  cubic  at  the  center  of  the  segment  to  be  the  smoothed  value  at  that  point.  The 
segment  is  called  a  “bin”  and  is  moved  point  by  point  along  the  set  of  data,  producing 
a  smoothed  value  at  the  center  point  of  the  bin.  The  user  specifies  1/2  the  number  of 
data  points  in  the  bin.  The  bin  size  should  never  be  less  than  10  and  frequently  may 


159 


have  to  be  larger  to  obtain  adequate  smoothing.  The  1  /2  bin  size  is  therefore  always  5 
or  greater. 

JUMPPT,  POLYS  and  the  next  program  HARDEN  use  double  precision  IMSL  sub¬ 
routines  to  perform  their  relative  tasks.  The  correct  IMSL  library  obviously  must  then 
be  specified  when  linking  either  of  these  two  programs. 

JUMPPT  requests  the  following  input: 


1.  Name  of  input  data  file  (TRUECON  output). 

2.  First  sample  number  in  first  range  of  data  where  strain  rate  is  to  be  determined. 

3.  Last  sample  number  is  first  range  of  data  where  strain  rate  is  to  be  determined. 

4.  First  sample  number  in  second  range  of  data  where  strain  rate  is  to  be  determined. 

5.  Last  sample  number  in  second  range  of  data  where  strain  rate  is  to  be  determined. 
(Two  ranges  are  specified  so  that  strain  rates  before  and  after  a  strain  rate  jump 
may  be  determined  if  desired.) 

6.  Data  samplinb  interval  in  seconds. 

7.  Name  of  output  data  file. 

8.  First  sample  number  in  range  for  strain  smoothing. 

9.  Last  sample  number  in  range  for  strain  smoothing. 

10.  1/2  bin  size  for  smoothing. 


JUMPPT  produces  the  following  output  file: 


First  row: 
Subsequent  rows: 

Column 

1 

2 

3 

4 

5 

6-8 

9 


Same  as  raw  data  file 
Data  arranged  in  columns 
Contents 
Sample  number 
Load  in  Newtons 
Displacement  in  mm 
True  stress  in  MPa 
True  strain  (smoothed) 
Null  data 
Strain  rate 
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POLYS  requests  the  following  input: 


1.  Name  of  input  file  (JUMPPT  output). 

2.  1/2  smoothing  bin  size. 

3.  First  sample  number  in  range  for  stress  smoothing. 

4.  Last  sample  number  in  range  for  stress  smoothing. 

5.  Name  of  output  file. 


POLYS  produces  the  following  output  file: 


First  row: 
Subsequent  rows: 

Column 

1 

2 

3 

4 

5 

6 

7-8 

9 


Same  a s  raw  data  file 
Data  arranged  in  columns 
Contents 
Sample  number 
Load  in  Newtons 
Displacement  in  mm 
True  stress  in  MPa 
True  strain  (smoothed) 
True  stress  (smoothed) 
Null 

Strain  rate 


The  user  should  plot  the  smoothed  and  unsmoothed  data  after  each  smoothing  operation 
to  confirm  that  the  smoothing  does  not  distort  the  actual  data. 

A. 3  Plastic  Hardening  Determination:  (HARDEN) 

The  next  step  required  is  to  determine  the  rate  of  hardening,  i.e.,  the  rate  of  change 
of  stress  with  strain,  for  each  constant  true  strain  rate  test.  This  is  accomplished  with 
program  HARDEN.  HARDEN  produces  a  new  data  file  of  the  original  data  from  POLYS 
with  two  new  columns: 


1.  Column  7:  Plastic  hardening 

2.  Column  8:  Plastic  strain 
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HARDEN  requests  the  following  input: 


1.  File  name  of  input  data  (POLYS  output). 

2.  First  sample  in  range  of  data  for  hardening  calculation. 

3.  Last  sample  in  range  of  data  for  hardening  calculation. 

4.  1/2  smoothing  bin  size. 

5.  File  name  to  contain  output  data. 

6.  Young’s  modulus  in  GPA  at  test  temperature.  (Used  to  calculate  plastic  strain. 


HARDEN  calculates  the  hardening,  or  slope  of  the  stress  versus  plastic  strain  curve, 
by  selecting  a  bin  of  data  (like  POLYS)  fitting  a  cubic  polynomial  to  the  data  in  the 
bin,  then  using  the  derivative  of  the  polynomial  at  the  center  point  of  the  bin  to  be  the 
hardening  rate  at  that  point.  The  bin  moves  through  the  data,  sequentially  fitting  the 
data  to  the  polynomial  and  taking  the  derivative  at  the  center  point.  The  size  of  the 
bin,  in  data  points,  is  specified  by  the  user  each  time  HARDEN  is  run.  A  good  range 
seems  to  be  between  5  and  10,  although  this  depends  on  the  smoothness  of  the  data 
being  analyzed.  The  bin  size  may  have  to  be  increased  if  there  are  many  data  points 
>  100  which  are  closely  packed  or  if  the  data  is  particularly  noisy. 

To  summarize,  HARDEN  produces  the  following  output  file: 


First  row: 
Subsequent  rows: 

Column 

1 

2 

3 

4 

5 

6 

7 

8 
9 


Same  as  raw  data  file 
Data  arranged  in  columns 
Contents 
Sample  number 
Load  in  Newtons 
Displacement  in  mm 
True  stress  in  MPa 
True  strain  (smoothed) 
True  stress  (smoothed) 
Plastic  strain 
Plastic  hardening 
Strain  rate 


Once  again,  it  is  a  good  idea  to  plot  the  hardening  (column  8)  versus  stress  to  confirm 
that  HARDEN  produced  reasonable  results.  It  is  also  a  good  idea  to  try  several  different 
bin  sizes  to  examine  the  variation  in  predicted  hardening. 
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w*non  n  noonnnnnonnnonoonnnonoonnooonoonoonnooonon 


TRUECDN 


CONVERTS  LOAD  VERSUS  DISPLACEMENT  DATA  FILES 
INTO  TRUE  STRESS  VERSUS  TRUE  STRAIN  DATA  FILES. 

DATA  IS  READ  FRCM  FORTRAN  DATA  FILE  SPECIFIED  BY 
ASSIGN  FOR025  [  INFILE]  STATEMENT.  DATA  IS  OUTPUT 
TO  FOR026  [OUTFILE] .  THE  INPUT  FILE 


FORMAT  IS: 

NLM 

NUMBER  OF  DATA  POINTS 

ARRAY 

SAMPLE,  LOAD.  DISPLACEMENT 

VARIABLES: 

AREA 

SPECIMEN  AREA  IN  MM* *2 

TEMP 

TEST  TEMPERATURE  IN  CELSIUS 

ALPHA 

OOEF  OF  THERMAL  EXPANSION 

CCMP 

TEST  MACHINE  COMPLIANCE 

DISP 

STARTING  DISPLACEMENT 

LI 

INITIAL  SPECIMEN  LENGTH 

LS 

SPECIMEN  LENGTH  ADJUSTED 

FOR  TEMPERATURE 

LQADCF 

CONVERSION  FACTOR  FOR  LOAD, 

(  NEWTONS  AULT  ) 

DISPCF 

CONVERSION  FACTOR  FOR  DISPLACE¬ 
MENT  (MM/VOLT) 

NTESTID 

IDENTIFICATION  NLMBER  FOR 

TEST  (INTEGER) 

TEMPT 

TEST  TEMPERATURE  IN  CELSIUS 

RATE 

TEST  STRAIN  RATE 

NX 

NLMBER  OF  DATA  POINTS 

INFILE 

INPUT  FILE  NAME 

OUTFILE 

OUTPUT  FILE  NAME 

ARRAY(I.J)  DATA  FILE 

SIZE 

TEMPERATURE  SCALING  FACTOR 

FOR  SPECIMEN  DIMENSIONS 

STUART  BROvW  1-3-84 


INITIALIZE  DATA 
REAL  ARRAY2(2000,5),ARRAY(l,3) 

REAL  AREA, TEMP, ALPHA, OCMP, DISP, LI  ,LS, TEMPT, RATE, SIZE, LQADCF 
REAL  DISPCF 

CHARACTER«20  INFILE .OUTFILE 


READ  DATA  AND  CONVERT  TO  TRUE  STRESS  AND  STRAIN 


FORMAT(A20) 
FORMAT  (F10.S) 


PRINT 
READ  ( 
PRINT 
READ  ( 
PRINT 
READ  ( 

print' 

READ  ( 
PRINT 
READ  ( 
PRINT 


."INPUT  SPECIMEN  AREA  IN  NM*NM 
,5)  AREA 

."INPUT  TEMP  IN  DEGREES  CELSIUS 
,5)  TEMP 

."INPUT  COEF  THERMAL  EXPANSION  * 
,5)  ALPHA 

!  "INPUT*  OCMPLIANCE  IN  KNAM 
,5)  OCMP 

."INPUT  STARTING  DISPLACEMENT 
,5)  DISP 

."INPUT  STARTING  LENGTH  IN  Mvl 


1(T6 
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c 


c 

10 

c 

c 

c 


c 

c 

20 


READ  (*  ,5)  LI 

PRINT  *, "INPUT  RAW  DATA  FILE  NAME 
READ  (*,4)  INFILE 

PRINT  *,  "INPUT  LOAD  CONVERSION  FACTOR  ( NEWTONS  /VOLT ) 
READ  (*  5)  LCADCF 

PRINT  •! "INPUT  DISP  CONVERSION  FACTOR  (Mvl/VOLT) 

READ  (*,S)  DISPCF 

PRINT  INPUT  OUTPUT  FILE  NAME 

READ  (*,4)  OUTFILE 

OPEN  (25, FILE-INFILE, STATUS-’OLD’ 

1  ,RECFM-’DS’) 

READ  (25,*)  NIM.TESTID, TEMPT, RATE 
DO  10  I  -  1 ,NLM 


READ  (25,*)  ARRAY( 1,1) ,ARRAY( 1,3) ,ARRAY(1 ,2) 

SIZE  -  1.  +  (TEMP-20.)  *  ALPHA  *  l.E-06 
LS  -  LI  *  SIZE 
SS  -  ARRAY(l ,3) 

ARRAY(1 ,3)  -  DISPCF  «  ARRAY(l,2)  /  204.7 
ARRAY(l ,2)  -  LQADCF  *  SS  /  204.7 

ENG  -  (ARRAY(l ,  3)  -  ARRAY(l  ,2)/(CCMP*1000)  -  DISP)/LS 
ARRAY2( I ,5)  -  -ALOG(l.+  ENG) 

ARRAY2(I,4)  -  -ARRAY( 1,2)  *  (l.  +  ENG)/(AREA*SIZE**2) 
ARRAY2(I ,1)  -  ARRAY(l.l) 

ARRAY2(I ,3)  -  ARRAY( 1,3) 

ARRAY2( 1 ,2)  -  ARRAY(1 ,2) 


CONTINUE 


STORE  TRUE  STRESS  AND  STRAIN  DATA 

OPEN  (  2  6,  FILE-WIT  I LE,  STATUS-' NEW* , 
1  RECFM-'DS’) 

VRITE  (26,*)  NIM.NTESTID, TEMPT, RATE 
DO  20  I-l.NLM 

WilTE  (26,*)  (ARRAY2( I , J) , J-l ,5) 

CONTINUE 
CLOSE  (UNIT-25) 

CLOSE  (UNIT-26) 

STOP 

END 
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ooo  oui  non  ooooooooonooooonnooo 


JIMPPT 


STU  BROW  2-16-86 

DETERMINES  SM30THED  STRAIN  AND  INTERSECTION  STRAIN 


VARIABLES 


ARRAYT(1200,9) 

1 

2 

3 

4 

5 

6 

7 

8 
9 


SAMPLE  NIMBER 
LOAD 

DISPLACEMENT 
TRUE  STRESS 
TRUE  STRAIN 
SMOOTHED  STRESS 
PLASTIC  STRAIN 
HARDENING 
STRAIN  RATE 


DOUBLE  PRECISION  LPT1  .LPT2  ,UPT1  ,UPT2 ,ARRAYT(1200 ,9)  .SAMRATE 
DOUBLE  PRECISION  ARRAYl(l200,2)  ,ARRAY2(1200,2) ,  JLMPPT 
DOUBLE  PRECISION  ALBAP(3) ,DES(5)  ,ANOVA(l4) ,STAT(9) ,PRED(l ,7) 
CHARACTER‘20  INFILE.OUTFILE.ANSWER 


SMOOTHING  VARIABLES 


DOUBLE  PRECISION  SPT1  ,SPT2  ,W)RK(2 , 1200)  .X(100)  ,Y(  100)  ,RSQ,P(200) 
DOUBLE  PRECISION  T(l6) ,C(6) ,S(6)  ,A(3) ,B(3) ,Z(4) 

REAL  RATE(l200) ,STRAIN(1200) ,STRAINL(1200) 

F0RMAT(15A) 


PRINT 
PRINT 
READ  ( 
PRINT 
PRINT 
READ  ( 
PRINT 
READ  ( 
PRINT 
READ  ( 
PRINT 
READ  ( 
PRINT 
READ  ( 


,  ’’PROGRAM  JIMPPT” 

,  "ENTER  NAME  OF  INPUT  DATA  FILE 
,5)  INFILE 

."STRAIN  RATE  DETERMINATION  SECTION 
.’’INPUT  FIRST  POINT  FOR  LOWR  RANGE 
.*)  LPT1 

.’’INPUT  SECOND  POINT  FOR  LOWER  RANGE 
,*)  LPT2 

."INPUT  FIRST  POINT  FOR  UPPER  RANGE 
,*)  UPT1 

."INPUT  SECOND  POINT  FOR  UPPER  RANGE 
,*)  UPT2 

INPUT  DATA  SAMPLING  INTERVAL  (SEC) 
,*)  SAMRATE 


w 


ft 


READ  DATA 


NPT1  -  0 
NPT2  -  0 

OPEN(24,  FILE- INFILE,  RECFM-’DS’) 

OPEN( 12 .FILE- ’ERROR. LS ’ ) 

READ  (24,*)  NPT 
DO  10  I  -  1,  NPT 

READ(24,*)  (ARRAYT( I , J) , J-l ,5) 

ARRAYTd  ,6)  -  0. 

ARRAYT( 1,7)  -  0. 

ARRAYTd, 8)  -  0. 

ARRAYTd, 9)  -  0. 

IF  ( (ARRAYTd,  1  ).GE.LPT1). AND. (ARRAYTd,  1).LE. LPT2))  THEN 


165 


nnnnno 


c 


c 

c 

10 

c 

c 

c 


c 


c 


NPT1  -  NPT1  +  1 
ARRAYl(NPTl.l)  -  ARRAYTC  1,1) 

ARRAYl(NPTl ,2)  -  ARRAYTCl ,5) 

ELSE  IF  ((ARRAYTCl,  1)  .GE.UPT1)  .AND.  (ARRAYTCl,  l).LE.UPT2)) 

1  THEN 

NPT2  -  NPT2  +  1 

ARRAY2 (NPT2 , 1 )  -  ARRAYTCl, 1) 

ARRAY2(NPT2,2)  -  ARRAYTCl, 5) 

END  IF 

CONTINUE 

CLOSEC24) 

PERFORM  LINEAR  REGRESSION 

IX  -  1200 
I  MX)  -  0 
IPRED  -  0 
ALBAP(l)  -  .05 
ALBAPC2)  -  .05 
IP  -  1 
NN  -  0 

CALL  RL0NECARRAY1 ,  IX.NPT1 , 1  MX),  IPRED,  ALBAP,  DES  .ANOVA.STAT, 

1  PRED, IP.NN, IER) 

PRINT  «, "REGRESSION  COEF  FOR  LOWER  RANGE  -  " , (STAT(l)/SAMRATE) 

A1  -  STAT(l) 

B1  -  STATC5) 

PRINT  * , "IER  (ERROR  PARAMETER)  -  ",IER 

IX  -  1200 
IMX)  -  0 
IPRED  -  0 
ALBAP(l)  -  .05 
ALBAPC2)  -  .05 
IP  -  1 
NN  -  0 

CALL  RLONECARRAY2,  IX.NPT2,  IMDD,  IPRED.ALBAP, MS  .ANOVA.STAT, 

2  PRED, IP, NN, IER) 

PRINT  *, "REGRESS ION  COEF  FOR  UPPER  RANGE  -  ”,(STAT(l)/SAMRATE) 
A2  -  STATC 1 ) 

B2  -  STATC 5) 

PRINT  * ,  "IER  (ERROR  PARAMETER)  -  ",IER 


CALCULATE  INTERSECTION 


BEGIN  SMOOTHING  SECTION 

PRINT  • , "  ” 

PRINT  *, "SMOOTH  STRAIN  AND  DETERMINE  STRAIN  RATE?  (Y  OR  N)” 

READ  (*  ,5)  ANSWER 

IF  (ANSWER. EQ. ’N’ )  GO  TO  500 

PRINT  *, "ENTER  NAME  OF  OUTPUT  FILE 

READ  ( * , 5)  OUTFILE 

OPENC25 , FILE-OUTFILE .RECFM- ’DS  ’  ) 

rivxni  '  t  xliiiilA.  im  anivxrx-ii  m-iviuiixv 

READ  (*,*)  SPT1 

PRINT  *, "ENTER  FINAL  SAMPLE  NUMBER 
t>n*r>  ( *  * ")  <;pt9 

PRINT  *\ "ENTER  1/2  BIN  SIZE 
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c 

c 

c 


c 

c 


c 


15 

C 

c 

c 


20 

C 

C 


30 


33 

C 

C 

C 


C 

C 


40 

C 

c 

45 

C 


READ  (*,*)  NBIN 


SELECT'  DATA  ARRAY 


NBEG  -  NPT 

NSAM  -  0 

DO  15  I  -  l.NPT 

IF  ((arraytCi ,l).GE.SPTl).AND. (ARRAYTCI, l).LE.SPT2))  THEN 

NSAM  -  NSAM  +  1 
ARRAY2(NSAM,  1 )  -  ARRAYTC 1 , 1 ) 

ARRAY2(NSAM, 2)  -  ARRAYTC 1, 5) 

NBEG  -  MINCNBEG.I) 

END  IF 
CONTINUE 

SET  UP  WORK  ARRAY 

DO  20  I  -  l.NBIN 

WORKC 1,1)  -  ARRAYT(NBEG-NBIN+I-1,1) 

W0RKC2, I)  -  ARRAYT(NBEG-NBIN+I-1 ,5) 

W0RK(2  .NBIN+NSAM+-I )  -  2*ARRAY2(NSAM,2)-ARRAY2(NSAM-1 ,2) 
WORK(  1  ,NBIN+NSAMf I )  -  2*ARRAY2(NSAM, l)-ARRAY2(NSAM-I , l) 
CONTINUE 

DO  30  I  -  l.NSAM 

W)RK( 1 ,NBIN+I )  -ARRAY2(I,1) 

W»K(2,NBIN+I)  -  ARRAY2C 1 ,2) 

CONTINUE 

0PEN(27,FILE-,ZO7r*  .RECFM-’DS’) 

DO  33  I  -  1 ,NSAMf2*NBIN 

WRITE(27,*)  W)RK(1,I),W)RK(2,I) 

CONTINUE 


CURVE  FIT 


RSQ  -  99.99D0 
MD  -  3 

N  -  2*NBIN  +  1 
DO  50  I  -  l.NSAM 
VKITEC27,*)  I 

DO  40  J  -  l.N 

X(J)  -  W»K(1,I+J-1) 

Y(J)  -  WDRKC2.I+J-1) 

VrRITE(27 , *)  X(J),Y(J) 

CONTINUE 

DO  45  K  -  1,6 

C(K)  -  0.0 
CONTINUE 

CALL  RLFCTIH(X,Y,N,RSQ,MD,  ID.P.C.S  ,A,B,  J.ER) 

CALL  RUDOFMCC,  ID.A.B.T) 

STRAIN(I)  -  SNGL(CCl)  +  C(2)*ARRAY2(I .1)  + 

C(  3 ) *ARRAY2  ( 1 ,1)**2  +  CC 4 ) *ARRAY2 C I ,l)**3) 
RATE(I)  -  SNGL((C(2)  +  2*C(3)*ARRAY2( I ,1)  + 
3*C(4)*ARRAY2(1 ,1)**2)/SAMRATE) 
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no  noon  non 


STRAI NL( I ) -  SNGL( STAT( 1 ) *ARRAY2 (1,1 )+STAT( 5 ) ) 

50  CONTINUE 

CLOSE(27) 

OUTPUT  RESULTS 

VRITE(25,»)  NPT 
DO  60  I-  l.NPT 

IF(  (ARRAYTC  1,1)  .LT.  SPT1 )  .OR.  (ARRAYT(  1,1)  .GT.SPT2))  THEN 
WUTEC2S.*)  (ARRAYTC I,J),J— 1,9) 

ELSE 

VRITE( 25 , * )  (ARRAYTC I , J ) , J-l , 4) , STRAINC I -NBEG+l) , 

1  STRAI NL( I -NBEG+l ) , (ARRAYTC I , J ) , J-7 , 8) ,RATE( I -NBEG+l ) 

END  IF 


60  CONTINUE 

CLOSE(2S) 
500  STOP 

END 
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nnnw  non  nnnnnonnnnnnnoonnnnnnnnonnnnnnnnnnnn 


POLYS 


DATA  SMOOTHING  PROGRAM  FOR  STRESS/STRAIN  CURVES 
STUART  BROAN  5-29-85/MOD  MEHRDAD  HAGHI  7-16-85 
VARIABLES 


NX 

ADA.TA(I  ,NX) 
RCW  1 
RCW  2 
RCW  3 
RCW  4 
RCW  5 
RCW  6 
RCW  7 
RCW  8 
RCW  9 
NPT 

W0RK(2,NX) 
X(  I ) 

Y(I) 

P(I) 

C(6) 

S(6) 

A(3) 

B(3) 


NUMBER  OF  DATA  POINTS 
ORIGINAL  DATA  ARRAY 
SAMPLE  NLMBER 
LOAD 

DISPLACEMENT 
TRUE  STRESS 
TRUE  STRAIN 
SMOOTHED  STRESS 
PLASTIC  STRAIN 
HARDENING 
STRAIN  RATE 

1/2  SMOOTHING  FRAME  SIZE 
YORK  ARRAY  FOR  SMOOTHING 
X  AXIS  TORK  ARRAY 
SMOOTHED  VARIABLE  WORK  ARRAY 
WORK  VT-nrOR 

REGRESSION  COEFFICIENTS 
SCALING  COEFFICIENTS 

n  »* 

n  n 


I/O  ASSIGNMENTS 

FQR025  INPUT  FILE,  PER  ADATA( I , NX)  ABOVE 
FOR026  OUTPUT  FILE,  PER  ADATA(I.NX)  ABOVE 


PROGRAM  POLYS 


INITIALIZE  DATA 


REAL  ADATA(9 ,2000)  ,W0RK(2 ,2000)  ,SSTRESS(2000)  .ZERO 
CHARACTER*20  INFILE.OUTFILE 

DOUBLE  PRECISION  P(lOO) ,T(16) ,X(50) ,Y(50) ,C(6) ,S(6) ,A(3) ,B(3) , 
1  Z(4),RSQ 

RSQ  -  99.99D0 
MD  -  2 
FORMATC15A) 


READ  DATA 


PRINT  * , "VHAT  IS  INPUT  FILENAME? 

READ  (*  ,5)  INFILE 

PRINT  *, "INPUT  1/2  BIN  SIZE 

READ  (*,*)  NPT 

PRINT  * ,  ”VHAT  IS  INITIAL  SAMPLE  NIMBER? 
READ  (*,*)  NSAM1 

PRINT  *,"YHAT  IS  FINAL  SAMPLE  NIMBER?  " 
READ  (*,0  NSAM2 

PRINT  *,"v*IAT  IS  GCarUr  FILtNAMn?  “ 

READ  (* ,5)  OUTFILE 

OPEN  (12,FILE-"ERROR.LS") 

OPEN  (25 .FILE- INFILE, RECFM* ’DS ’ ) 

READ  (25,*)  NX 
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NSAM  -  0 
DO  10  I  -  1, NX 

?}FAn  *)  (ATlATAf  J  I)  9) 

IF  ((ADATA(  1,1)  .GE.NSAMl )  .AND.  (ADATA(  1,1)  .LE.NSAM2) )  THEN 
NSAM  -  NSAM  +  1 
END  IF 

IF  (ADVTA(l.I).EQ.NSAMl)  NBEG-I 


C 

10 

C 

C 

C 

C 


C 


20 

C 

c 


30 

c 

c 

c 

c 


c 

c 


40 

C 


50 

C 

C 

C 

C 


C 

C 

C 


CONTINUE 
CLOSE  (UNIT-25) 


SETUP  MORK  ARRAY 

NEND  -  NBEG  +  NSAM  -  1 
DO  20  I  -  l.NPT 

MORK(l.I)  -  2 *ADATA(  1 , NBBG ) -ADATA(  1 , NBBG+NPT- 1  + 1 ) 

MORK(2,I)  -  2*AEATA( 4 ,NBEG) -AEATA(4 .NBEG+NPT- 1+1 ) 
VORK(2,NSAM+NPT+I )  -  2 *ADATA( 4 , NEND) -  AQATA( 4 , NEND- 1 ) 

MORK(  1 ,  NSAMf NPT+ 1 )  -  2  *AEATA(  1 ,  NEND)  -ADATA(  1 ,  NEND- 1 ) 

CONTINUE 

DO  30  I  -  l.NSAM 

MORKCl.I+NPT)  -  ADATA(  1 , 1 +NBEG- 1 ) 

WDRK( 2 , 1 +NPT )  -  AEWTA(4, 1+NBBG-1) 

CONTINUE 


CURVE  FIT 

N  -  2*NPT  +  1 
DO  50  I  -  l.NSAM 

DO  40  J  -  1,  2*NPT+1 

X(J)  -  DBLE(WRK(l,I+J-l)) 

Y(J)  -  DBLE(WRK(2,I+J-l)) 

CONTINUE 

CALL  RLFOTHCX.Y.N.RSQ.MD.ID.P.C.S.A.B.IER) 

CALL  RLDOPM(C,ID,A.B,T) 

SSTRESSO  )  -  SNGL(C(  1  )+C(2)*ADATA(  1 , 1+NBEG-1 ) 

1  +C(3)*ADATA(1 , I+NBBG-l)**2) 

CONTINUE 


OUTPUT  RESULTS 

OPEN  (26.FILE-OUTFILE) 

VrRITE  (26,*)  NX 
ZERO  -  0.0 
DO  60  I  -  1, NX 


I F  ( (AQATA(  !  ,  I ) .  LT .  NSAMl ) .  OR .  (ADATA(  l.D.GT.  NSA\2  )  )  THEN 


V*ITE  (26,*)  (AOATA(  J , I ) , J-l , 5) .ZERO, (ADATA(K, I ) ,K«7 ,9) 


ELSE 
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U'OU 


c 


WRITE  (26,*)  (ADATA(  J  ,  I )  ,  J-l , 5)  ,SSTRESS( I -NBEG+1) , 
1  (ADATA(K,  I  )  ,K-7 ,9) 

END  IF 

CONTINUE 

CLOSE  (26) 

STOP 

END 
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c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


5 

c 

c 


HARDEN 

HARDENING  PROGRAM  FOR  JIMP  TESTS 
STUART  BROWN  7-18-86 
VARIABLES 
NX 

adataCi  ,nx) 

ROW  1 
R 05V  2 
ROW  3 
RCW  4 
ROW  5 
ROW  6 
ROW  7 
ROW  8 
ROW  9 
NPT 

WORK(  2 ,  NX) 

X(  I ) 

Y(  I ) 

P(I) 

C(6) 

S  (  6  ) 

A(3) 

B(3) 


I/O  ASSIGMvENTS 

FOR025  INPUT  FILE,  PER  ADATA(I.NX)  ABOVE 
FOR026  OUTPUT  FILE,  PER  ADATA(  I  ,NX)  ABOVE 

PROGRAM  POLYS 

INITIALIZE  DATA 

REAL  ADATAC  9 , 2000  )  ,WORK(  2,2000)  ,HARDEN( 2000 ) ,  ZERO ,  E 
CHARACTER* 20  INFILE, OUTFILE 

DOUBLE  PRECISION  P(100) ,T(16) ,X(S0) ,Y(S0) ,C(6) ,S(6) ,A(3) ,B(3) , 

1  Z(4),RSQ 

RSQ  ~  99.S9D0 
MD  -  2 
F0RMAT(15A) 

READ  DATA 


NUMBER  OF  DATA  POINTS 
ORIGINAL  DATA  ARRAY 
SAMPLE  NUMBER 
LOAD 

DISPLACEMENT 

STRESS 

STRAIN 

SMOOTHED  STRESS 
PLASTIC  STRAIN 
HARDENING 
STRAIN  RATE 

1/2  SMOOTHING  FRAME  SIZE 
WORK  ARRAY  FOR  SMOOTHING 
X  AXIS  WORK  ARRAY 
SMOOTHED  VARIABLE  WORK  ARRAY 
WORK  VECTOR 

REGRESSION  COEFFICIENTS 
SCALING  COEFFICIENTS 

H  ft 


PRINT  *,*V«AT  IS  INPUT  FILENAME? 

READ  ( * , 5)  INFILE 

PRINT  * , "WHAT  IS  INITIAL  SAMPLE  NUMBER?  ” 
READ  (*,*)  NSAM1 

PRINT  * ,  ”WHAT  IS  FINAL  SAMPLE  NUMBER? 

READ  (*,*)  NSAM2 
PRINT/, "ENTER  1/2  BIN  SIZE 
READ  (s,s)  Ni'i 

PRINT  *,”WHAT  IS  OUTPUT  FILENAME? 

READ  (*  , 5)  OUTFILE 

PRINT  *, "ENTER  YOUNG’S  MXIULUS  IN  GPA 
READ  (*,*)  E 

OPEN  (12,FILE-”ERROR.LS") 
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OPEN  (26,FILE-’ZOUT’ , RECFM- ’ DS ’  ) 

OPEN  (25, FILE-INFILE, RECFM-’DS’) 

READ  (25,*)  NX 
NSAM  -  0 
DO  10  I  -  1 ,NX 
C 

READ  (25,*)  (ADATA(  J  ,  I ) ,  J-l  ,9) 

ADATA(  7,1)  -  ADATA( 5,1)  -  ADATA(6, 1  )/(E*1000. ) 

IF  ( (ADATA(  1,1)  .GE.NSAM1 )  .AND.  (ADATA(  1,1)  .LE.NSAM2) )  THEN 
NSAM  -  NSAM  +  1 
END  IF 

IF  (ADATA(1 , 1 )  .EQ.NSAM1 )  NBEG-I 
C 

C 

10  CONTINUE 

CLOSE  (UNIT-25) 

C 

C 

C  SETUP  WORK  ARRAY 

C 

NEND  -  NBEG  +  NSAM  -  1 
DO  20  I  -  l.NPT 
C 

WORK(l.I)  -  ADATA(7,NBEG+I-NPT-l) 

IF  (WORK(  1 , 1 )  . LT .WORK(  1,1-1))  V£)RK(  1 , 1  )-WORK(  1 , 1  - 1 )+ . 000001 
VORK(2,I)  -  ADATA(6 .NBEG+I -NPT-l) 

WORK(2,NSAMfNPT+l)  -  2 *ADATA( 6 , NEND) -ADATA( 6 , NEND- 1 ) 

WORK(l  ,NSAMfNPT+I)  -  2  * ADATA(  7 ,  NEND)  - ADATA(  7 ,  NEND- 1 ) 

20  CONTINUE 

C 

DO  30  I  -  l.NSAM 
C 

WORK(l ,  I+NPT)  -  ADATA(7, 1+NBEG-l) 

IF  (WORK(l ,  I+NPT)  .LT.W)RK(l ,  I+NPT-1))  WORK(  1 ,  I+NPT)  - 
1  W)RK( 1 , I+NPT- 1 )+ . 000001 

WORK(2, I+NPT)  -  ADATA(6, I+NBEG-1) 

30  CONTINUE 

DO  31  J  -  1  ,NSAM4-2*NPT 

WRITE(26,*)  WDRK(  1 ,  J )  ,WORK( 2 ,  J ) 

31  CONTINUE 
CLOSE(26) 

C 

C 

C  CURVE  FIT 

C 

N  -  2*NPT  +  1 
DO  50  I  -  l.NSAM 
C 

DO  40  J  -  1,  2*NPT+1 
C 

X(J)  -  DBLE(WORK(l,I+J-l)) 

Y(J)  -  DBLE(WORK(2, I+J-l)) 

40  OWTINUE 

C 

CALL  RLFOTO(X,Y,N,RSQ,KD,  ID,P,C,S,A,B,  IER) 

CALL  RLDOPM(C,  ID,A,B,T) 

HARDEN(I)  -  SNGL(C(2)+2*C(3)*ADATA(7, 1+NBEG-l)) 

50  CONTINUE 

C' 

C 

C  OUTPUT  RESULTS 

C 
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O  On  o 


OPEN  (26.FILE-CUTFILE) 

WRITE  (26,*)  NX 
ZERO  -  0.0 
DO  60  I  -  1, NX 
C 

IF  ( (ADATA(  1 , 1 )  .LT .NSAM1 )  .OR.  (ADATA(l ,  I )  .GT.NSAM2))  THEN 
C 

WRITE  (26,*)  (AEATA(J,I),J-l,7),ZERO,AnATA(9,I) 

C 

ELSE 

C 

WRITE  (26,*)  (AEATA(  J ,  I )  ,  J-l , 7)  ,HARDEN(I -NBEG+1 )  ,ADATA(9 , 1 ) 
C 

END  IF 

CONTINUE 

CLOSE  (26) 

STOP 
END 
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APPENDIX  B 


Material  Parameter  Determination  Programs 
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c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


c 


7 

c 

c 

c 

c 

c 


c 

c 

c 

10 


20 

c 

c 

c 

c 


QNON 


FIT  OF  HARDENING  DATA  BETWEEN  SPECIFIED  VALUES  OF 
STRESS 

STUART  BROWN  8-13-85 
VARIABLES 
NX 

ADATA(NX,7) 

R OW  5 
6 
7 

STRESS(NX) 

STRAIN(NX) 

SL 

su 

Q,QD 
CHO 

SIGS.SIGSD 
SO 

SOCALC 

PROGRAM  QNON 
EXTERNAL  FITQ1  ,FITQ2 

REAL  AQATA(450,7)  .MAXSTRAIN 
DOUBLE  PRECISION  Q.CHO.SIGS 

DOUBLE  PRECISION  X(2)  ,F(450)  ,XJAC(450,2)  ,XJTJ(3)  ,XY(450,2) 

1  .WRKC918)  ,PARM(4)  ,STRAIN(450)  ,STRESS(450)  ,HARDEN(450) 

DOUBLE  PRECISION  EPS,  DELTA,  SSQ,SU,SL,  STS  SUM  .SOCALC 
OGM4DN  STRAIN, STRESS.XV.SIGS.Q, CHO 
CHARACTER‘20  FNAME.PNAME 
FORMAT(  20A) 

INITIALIZE  DATA 
ZXSSQ  PARAMETERS 

N  -  1 
M  -  1 
NSIG  -  6 
EPS  -  0.0D0 
DELTA  -  0.0D0 
MAXFN  -  5000 
I  OPT  -  1 
IXJAC  -  450 
N1  -  1 

READ  STRESS  RANGE 

PRINT  10 

FORMAT  ( *  INPUT  LONER  STRESS  LIMIT  FOR  FIT  ' ) 

READ  * ,SL 
VrRITE  (6,20) 

FORMAT  V  INPUT  UPPER  STRESS  LIMIT  FOR  FIT  ’) 

READ  *  ,avJ 


NUMBER  OF  DATA  POINTS 
DATA  ARRAY 

SMXJTHED  STRESS 
HARDENING 
PLASTIC  STRAIN 
SMOOTHED  STRESS 
PLASTIC  STRAIN 

LOWER  STRESS  DATA  LIMIT  FOR  FIT 
UPPER  STRESS  DATA  LIMIT  FOR  FIT 
EXPONENT  IN  HARDENING  LAW 
C  •  HO 
SIGMA  STAR 
SIGMA  ZERO 

SIGMA  ZERO  FRCM  FIRST  DATA  POINT 


READ  DATA  AND  SET  UP 
ZXSSQ  FIT 
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non  oooo  o^o 


PRINT  * ,  "NAME  OF  DATA  FILE 
READ(* ,7)  FNAME 

OPEN  (25,FILE-FNAME,RECFM-’DS’ ) 

READ(25,*)  NX 

PRINT  * , "READING  DATA” 

DO  30  I  -  1, NX 
C 

READ(25,*)  (ADATA( I,J),J-1,7) 

STRESS(I)  -  DBLE(ADATA( 1,5)) 

STRAINCI)  -  DBLE(ADATA(  1,7)) 

HARDEN(I)  -  DBLE(ADATA(  1,6)) 

IF  ( STRESS ( I ).(7T.SIGS)  THEN 
C 

SIGS  -  STRESS(I) 

MAXS TRAIN  -  STRAIN(I) 

END  IF 
C 

30  CONTINUE 

C 

C  SELECT  DATA  FOR  FIT 

C 

DO  33  I  -  1,NX 
C 

IF  ((STRESSO  )  .GT.SL)  .AND.  (STRESSCl )  .LT.SU)  .AND. 

1  (STRAIN(I)  .LT.MAXSTRAIN))  THEN 

C 

XY(N1,1)  -  STRESSCl) 

XY(N1 ,2)  -  HARDEN(I) 

N1  -  N1  +  1 
END  IF 

CONTINUE 

N1  -  N1  -  1 
CLOSE  (UNIT-25) 

PRINT  «, "FINISHED  READING  DATA" 

CALCULATE  SIGS  AND 
TRANSFORM  DATA 

SIGS  -  SIGS  +  1.0D-3 

HARDENING  FIT 

TOUTE  (6,35) 

35  FORMAT  ( *  INPUT  ESTIMATE  OF  CHO  • ) 

READ  *,X(1) 

TftftITE  36) 

36  FORMAT  ( 5  INPUT  ESTIMATE  OF  Q  * ) 

READ  * ,Q 

PRINT  *, "STARTING  FIRST  ZXSSQ" 

CALL  ZXSSQCFITQl ,N1 , M.NSIG.EPS .DELTA, MAXFN, IOPT.PARM.X, 

1  SSQ,F,XJAC,IXJAC,XJTJ,WRK,INFER,IER) 

PRINT  *, "ZXSSQ  FINISHED” 

CHO  -  X(l) 

VRITE(6,40)  Q 

40  FORMAT  (•  Q  -  ’.D10.5) 

V>RITE(6,50)  SIGS 

r a  TV\m/AT  (  •  r  iro  »  a 

ivxurm  V  o  ;Uo  ••  tuiu  •  o  j 

WITEC6.55)  CHO 

55  FORMAT  (’CHO  -  ’.D10.5) 

W*nE(6,60)  IER 
60  FORMAT  (’IER  -  ’,15) 
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c 

c 

c 

c 

c 


c 


100 

c 

c 

c 

c 


75 

80 


90 

C 

C 

C 

C 

C 


C 


120 

c 


SET  UP  PLOT  FILE 
OF  DATA  VERSUS 
FIT  CURVE 

PRINT  * ,  "NAME  OF  HARDENING  PLOT  FILE 
READ  (*,7)  PNAME 
OPEN(UNIT-24 . FILE-PNAME) 

WRITE(24, *)  NX 
DO  100  I  -  1, NX 

HARD  -  REAL(CH0*(l  .DO  -  (STRESSC I )/S  IGS))*  *Q) 

STRESSR  -  REALC STRESSC I)) 

WRITE  (24,*)  STRESSR, ADATA(  1, 6), HARD 
CONTINUE 

CL0SE(24) 

SET  UP  FOR  FIT  FOR  SO 

PRINT  *, "INPUT  ESTIMATE  OF  SIGMA  ZERO 
READ  (*,*)  X(l) 

PRINT  * ,  "ZXSSQ  STARTING” 

CALL  ZXS SQCFITQ2, NX, N.NSIG.EPS, DELTA, MAXFN, IOPT,PARM,X, 

1  SSQ,F,XJAC,IXJAC,XJTJ  .WORK,  INFER,  IER) 

PRINT  *, "ZXSSQ  FINISHED" 

WRITE(6,75)  INFER, IER 

FORMAT  (’  INFER  -  ’,17,’  IER  -’,15) 

WRITE(6,80)  X(l) 

FORMAT  (*  SO  -  ’  ,D10.5) 

SOCALC-S IGS  -  (S IGS  -DBLE(ADATA(  1,1))) 

1  *DEXP(DBLE(ADATA(l  ,3))*CE0/SIGS) 

WRITE(6,90)  SOCALC 
FORMAT  (*  SOCALC  -  ’.D10.5) 

SET  UP  PLOT  FILE 
OF  DATA  VERSUS 
FIT  CURVE 

PRINT  * ,  "NAME  OF  STRESS/STRAIN  PLOT  FILE 
READ  (*,7)  PNAME 
OPEN(UNIT-25. FILE-PNAME) 

WRITE(25 , * )  NX 
DO  120  I  -  1, NX 

STRESSC  -  REAL(SIGS-((SIGS-X(l))**(l-Q)+ 

1  (Q-l)*CHO*STRAIN(I)/(SIGS**Q))**(l/(l-Q))) 

STRAINR  -  REAL ( STRAIN ( I )) 

STRESSR  -  REAL(STRESSU)) 

WRITE  (25,*)  STRAINR, STRESSR, STRESSC 
CONTINUE 

CL0SEC25) 

STOP 

END 
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FITQi 


EVALUATION  SUBROUTINE  FOR  QNON 
STUART  BROW  8-14-85  ' 

SUBROUTINE  FITQ1(X,NX,N,F) 

DOUBLE  PRECISION  X(2) ,F(450) ,STRESS(450) ,STRAIN(450) ,SIGS,Q 
1  ,XY(450,2),CH0 

CCM4DN  STRAIN, STRESS ,XY,SIGS  ,Q,CH0 

PRINT  *,X 

DO  10  I-l.NX 

F(I)  -  XY( I , 2)  -  X(1)*((1.D0-(XY(I,1)/SIGS))**Q) 

CONTINUE 

RETURN 

END 
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FITQ2 


EVALUATION  SUBROUTINE  FOR  QNON 
STUART  BR05VN  8-14-85 
SUBROUTINE  FITQ2<'X,NX,N,F) 

DCXJBLE  PRECISION  X(2) ,F(450) ,STRESS(450) ,STRAIN(450) .SIGS.Q 
1  ,XY(450,2) ,CHO 

CCMvON  STRAIN, STRESS  ,XY,SIGS,Q,CHO 

PRINT  *,X 

DO  10  I-l.NX 

F( I )  -  STRESS(I)  -  SIGS  +  ((SIGS-X(l))*«(l-Q) 

1  +  (Q-1)*CH0*STRAIN(I)/(SIGS**Q))**(1/(1-Q)) 


CONTINUE 

RETURN 

END 
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ISCXDN 


DETERMINES  MATERIAL  CONSTANTS  FOR  CONSTANT 

STATE  DATA  USING  NONLINEAR,  LEAST  SQUARES,  IMSL  ROUTINE 

ZXSSQ.  (DOUBLE  PRECISION  VERSION) 

STUART  BRCMN  4-4-86 


VARIABLES 


DATA(NX,6) 
COLIMN  1 
2 

3 

4 

5 

6 

TEMP(25) 

RATEI(25) 

RATEFC25) 

CURVE (25) 

STRESSF(25) 

STRESS I (25) 

X(3) 

F(2S) 

N 

NSIG 


DATA  FILE 

TEST  ID  NUMBER 
TEMPERATURE 
INITIAL  STRAIN  RATE 
FINAL  STRAIN  RATE 
FINAL  STRESS 
INITIAL  STRESS 
TEMPERATURE  (KELVIN) 

INITIAL  STRAIN  RATE 
FINAL  STRAIN  RATE 
ID  NUMBERS  OF  CURVES 
STRESS  IMMEDIATELY  AFTER  JUMP 
STRESS  IMMEDIATELY  BEFORE  JUMP 
EQUATION  PARAMETERS 
RESIDUE  FRCM  FIT 
NUMBER  OF  PARAMETERS 
SIGNIFICANT  DIGITS  REQUIRED  FOR 
PARAMETERS 


FILE  ASSIGNMENTS 

FOR025  INPUT  FILE 


SUBROUTINES 

FITISO  USED  IN  ZXSSQ 

ZXSSQ  FITTING  ROUTINE 


PROGRAM  ISOOON 


INITIALIZE  DATA 


EXTERNAL  FITISO 

DOUBLE  PRECISION  X(5) ,F(99) ,XJAC(99,5) ,XJTJ(15) ,RATEI(99) 

1  ,W)RK(238)  ,PARM(4)  ,RATEF(99)  ,STRESSF(99)  ,TEMP(99)  ,CURVE(99) 

2  ,STRESSI(99) 

REAL  DATA( 99,6) 

DOUBLE  PRECISION  EPS  ,DELTA,SSQ,Q,A,STXI  ,MP 
CHARACTER*20  INFILE 

OCMvON  RATE  I  ,RATEF,TEMP,STRESSF,Q,A,STXI  ,MP,  STRESS  I 

N  -  2 

NSIG  -  8 

EPS  -  0.0D0 

DELTA  -  O.ODO 

MAXFN  -  5000 

IOPT  -  0 

IXJAC  -  99 


READ  DATA 


FORMAT  (20A) 

PRINT  *, "ENTER  NAME  OF  DATA  FILE: 
READ  (*  ,6)  INFILE 
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OPEN  (25,  FILE- INFILE,  RECFM-’DS*  ) 

OPEN  (12,FILE-’ERR0R.LS’ ) 

READC25,*)  NX 
DO  10  I  -  1. NX 
C 

READC25,*)  (DATA( I ,  J)  ,  J-l  ,6) 

RATEI(I)  -  DBLE(DATA(  1,3)) 

STRESSF(I)  -  DBLE(EATA(I ,5)) 

STRESSI(I)  -  DBLE(DATA(I ,6)) 

TEMPO )  -  DBLE(DATA(  1,2))  +  273.14D+0 
RATEF(I)  -  DBLE(EATA(I,4)) 

10  CONTINUE 

C 

CLOSE  (UNIT-25) 

C 

C  ESTIMATE  CONSTANTS 

C 

PRINT  *,” INPUT  ESTIMATE  OF  A-Q  TERM  :  ” 

READ  (*  *)  X(l) 

PRINT  INPUT  ESTIMATE  OF  M  :  M 

READ  (*,*)  X(2) 

C 

C  PERFORM  FIT 

C 

CALL  ZXSSQ(FITI SO, NX, N.NSIG.EPS , DELTA, MAXFN, IOPT.PARM.X, 

1  SSQ,F,XJAC,  IXJAC.XJTJ  ,W)RK,  INFER,  IER) 

V.RITE  (*,70)  INFER, IER 

70  FORMAT  (’  INFER  -  ’,17,*  IER  -  \I5) 

TOIITE  (*,80)  X(l) 

80  FORMAT  ( ’  A/Q  -  ’.D12.5) 

WYE  (*,81)  X(2) 

81  FORMAT  (’M-  ’.D12.5) 

PRINT  * , (F(J) , J-l ,NX) 

STOP 

END 
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FITISO 


EVALUATION  SUBROUTINE  FOR  ISOOQN 
STUART  BROW  4-4-86 
SUBROUTINE  FITISO(X,NX,N,F) 

DOUBLE  PRECISION  X(5) ,F(99) ,RATEI (99) ,STRESSF(99) ,TEMP(99) 

DOUBLE  PRECISION  ZI  ,ZF,ARGNLWER,AR<Z£NCM,Q,A,STXI  ,MP,RATEF(99)  ,A1 

DOUBLE  PRECISION  STRESS  I (99) ,ARG1 ,ARG2 

COAON  RATE  I  ,RATEF,TEMP  ,STRESSF,Q,A,  STXI  ,MP,  STRESS  I 

PRINT  * ,X(  1 ) ,X(2) 

DO  10  I  -  1, NX 

ARGNLMER  -  X(l)  *  RATEF( I )**X(2) 

ARGDENCM  -  X(l)  *  RATElO  )**X(2) 

ARG1  -  DL0G(  ARGNLMER  +  DSQRT(  1.0D0  +  ARGNLMER  ‘ARGNLMER  )) 

ARG2  -  DL0G(  ARGDENCM  +  DSQRT(  1.0D0  +  ARGDENCM*  ARGDENCM  )) 

F(I)  -  (STRESSF( I )/STRESSI ( I ) )  -  (ARG1/ARG2) 

CONTINUE 

RETURN 

END 
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SATOONS 


DETERMINES  TWO  MATERIAL  CONSTANTS  FOR  SATURATION 
STRESS  DATA  USING  NONLINEAR,  LEAST  SQUARES,  IMSL  ROUTINE 
ZXSSQ.  (DOUBLE  PRECISION  VERSION) 

STUART  BROW  11-22-86 


VARIABLES 


DATA(NX,6) 
COLUMN  1 
2 

3 

4 

5 

6 

TEMP(25) 

RATE(25) 

SIGSTARC25) 

CURVE(25) 

CHO(25) 

X(3) 

F(2S) 

N 

NSIG 


DATA  FILE 

TEST  ID  NUMBER 
TEMPERATURE 
STRAIN  RATE 
C  *  HO 
SIGMA  STAR 

7PPO 

TEMPERATURE  (KELVIN) 

TEST  STRAIN  RATES 
TEST  SIGMA  STARS 
ID  NUMBERS  OF  CURVES 
C  *  HO  VALUES 
EQUATION  PARAMETERS 
RESIDUE  FRCM  FIT 
NUMBER  OF  PARAMETERS 
SIGNIFICANT  DIGITS  REQUIRED  FOR 
PARAMETERS 


FILE  ASSIGNMENTS 


FOR025  INPUT  FILE 


SUBROUTINES 


FITSAT  USED  IN  ZXSSQ 

ZXSSQ  FITTING  ROUTINE 


PROGRAM  SATCONS 


INITIALIZE  DATA 


EXTERNAL  FITSAT 

DOUBLE  PRECISION  X(3)  ,F(99)  ,XJAC(99 , 3)  ,XJTJ(10) 

1  ,W*K(219),PARM(4),RATE(99),SIGSTAR(99),TEMP(99),CHO(99) 

REAL  DATA(99  6) 

DOUBLE  PRECISION  EPS  .DELTA, SSQ.Q, XI  ,FM,A 
CHARACTER*20  INFILE 

COMMON  RATE,  SIGSTAR,  TEMP  .QD.Q,  XI  ,FM,A 
N  -  3 
NSIG  -  6 
EPS  -  0.0D0 
DELTA  -  0.0D0 
'fMAXFN  -  5000 
I OPT  -  0 
IXJAC  -  99 


READ  DATA 


FO«iviAT(20A) 

PRINT  * ,  "ENTER  NAME  OF  DATA  FILE: 
READ  (* ,6)  INFILE 
OPEN  ( 25, FILE- INFILE, RECFMk ’DS ' ) 
OPEN  (12 ,FILE- ’ERROR. LS ’ ) 
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READC25,*)  NX 
DO  10  I  -  1 ,NX 

READ(25,»)  (DATA(I , J ) , J-l ,6) 

RATE(I)  -  DBLE(DATA(I ,3)) 

SIGSTAR(I)  -  DBLE(DATA(  1,5)) 

TEMP(I)  -  DBLE(DATA( 1,2))  +  273.14D+0 
CHO(I)  -  DBLE(DATA(  1,4)) 

CONTINUE 


CLOSE  (UNIT-25) 

PRINT  * 

, "INPUT  ESTIMATE  OF  STILDA: 

ESTIMATE  CONSTANTS 

READ  (• 

,*)  X(l) 

PRINT  * 

, ’’INPUT  ESTIMATE  OF  N 

•  H 

READ  (* 

,*)  X(2) 

PRINT  * 

,”  INPUT  ESTIMATE  OF  QDS 

•  « 

READ  (* 

,*)  X(3) 

PRINT  * 

, "INPUT  VALUE  OF  XI 

M 

READ  (* 

,*)  XI 

PRINT  » 

."INPUT  VALUE  OF  A 

tf 

READ  (» 

,*)  A 

PRINT  * 

,”  INPUT  VALUE  OF  M 

H 

READ  (• 

,*)  PM 

PRINT  » 

."INPUT  VALUE  OF  Q 

M 

READ  (* 

.*)  Q 

PERFORM  FIT 

CALL  ZXSSQ(FITSAT,NX, N.NSIG.EPS, DELTA, MAXFN,  IOPT.PARM.X, 
1  ,  SSQ.F.XJAC,  IXJAC.XJTJ  ,VORK,  INFER,  IER) 

VRITE  (*,70)  INFER, IER 
70  FORMAT  (’  INFER  -  ’,17,’  IER  -  ’,15) 

VRITE  (*  ,  80)  (X(I), 1-1,3) 

80  FORMAT  (’  STILDA  -  ’.D12.5,/’  N  -  \D12.4, 

1  /’  QDS  ’ ,D12 .5) 

PRINT  * , (FCJ ) , J-l ,NX) 

STOP 

END 
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FITSAT 


EVALUATION  SUBROUTINE  FOR  SATCDNS 
STUART  BROW  11-22-86 
SUBROUTINE  F1TSAT(X,NX,N,  F) 

DOUBLE  PRECISION  X(3)  ,F(99)  ,RATE(99)  ,SIGSTAR(99)  ,TEMP(99)  ,CHO(99) 
DOUBLE  PRECISION  Q,XI  ,PM,A,ARG,C 
CCMON  RATE, SIGSTAR, TEMP, CHO,Q, XI  ,PM,A 

PRINT  *,X(1),X(2) 

DO  10  I  -  1,NX 

ARG  -  ((RATE(I)/A)*DEXP(Q/(8.314D-3*TEMP(I))))**PM 
C  -  ( 1/XI )  *DL0G(ARG  +  DSQRTC  1+ARG*ARG)  ) 

F(I)  -  DLOG(SIGSTARCI))  -  DLOG(C*XCl)*((RATE(l  )/A)* 

1  DEXP(X(3)/(8.314D-3»TEMP(I))))**X(2)) 

CONTINUE 

RETURN 

END 
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HOCAIC 


DETERMINES  HO  AND  INITIAL  STATE  VALUES  USING  MATERIAL 
CONSTANTS  DETERMINED  FRCM  SATOONS,  AND  ZEROON  PROGRAMS. 

STUART  BRCWN  12-13-85 


VARIABLES 


DATA(NX,6) 
COLUMN  1 
2 

3 

4 

5 

6 

TEMPC99) 
RATEC99) 
SIGSTARC99) 
CURVE(99) 
CHOC  99) 


DATA  FILE 

TEST  ID  NUMBER 
TEMPERATURE 
STRAIN  RATE 
C  *  HO 
SIGMA  STAR 
SIGMA  ZERO 
TEMPERATURE  (KELVIN) 
TEST  STRAIN  RATES 
TEST  SIGMA  STARS 
ID  NUMBERS  OF  CURVES 
C  *  HO  VALUES 


FILE  ASSIGNMENTS 


FOR025  INPUT  FILE 

F0R026  OUTPUT  FILE 


PROGRAM  HOCALC 


INITIALIZE  DATA 


REAL  RATEC99)  ,SIGSTAR(99)  ,TEMP(99)  ,CH0(99)  ,SIGZER0(99) 
REAL  DATA(99.6),X,C(99),S0(99),H0(99) 

REAL  ALPHA.A.Q.M.HOAVE.HOSUM 
CHARACTER‘20  INFILE, OUTFILE 

INPUT  CONSTANTS 

PRINT  * , "INPUT  VALUE  OF  A  :  " 

READ  (*  *)  A 

PRINT  •! "INPUT  VALUE  OF  ALPHA  :  ” 

READ  (*,*)  ALPHA 

PRINT  *, "INPUT  VALUE  OF  Q  s  " 

READ  (*,*)  Q 

PRINT  », "INPUT  VALUE  OF  M  :  ” 

READ  (*,*)  M 

READ  DATA 


F0RMATC20A) 

PRINT  *, "ENTER  NAME  OF  INPUT  DATA  FILE: 
READ  (* ,6)  INFILE 
OPEN  (25, FILE-INFILE .RECFM-’DS’) 

PRINT  *, "ENTER  NAME  OF  OUTPUT  DATA  FILE:  " 
READ  (* ,6)  OUTFILE 
OPEN  (26  .FILE-OUTFILE  .RECFM-  ’DS  ’ ) 

READ(25,«)  NX 

T»T\T wn  f  a  /  *  \ 

vVr.1  xnvxo,  -)  1XA. 

c 

HOSUM  -  0. 

C 

DO  10  I  -  1, NX 
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c 


READ(25 , * )  (DATA( I ,  J) ,  J-l  ,6) 

RATE(I)  -  DATA(I,3) 

SIGSTAR(I)  -  EATA(I,5) 

TEMP(I)  -  DATA(I,2)  +  273.14E+0 
CHO(I)  -  DATA(I,4) 

SIGZERO(I)  -  DATA(  1,6) 

X  -  (RATE( I ) *EXP(Q/8. 314E-3/TEMP( I  ))/A)**M 
C(I)  -  (1 /ALPHA) *AL0G(X+SQRT(1+X*X)) 

S0(I)  -  SIGZEROC I )/C( I ) 

H0( I )  -  CH0( I )/C(  I ) 

HOSLM  -  H0( I)  +  HOSLM 

WIITE  (26,9)  (DATA( I , J ) , J— 1 ,6) ,S0( I ) ,C( I ) ,H0( I ) 
FORMAT(2F9.0,F6.3,4F10.4,F9.S,F9.0) 

CONTINUE 

CLOSE  (UNIT-26) 

CLOSE  (UNIT-25)  ■ 


PERFORM  FIT 


HOAVE  -  HOSIM/NX 
VRITE  (*,70)  HOAVE 
70  FORMAT  (’HO  -  ’.F9.3) 

STOP 
END 
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APPENDIX  C 

Constitutive  Model  Integration  Routines 
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c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


5 

c 


c 

c 

c 


ISOINT 

INTEGRATION  OF  ISOTHERMAL,  CONSTANT  RATE  TESTS 
STUART  BROWN  10-29-86 


VARIABLE 

DESCRIPTION 

CP 

STRESS/STATE  PROPORTIONALITY 

THDLD 

HOLD  TIME 

EMDD 

YOUNG’S  MXJULUS  IN  MPA 

TOUR 

DURATION  OF  TEST 

RATE 

STRAIN  RATE 

N 

NUMBER  OF  EQUATIONS 

IND 

SUBROUTINE  OPTION  FLAG 

TOL 

NW 

LOCAL  CONVERGENCE  TOLERANCE 

X 

TIME 

Y(l) 

C(  21) 

1ER 

W(l,9) 

VALUE  OF  STATE  VARIABLE 

PROGRAM  ISOINT 
EXTERNAL  EVOL 


REAL  STRESS.  STRAIN, THDLD,W(1, 9  ),C(  24)  ,Y(l) 

CHARACTER* 20  OUTFILE,  INFILE 

roScA?1*’  HO,  SSTAR,  P,  A,  SO,  B,  QR,  SMDD 


PRINT  * , "ENTER  INPUT  FILENAME 
READ  (• ,5)  INFILE 
PRINT  *, "ENTER  OUTPUT  FILENAME 
READ  ( * ,  5)  OUTFILE 


OPEN(24,FILE 

1 

*—* 

*— * 

E 

| 

£ 

READ  (24. 

•) 

RATEl 

READ  (24, 

*) 

RATE2 

READ  (24, 

*) 

SJIMP 

READ  (24, 

•) 

TEMP 

READ  (24, 

*) 

SFIN 

READ  (24, 

*) 

NSTEPLD1 

READ  (24, 

«) 

NSTEPLD2 

READ  (24, 

*) 

Y(l) 

READ  (24, 

*) 

EMDD 

READ  (24, 

*) 

SMDD 

READ  (24, 

*) 

XI 

READ  (24, 

•) 

Q 

READ  (24, 

*) 

EM 

READ  (24, 

*) 

AP 

READ  (24, 

*) 

HO 

READ  (24, 

*) 

A 

READ  (24, 

*) 

STILDE 

READ  (24, 

») 

EN 

tJCArV  ( A 

\  At  *T  » 

READ  (24, 

*) 

rrvc 

yuu 

*) 

B 

READ  (24. 

*) 

QR 

READ  (24,*) 
CLOSE  (24) 

P 

READ  DATAFILE 

.STATUS- 'OLD’) 

STRAIN  RATE  ONE 
STRAIN  RATE  TWO 
JIMP  STRAIN 

TEMPERATURE  IN  CELSIUS 
FINAL  VALUE  OF  STRAIN 

#  STEPS  IN  LOADING  PERIOD  1 

*  STEPS  IN  LOADING  PERIOD  2 
INITIAL  VALUE  OF  S  (MPA) 
TOUNGS  MODULUS  (MPA) 

SHEAR  MODULUS  (MPA) 
STRESS/STATE  SCALING 
RATE  ACTIVATION  ENERGY 
RATE  EXPONENT 
PRBEXPONENTIAL 
ATHERMAL  HARDENING  RATE 
EXPONENT  IN  HARDENING  LAW 
SATURATION  SCALING 
SATURATION  EXPONENT 
SATURATION  Aj^x  IVATION  wnukGY 
STATIC  REC  PREMULTIPLIER 
RECOVERY  ACTIVATION  ENERGY 
RECOVERY  TERM  EXPONENT 
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SO  -  Y(l) 

ARG  -  ((RATE1/AP)*EXP(Q/((TEMP+273.16)*8.314E-3)))**EM 
CP  -  ( 1 .  /XI )  *ALOG(ARG+SQRT(  1 .  +ARG*ARG)  ) 

SSTAR-STILDE*  (  (RATE1/AP )  »EXP(QDS/ (  (TEMP+273 . 16  )  *  8 . 314E- 3  )  )  )  *  *EN 

TINT1  -  S JUMP /RATE 1 

TINT2  -  (SFIN-SJLMP)/RATE2 

TLD1  -  TINT1/NSTEPLD1 

TLD2  -  TINT2/NSTEPLD2 

OPENCUNIT-25 ,  FILE-CUTFILE) 

V*RITE  (25,*)  NSTEPLD1  +  NSTEPLD2 

INITIAL  LOADING  INTEGRATION  LOOP 


X  -  0.0 
RATE  -  RATE1 
IDO  100  I  -  1,  NSTEPLD1 
C 

N  -  1 

TEND  -  X  +  TLD1 
IND  -  1 
m  -  i 

TOL  -  .0010 

CALL  DVERK(N, EWL.X.Y, TEND, TOL,  IND, C.NW.W.IER) 

STRESS  -  CP*Y(l) 

STRAIN  -  (  STRESS  /  EMX>  )  +(  X  *  RATE1  ) 

■WRITEC25,*)  STRAIN,  STRESS, TEND,  Y(l) 

C 

IF  ((IND.LT.O)  .CXI.  (IER.GT.O. ))  THEN 

C 

PRINT  *,”IND  -  ",IND,”  AND  IER  -  ",IER 
GO  TO  400 
C 

END  IF 
C 

100  CONTINUE 
C 

X  -  0.0 
RATE  »  RATE2 

ARG  -  ((RATE2/AP)*EXP(Q/((TEMP+273.16)*8.314E-3)))**EM 
CP  -  ( 1 .  /XI )  *ALOG(ARG+SQRT(  1 .  +ARG*ARG)  ) 

SSTAR-STILEC*  (  (RATE2/AP)  *EXP(QDS/ (  (m<P+273 . 16  )  *  8 . 314E- 3  )  )  )  *  *EN 
DO  200  I  -  1 ,  NSTEPLD2 
C 

N  -  1 

TEND  -  X  +  TLD2 
IND  -  1 
m  -  i 
TOL  -  .0010 

CALL  IMDUCCN, EWL,X,Y, TEND, TOL,  IND, C,NW,W,  IER) 

STRESS  -  CP*Y(1) 

STRAIN  -  (  STOESS  /  EMDD  )  +(  X  *  RATE2)  +  SJUMP 
WITEC25,*)  STRAIN,  STRESS, TENDfTINTl.Y(l) 

C 

IF  ((IND.LT.O) .OR. (IER.GT.O. ))  THEN 
C 

i'Kim  *,  inn  -  ,  nxD,  aixu  in*.  -  ,  lov 

GO  TO  400 
C 

END  IF 
C 

200  CONTINUE 
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C 

400  CLOSE(25) 
STOP 
END 


EVOL 

EVALUATION  SUBROUTINE  FOR  ISOINT 

STUART  BRCSVN  10-28-86 

SUBROUTINE  ETOL(N,X,Y, YPRIME) 

REAL  Y(N)  ,YPRIME(N)  .TODLD 

OCMON  RATE,H0,SSTAR,P,A,S0,B,QR,SMX> 

YPRIME(l)  -  HO  •  C ( 1  - (Y(  1 ) /SSTAR) ) * *A)  *  RATE  - 
1  (B*(Y(l)/SMX>)**P)»SM3D*EXP(-QR/((TEMP+273.l6)*8.314E-3)) 


RETURN 

END 


192 


RECINT 

INTEGRATION  OF  RECOVERY  TESTS 
STUART  BR05W  10-29-86 


VARIABLE 

DESCRIPTION 

CP 

STRESS/STATE  PROPORTIONALITY 

THOLD 

HOLD  TIME 

EMX) 

YOUNG’S  MODULUS  IN  MPA 

TOUR 

DURATION  OF  TEST 

RATE 

STRAIN  RATE 

N 

NIMBER  OF  EQUATIONS 

IND 

SUBROUTINE  OPTION  FLAG 

TOL 

LOCAL  CONVERGENCE  TOLERANCE 

m 

X 

TIME 

y(i) 

VALUE  OF  STATE  VARIABLE 

C(2l) 

IER 

Wd,9) 

PROGRAM  RECINT 

EXTERNAL  EVOL 

REAL  STRESS. STRAIN, THDLD.Wd, 9), C(24),Y(1) 

CHARACTER*20  OUTFILE, INFILE 

CCNMON  THOLD,  RATE,  HO,  SSTAR,  B,  EMJ,  TEMP,  QR,  P,  THB.  A,  SO 
FORMAT (A) 

PRINT  •, "ENTER  INPUT  FILENAME 
READ  (* ,5)  INFILE 
PRINT  *, "ENTER  OUTPUT  FILENAME 
READ  (* ,5)  OUTFILE 


READ  DATAFILE 


OPEN(24, FILE-INFILE, RECFM-’DS’  .STATUS-’ OLD’  ) 


READ 

(24, 

*) 

RATE 

I  STRAIN  RATE 

READ 

(24. 

*) 

SHOLD 

1  HOLD  STRAIN 

READ 

(24, 

*) 

THOLD 

!  TIME  HELD  FOR  RECOVERY 

READ 

(24, 

*) 

SPIN 

!  FINAL  VALUE  OF  STRAIN 

READ 

(24, 

*) 

NSTEPLD1 

1  #  STEPS  IN  LOADING  PERIOD 

READ 

(24. 

*) 

NSTEPHD 

i  #  STEPS  IN  HOLD  PERIOD 

READ 

(24, 

*) 

NSTEPLD2 

1  *  STEPS  IN  RELOADING  PERIOD 

READ 

(24, 

•) 

Y(l) 

!  INITIAL  VALUE  OF  S  (MPA) 

READ 

(24, 

*) 

EMCD 

I  YOUNGS  MDDUIUS  (MPA) 

READ 

(24, 

*) 

CP 

!  STATE/ STRESS  PROPORTIONALITY 

READ 

(24, 

•) 

HO 

1  ATHERM4L  HARDENING  RATE 

READ 

(24, 

*) 

A 

I  EXPONENT  IN  HARDENING  LAW 

READ 

(24, 

*) 

SSTAR 

!  SATURATION  VALUE  OF  S 

READ 

(24, 

•) 

B 

I  PREMJLTIPLIER  IN  RECOVERY  TERM 

READ 

(24, 

•) 

EMJ 

I  SHEAR  MDDULUS  IN  GPA 

READ 

(24, 

•) 

QR 

!  RECOVERY  ACTIVATION  ENERGY 

READ 

(24, 

*) 

TEMP 

1  TEMPERATURE  IN  KELVIN 

READ 

(24, 

*) 

P 

1  RECOVERY  TERM  EXPONENT 

CLOSE  (24) 


SO  -  Y(l) 

THB  -  SHDLD/RATE 


c 

c 

c 


c 


c 

c 


c 

c 

100 

c 

c 

c 

c 


c 

c 


c 


TLD1I  -  THB/NSTEPLD1 
TLD2  -  (SFIN-SHOLD)/RATE 
TLD2I  -  TLD2/NSTEPLD2 
THDI  -  THOLD/NSTEPHD 
OPEN (UNIT-25 , FI LE-OUTF I LE) 

ViRITE  (25,*)  NSTEPLD1 +NSTEPHDf NSTEPLD2+ 1 

INITIAL  LOADING  INTEGRATION  LOOP 


X  -  0.0 

DO  100  I  -  1,  NSTEPLD1 
N  -  1 

TEND  -  X  +  TLD1I 
IND  -  1 
m  -  i 
TOL  -  .0010 

CALL  DVERK(N, EVOL.X.Y, TEND, TOL, IND, C,Ntf,W,IER) 
STRESS  -  CP*Y(l) 

STRAIN  -  (  STRESS  /  IM3D  )  +(  X  *  RATE  ) 
WRITE(25, *)  STRAIN, STRESS ,TEND,Y( 1 ) 

IF  ((IND.LT.O).OR.(IER.CT.O.))  THEN 

PRINT  * , "IND  -  " , IND, "  AND  IER  -  ",IER 
GO  TO  400 

END  IF 


CONTINUE 


HOLD  INTEGRATION  LOOP 


DO  200  I  -  1,  NSTEPHD 


N  -  1 

TEND  -  X  +  THDI 
IND  -  1 
m  -  i 
TOL  -  .001 

CALL  DVERK(N, EVOL.X.Y, TEND, TOL,  IND, C.NW.W,  IER) 

STRESS  -  0.0 

STRAIN  -  THB  •  RATE 

WITE(25,*)  STRAIN,  STRESS,  TEND,  Y(l) 

IF  ((lND.LT.O).C».(IER.Crr.O.))  THEN 

PRINT  * , "IND  -  M,IND,"  AND  IER  -  ",IER 
GO  TO  400 

END  IF 


C 

C 


CONTINUE 


RELOADING  INTEGRATION 


STRESS  —  CP  *  Y( 1 ) 

STRAIN  -  (STRESS/EMOD)  +  (X-THDLD) ‘RATE 
■WRITE  (25,*)  STRAIN, STRESS, X,Y(1) 

DO  300  I  -  1,  NSTEPLD2 
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N  -  1 

TEND  -  X  +  TLD2I 
IKD  -  1 

m  -  i 
TOL  -  .001 

CALL  DVERK(N, EWL,X,Y, TEND, TOL,  IND,C,NW»W,  IER) 
STRESS  “  CP  *  Y(l) 

STRAIN  -  (  STRESS  /  EMD  )  +  (X-THDLD)  *  RATE 
WRITE(25, *)  STRAIN,  STRESS, TEND,  Y(l) 

C 

IF  ((IND.LT.O)  .OR.  (IER.GT.O.  ))  THEN 
C 

PRINT  * , ” IND  -  ”, IND,”  AND  IER  -  ", IER 
GO  TO  400 
C 

END  IF 
C 

300  CONTINUE 

400  CLOSE(25) 

STOP 

END 


EVOL 

EVALUATION  SUBROUTINE  FOR  RECINT 

STUART  BROW  10-28-86 

SUBROUTINE  EVOL(N,X,Y,YPRIME) 

C 

REAL  Y(N)  ,YPRIME(N)  ,TH3LD 

OOM4DN  THOLD, RATE, HO, SSTAR,B,EMJ, TEMP, QR,P,THB,A,SO 
C 

IF  ((X.LT.THB) .(XI. (X.GE. OlHB+THOLD)) )  THEN 
C 

YPRIME(l)  -  HO  *  ( ( 1  - (Y(  1 ) /SSTAR) ) * *A)  *  RATE 
1  -  B*EXP(-QR/8.314E-3/TEMP)*EMJ  *  (Y( l)/EMJ)«*P 

C 

ELSE 

YPRIME(l)-  -B*EXP(-QR/8.314E-3/TEMP)«EMJ*(Y(l)/EMLJ)**P 
C 

END  IF 
C 

RETURN 

END 
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SUBROUTINE  IM\T( STRESS ,  STATEV , EDSDDE ,  SSE ,  SPD ,  SCD ,  STRAN , DSTRAN , 

1  TIME, DTIME, TEMP .DTEMP.PREEEF.DPRED.^TERL.NDI , NSHR , NTENS , 

2  NSTATV ,  PROPS ,  NPROPS ,  COORDS ) 

C** *****************  ***************************************  ************ 

C  Isotropic  Thermo-Elasto-Viscoplasticity  with  pressure  sensitive 

C  plastic  flow  and  plastic  dilatancy. 

c,«, ,,.****************•****************.»********************». ******* 

C  This  IMAT  version  interfaces  with  the  *VISCD  procedure  in  ARAQUS. 

C  Automatic  timestep  control  is  done  us ing  the  CEKftX  parameter.  The 
C  timestep  is  decreased  if  CEMAX.  exceeds  CETOL. 

Q*  *********** ********** ******* ****** *********** ************************ 

C  This  IMVT  version  is  not  for  use  in  plane  stress  or  any  other  cases 
C  where  more  strain  terms  than  stress  terms  are  used. 

Qt  ********************************************************************* 

C  State  Variables : 

C  STATEV(l)  -  S  (plastic  flow  resistance, tensile, suggested 

C  units  are  N/m2) 

C  STATEV(2)  -  TH  (temperature , suggested  units  are  Deg.K) 

C  Two  more  are  for  debugging  and  plotting  purposes. 

C  STATEV(3)  -  F  (plastic  shear  strain  rate) 

C  STATEV(4)  -  time  integral  of  F  (plastic  shear  strain) 

0*  ********************************************************************* 

C  Contents  of  PROPS  vector  in  this  version: 

C  J  PROPS(J) 

C . 

C  1  PLSLMT  --  limit  on  equiv.  plastic  tensile  strain  increment 
C  2  PHI  --  degree  of  implicitness  (ranges  from  0  to  1) 

C  3  CMEGA  --  fraction  of  plastic  work  going  into  adiabatic  heating 
C  4  SO  --  initial  value  for  internal  variable  S 

C  5  TO  --  initial  value  for  temperature 

C  6  At AJ  --  shear  modulus 

C  7  /1*APPA  --  bulk  modulus 
C  8  ALPHA  --  thermal  expansion  coefficient 
C  9  RHO  --  density 
C  10  C  - -  specific  heat 

C  11  A  --  pre-exponential  factor  in  plastic  shear  strain  rate 
C  12  Q/R  --  activation  energy  divided  by  universal  gas  constant 
C  13  XP  --  exponent  in  eqn  for  F  (1. /strain  rate  sensitivity) 

C  14  PALPHA  --  pressure  sensitivity  parameter 

C  15  HO  --  pre-multiplier  in  hardening  rate  equation 

C  16  SB  --  pre-multiplier  for  saturation  value  of  S 

C  17  AN  --  exponent  in  equation  for  saturation  value  of  S 

C  18  B  --  pre-multiplier  for  rate  of  static  restoration 

C  19  SA  --  annealed  value  of  S 

C  20  P  --  exponent  in  equation  for  RDOT 

C  The  above  list  of  properties  is  specific  to  the  functions  chosen  for 
C  modeling  one  particular  material.  For  other  materials  a  different  list 
C  may  be  desirable.  For  example,  additional  properties  would  be  needed 
C  to  define  functions  of  temperature  for  the  elastic  constants  and  for 
C  the  annealed  value  of  internal  variable  S. 

C* ********************************************************************* 

C  See  subroutines  UdPROP,  GAMXTT,  and  SDOT  for  suggested  units  for 
C  the  above  properties. 

c,  ***********.***************•******•*****•***•************************ 
C  The  values  input  for  S0,A,H0,SB,B,SA  should  correspond  to  tensile 
C  test  data.  These  values  are  converted  to  equivalent  shear  values 
C  within  the  subroutines. 

c******************************* *******  ******************************** 
C  The  parameter  PHI  controls  the  degree  of  impl icitness  of  the 
C  integration  procedure. 

C  PHI-0.0  -  explicit 

C  PHI-1.0  -  fully  implicit 

C  Suggested  value  of  PHI:  equal  to  or  greater  than  0.50 


196 


C*******t*******************************t****************************** 

C  The  parameter  CMEGA  control*  whether  the  problem  is  isothermal  or 
C  adiabatic. 

C  CMEGA-0.0  - isothermal 

C  CMEGA.- 1 .0  - adiabatic 

C  Suggested  value  for  fast  deformations  is  0.9 

C* ****** t****************************************** ***********  ********* 

IMPLICIT  R£AL*8(A-H,0-Z) 

c«. **...•*.*«*,«, ..•.?«**.••*.•*.•**.•******.,**.*****.***«*, .«*****.«* 

C  Common  blocks  CERROR  and  CONSTS  appear  here  exactly  as  they  exist  in 
C  ABAQUS  version  4-5-159.  They  will  generally  be  different  for  other 
C  ABAQUS  versions. 

C********************************************************************** 

CCMyDN/CERROR/RESMAX(  30  ) ,  JNREMXC  30  )  ,  ERRMAX(  2  )  ,  CETOL ,  CSLIM, 

1  CEMOt,PCIOL,TLIMIT,PSUBIN,RES\ON,IXMAX(30) ,  JNDCMX(30)  .ERRPRE, 

2  UIXiLSS  ,PTOL,AMTOL,EMCET,EMRETL,SIGTOL,DSIGMX,UTOL,IMAX,U4MAX, 

3  \MAX,V4MAX,AM\X,A4\AX,TMAX,EPPM«,RKftX,R4MAX,NGOPEN,NGaJOS, 

4  ROnTTOL.ROTFAC,  JRIKND.NINCCS  ,RIKUB,RIKLMX,RIKMU,RIKLAM,RIKDLA,  • 

5  RIKRO,RIKOLD,RIKLMX,CM\X,IXMAXP ,  STRRAT.PCUT.RIKDLO 

c*. *,«*««»*. ««***•***•**• . ***************************************** 

CCMCN/CDNSTS  /P I ,  S I N60 ,  COS  6  0 ,  KCROS  2(3)  ,KCROS  3(3),  ZERO ,  LZERO ,  LONE , 

1  ONE,  VfJO,  HALF,  ABIG,ASM\LL,BCBIG,L0CSHR(2 , 3),  THIRD,  PRECIS,  BLANK 

£,*«««*. *,*.,*,**«,.****»...«*.. ****** *,.*,**, ******* ****************** 

DIMENSION  STRESS(NTENS)  , STATEV(NSTATV)  , EDSDDE(NTENS , NTENS )  , 

1  STRAN( NTENS  )  ,  DSTRAN  ( NTENS ) ,  PREDEF  ( 1 )  ,  DPRED(  1 ) ,  PROPS  (NPROPS  ) , 

2  OOORDS(3) 

^******«»*s***s*c**«**s****s***********************«*s****************s 

SQART3-m)*SIN60 
NDIP1-NDI+1 
PHIDT-PROPS  (  2  )  *DTIME 

c*  «*»*»•*******»*»**«•*»««***•(***••*«*«*»*•»*»* »***•**•**«*****»*»*»*» 
C  Initialize  the  state  variables,  if  necessary. 

c.»,*.,,t,,»,««,,*,«***.t***s****. **»•»**•*»*»*•«»«***«*. ************** 

IF  ( STATEV( 1 ) . LE . ZERO)  THEN 
STATEV( 1 ) -PROPS ( 4 ) /SQART3 
STATEV( 2 )-PROPS (5) 

STATEV( 4 ) -ZERO 
END  IF 

C* ********************************************************************* 

C  Set  the  state  variables. 

c********************************************************************** 

S-STATEV(l) 

TH-STATEV(  2  ) 

c* ********•*•*«**•****••••••*•••*«*****•*****•*•***•**••••*********•*•* 

C  The  parameter  TAUTW.  is  the  minimum  value  allowed  for  TALJB.  This 
C  avoids  difficulties  that  could  arise  when  dividing  by  TAUB.  The  value 
C  used  corresponds  to  a  stress  level  at  which  F  is  considered 
C  negligible,  while  1 . /TAL7IOL**2  can  still  be  evaluated. 

c********************************************************************** 

TAUTOL-S’l  .D-6 

c*  ****************************** *************************************** 
C  Subroutine  IMPROP  determines  AMJ, AKAPPA, ALPHA, RHD, and  C  based  upon 
C  the  temperature  TH,  using  data  supplied  in  PROPS.  For  the  present 
C  case,  the  properties  are  assumed  constant  and  input  directly  in  PROPS 
C  In  other  cases,  additional  data  constants  defining  functions  of  TH  fo 
C  each  property  may  be  input. 

Ct*****.... ************* *******•••*•«*• t******************************* 

CALL  IMPROP  (AMU, AKAPPA, ALPHA, RHD, C.TH, PROPS, NPROPS) 
c>, ,**,*.«*«*******.•*•*****••*•*********••*, ************************** 

C  Pressure  PB  and  equivalent  shear  stress  TAUB 

c«. ...**«*****.*.»*«.****•*.*.».«*.*.***».*.*«****.**•**«***. ********** 

CALL  SINV  ( STRESS, SINV1.SINV2) 

TAUB-HvAXl(SINV2/SQART3,TAUIOL) 
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PB— SINV1 

C**«***«***t**«**tt«***t***«t«(»t**t************t*t*ttt**tt**t********* 

C  Subroutine  SBETA  determines  the  vslue  of  the  plastic  dilatancy 
C  factor  BETA.  Presently  set  to  zero. 

c*t**«**«.*««*s»**«*********«*«**, a*..*,.,*,,,,*.*,,*,,, ****** 

CALL  SBETA  (TAUB,PB,TH,S .BETA) 

C********************************************************************** 
IF  (PROPS(3) .GT.ZERO)  THEN 
CDN1 -PROPS  (  3  )*  (TAUB  -  BETA*  PB  )/(  RHD*C) 

ELSE 

CONI -ZERO 
END  IF 


Subroutine  GAMXTT  determines  the  equivalent  plastic  shear  strain  rate 
F  and  its  derivatives  PDA,PDB,PDC,PED  with  respect  to  TAUB,PB,TH,and 
S,  respectively.  To  make  the  subsequent  calculations  more  convenient, 
PDB.PDC,  and  PDD  are  returned  as: 

RATIOS  -  PDB/PEA 
RATIOC  -  PDC/PDA 
RATICD  -  PED/PEA 

CON2  is  returned  for  use  in  subroutine  SDOT. 


CALL  GAMX3T  (TADB , PB ,TH, S , PROPS  .NPROPS , SQART3 , 
1  Fl,PEA,RATIGB,RATIOC,RATIOD,<X>N2) 

AMLJB-AMJ*TADB/  (TAUB+AMLJ*PHIET*F1 ) 


Subroutine  SDOT  determines  the  hardening  rate  H  and  the  static 
restoration  rate  RDOT.  SA,  the  annealed  value  of  S  at  this 
temperature,  is  also  returned. 


CALL  SDOTCTAUB ,  PB  ,TO ,  S  ,AMJ , CON2 ,  PROPS , NPROPS ,  SQART3 ,  FI  ,H , RDOT ,  SA) 
DR-RDOT*DTIhE 


G-AMJ-  (  RATIOS  *AKAPPA*BETA+RATIOC*CONl+RATICD*H) 

V- PHIET*PEA*G 

VI- Fl*ETIME/(ONE+V) 

V2-PHIET*PEA/  (ONE+V) 

tt******************************************************************* 

Trace  of  strain  increment  --  EMX 


DVOL-ZERO 
DO  10  Kl-l.NDI 
10  EM3L-IMDL+DSTRAN  (K1 ) 


Convert  stress  to  deviatoric  stress. 


DO  20  Kl-l.NDI 
20  STRESS(K1  )-STRESS(Kl )+PB 

Deviatoric  stress  times  strain  increment  --  SDSI 


SDSI-ZERO 
DO  30  Kl-l.NTENS 

30  SDSI-SDSI+STRESS(K1  )*DSTRAN(Kl) 


Effective  plastic  shear  strain  increment. 
*****************««***•**•**•••••*•**•***•«*•*•*******»***********«** 

DQAMPB-V1+V2  •  (AMJ*  sds  i  /TAUB-AKAPPA*RATICS*EWX,) 


Increments  DS  and  DIH. 
C*  **************** ****** 
DS-H*DGAMPB-DR 
DTH-CON 1 *DGAMPB 
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C  Constant*  for  Jacobian  and  stress  increment. 

V3-TW)*AMJB 

VI  l-AKAPPA-m>*THIRD*AMJB 
V4-Vll*DVOL 

V5-AKAPPA*ALPHA‘DTH/THIRD 

V6-AKAPPA*BETA*DGAMPB 

V7-V4-V5-V6 

V8-AMJ*DGAMPB/TAUB-(AMLJ-AMJB)*SDSI/TAUB**2 

V9-AMUB 

VI 2-  (V2  *AMJ*  *  2+AMJB -AVU)  /TAUB  *  *  2 

VI 3— V2  *AKAPPA*RATI03*AMJ/TAUB 

V14-V2  ’AKAPPA*  (AMJ/TAUB)  *(BETA+OON1  *ALPHA7  THIRD) 

V15— V2MKAPPA**2*PATIOB*  (BETA4C0N1  *ALPHA/THIRD) 

V16-V11-V15 

to***************** 

C  Calculate  the  Jacobian,  which  is  nonsymnetric  unless  V13-V14.  This 
C  is  generally  true  only  if  CMEGA-0,  BETA-0,  and  PDB-0.  Otherwise. 

C  ABAQUS  will  use  only  the  symmetric  part  of  the  Jacobian  unless  an 
C  unsymnetric  Jacobian  has  been  called  for  on  the  title  card. 

C  Note  that  STRESS  used  here  is  the  deviatoric  stress. 

DO  40  Kl-1 .NTENS 
DO  40  K2-1 ,  NTENS 

40  EDSEDECKl  ,K2)— V12*STRESS(Kl)*STRESS(K2) 

DO  50  Kl-1  ,NDI 
DO  50  K2-1 .NTENS 

DOS  DDE  (K1  ,K2)-EDSEOE(Kl  ,K2)-V14*STRESS(K2) 

50  DDSEDE(K2  ,K1  )-EDSEDE(K2  ,K1 ) -VI 3*  STRESS  (K2) 

DO  60  Kl-l.NDI 

DDSEDE(K1  ,K1  )-DDSDDE(Kl  ,K1  )+V3 
DO  60  K2-1.NDI 

60  DDSH*(K1,K2)-EDSDDE(K1,K2)+V16 
IF  (NSHR.GT.O)  THEN 
DO  70  K1-NDIP1  .NTENS 
70  DDSEDECK1  ,K1  )-DDSDDE(Kl  ,K1  )+V9 
END  IF 


Calculate  complete  STRESS  at  the  end  of  the  increment  using  the 
deviatoric  stress  at  the  beginning  of  the  increment. 


DO  80  Kl-l.NDI 

80  STRESS  (Kl)-(ONE-VS  )  *  STRESS  (K1  )+V3  *DSTRAN(K1  )+V7  -PB 
IF  (NSHR.GT.O)  THEN 
DO  90  K1-NDIP1  .NTENS 

90  STRESS  (K1 )-  (ONE  -V8  )  *  STRESS  (K1 )  +V9  *DSTRAN(K1 ) 

END  IF 


Update  the  state  variables. 

CALL  SINV  (STRESS, SINV1.SINV2) 
TAUB-EWAX1(SINV2/SQART3  .TAUIOL) 

PB— S1NV1 
TH-TH+DIH 
S-S+DS 

IF  (S.LT.SA)  S-SA 

CALL  GAMDOT  (TAUB ,  PB  ,TH ,  S ,  PROPS ,  NPROPS ,  STVART3 , 
1  F2,PEA,RATIOB,RATIOC,RATI0d7C0N2) 

STATEV(l)-S 
STATEV(2)-TH 
STATEV(3 )-F2 

STATEV(  4)-STATEV(  4  )+DGAMPB 
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C  Comparison  of  the  plastic  strain  rates  before  and  after  the  time 
C  increment.  To  be  used  by  the  automatic  integration  scheme  of  ABAQUS. 
C  Note  that  the  factor  SQART3  is  used  to  convert  shear  strain  to 
C  tensile  strain. 

DIFF-DTIME*EABS(F1-F2)/SQART3 
CEMAX-EM4X1  (CEMAX.DIFF) 


C  Check  magnitude  of  plastic  strain  increment  against  a  reference  level. 
C  This  allows  the  automatic  timestep  control  in  ABAQUS  to  limit  the  size 
C  of  the  plastic  strain  increment  using  the  variable  CEh&X. 
£*************»***********«******************»************************* 
PLSLMT- PROPS  ( 1 )  *  SQART3 
IF  (PLSlMr.GTT.ZERO)  THEN 
PLSCHK-(DGAMPB/PLSLMT)  *cetol 
CEM^X-EMOCl  (CEM4X,  PLSCHK) 

END  IF 
RETURN 
END 

0*******«**********»****************«*****t********s**t**************** 


~  SUBROUTINE  LMPROP  (AMJ, AKAPPA, ALPHA., RHD.C.TH, PROPS .NPROPS) 
C  Determine  the  following  constants  for  the  material: 


C  AMJ  -  shear  modulus  (suggested  units:  N/m2) 

C  AKAPPA  -  bulk  modulus  (suggested  units:  N/m2) 

C  ALPHA.  --  thermal  expansion  coefficient  (suggested  units:  1/Deg, K) 
C  RHO  ----  mass  density  (suggested  units:  kg/m3) 

CC - specific  heat  (suggested  units:  Joules/kg/Deg.K) 


C  In  general,  these  properties  are  functions  of  temperature,  but  in 
C  the  present  case  they  are  input  directly,  assuming  no  temperature 
C  dependence.  Additional  entries  in  PROPS  could  be  used  to  define  them 
C  as  functions  of  temperature. 

C**i**«*s»***i<*********««»«*«»***t **»•*«*•**>*******•*****«****»****** 

IMPLICIT  REAL*8(A-H,0-Z) 


DIMENSION  PROPS (NPROPS) 

*******»*****«****•*****«**•*<************************•*«************ 

AMJ-PROPS(6) 

AKAPPA-PROPS ( 7 ) 

ALPHA-PROPS ( 8 ) 

RHD-PROPS(9) 

OPROPS(IO) 

*<*(********«s***«***s******a*****»********s************************* 

RETURN 

END 


C* ********************************************************************  * 

SUBROUTINE  SBETA  (TAUB ,  PB , TH , S ,  BETA) 

q*  a******************************************************************** 

C  Subroutine  SBETA  determines  the  plastic  dilatancy  factor. 
q* *«***«t*********t****a****(***************************<************** 

IMPLICIT  REAL*8(A-H,0-Z) 

BETA-O.  ODfO 

. . . . . 

RETURN 

END 

f'**««***«««**««**«**«*****«*t*«S*******S*********«*******«*****«******* 
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SUBROUTINE  GfflGOrr  ( TAUB ,  PB , TH ,  S ,  PROPS . NPROPS ,  SQART3 , 

1  F,PEA,RATI0B,RA.TI0C,RATI0D,Cm2) 

Subroutine  QA M30T  determinei  the  equivalent  plastic  shear  strain  rate 
and  its  derivatives  PDA,PDB,PDC,PDD  with  respect  to  TAUB, PB, TO, and  S, 
respectively.  Note  that  the  following  derivative  terms  are  returned: 
RATIOB  -  PDB/PDA 
RATIOC  -  PDC/PDA 
RATIOD  -  PCD/ PDA 

£**«•**»»»«*»**•****•*•*••*•*«•»**»•*•**«****«**»»»**»•«*•***«*******(* 
IMPLICIT  REAL*8(A-H,0-Z) 

C*  ********************************************************************* 

CChMWOONSTS/PI,SIN60,COS60,KCROS2(3),KCROS3(3), ZERO, LZERO, LONE, 
1  ONE.TOO,  HALF  ,ABIG,ASMALL,BCBIG,IjOCSHR(2, 3) , THIRD,  PRECIS,  BLANK 

Q*  *****««*•*•«•**•«•••*«•«••••««**•**«*«»*«********•«*»**»**»********** 

DIMENSION  PROPS (NPROPS) 

. . . 

C  Material  parameters  defining  the  equivalent  plastic  shear  strain  rate 
C  PROPS(ll)-A  —  pre-exponential  factor 
C  suggested  units  are:  l./sec 

C  PROPS(12)-Q/R  ---  activation  energy  divided  by  universal  gas  constant 
C  suggested  units  for  Q  are:  kJ/mol 

C  suggested  units  for  R  are:  8.314D-3  kJ/mol/Deg.K 

C  PROPSll3)-l . /AM  --  1. /strain  rate  sensitivity 
C  PROPS ( 14 )-PALPHA  --  pressure  sensitivity  parameter 
C  suggested  units  are:  l/(N/m2) 

C  Note:  PALPHA-2*(IMJ/DPB)/ML!0  ,  with 
C  MJO  --  shear  modulus  at  zero  pressure 

C  EMJ/DPB  --  pressure  dependency  of  shear  modulus 

. . . . . 

AL  -  3.250D0 

OON2-SQART3  *  PROPS (11) *DEXP ( - PROPS ( 1 2 ) /TO) 

C*** 

F-OON2*(DSINH(TAUB*AL/(S*(ONE+PROPS(14)*PB))))**PROPS(13) 


C  Note  that  a  tolerance  was  set  on  the  minimum  value  of  TAUB  allowed, 
C  to  avoid  calculat ional  difficulties. 


SP  -  S  * (ONE+PROPS ( 14 ) *PB ) 

PDA-F*PROPS(l3)»DCOSH(TAUB*AL/SP)*AL/DSINH(TAUB*AL/SP)/SP 
RATIOB-0 . 0 

RATIOC-PROPS  ( 1 2  )  *  S  «DTANH(TAUB  *AL/S  )  /  (AL*  PROPS  (13)  *TO*  *  2  ) 
RATIOD— TAUB/S 

C***»* . *.«»**«. . . 

RETURN 

END 

C** 

C*  **********************************************  ******************* 

SUBROUTINE  SDCTT  (TAUB,  PB,  TO,  S,AMJ,CDN2,  PROPS,  NPROPS,  SQART3.F, 
1  H.RDOT.SA) 

c*  ********•*««»••«»•**•*•*•*•*•*•**••*(•****•****•****************# 
C  'This  subroutine  determines  the  hardening  rate  H  and  the  static 
C  testoration  rate  RDOT. 


IMPLICIT  REAL*8(A-H,0-Z) 

. . a**************************************** 

OCMvON/OONSTS/PI  ,SIN60,COS60  ,KCROS2(3)  ,KCROS3(3)  , ZERO, LZERO, LONE, 
1  ONE, TOO, HALF  ,ABIG,A£MALL,BCBIG,IjOCSHR(2, 3), THIRD, PRECIS, BLANK 

DIMENSION  PROPS  (NPROPS) 
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C  Material  parameters  determining  the  rate  of  hardening: 

C  PROPS( 15)-H0  --  pre-multiplier  in  hardening  rate  equation 
C  suggested  units:  N/m2 

C  PROPS(16)-SB  --  pre-multiplier  for  saturation  value  of  S 
C  suggested  units:  N/m2 

C  PROPS(17)«AN  --  eiponent  in  eqn.  for  saturation  value  of  S 
C  PR0PS(18)-B  —  pre-multiplier  for  rate  of  static  restoration 
C  suggested  units:  N/m2 
C  PROPS(19)-SA  --  annealed  value  of  S. 

C  suggested  units:  N/m2 

C  PROPS! 20 )-P  ---  eiponent  of  ((S/SA)-1.)  in  RDCT 
C  The  following  constant  was  calculated  in  GAMX7T. 

C  COM2  -  A*DEXP(  -Q/CR'IlO  ) 

C  Note  that  HO.SB.B,  and  SA  should  be  from  tensile  test  data  because 
C  a  conversion  to  pure  shear  data  is  done  in  this  subroutine. 

C*  ********************************************************************* 
C  First  calculate  H. 

£******»*********•******•*•**«******************************«*******»**. 

SS-  ( PROPS (16) /SQ\RT3 ) * ( F/CDN2 ) * ‘PROPS (17) 

IF  (S.LE.SS)  THEN 

H-THIRD*PROPS( 15) * (ONE-S/SS )  “  1 .  SDO 

ELSE 

H-ZERO 

END  IF 

C  Determine  the  annealed  value  of  S.  Although  in  general  it  is  a 
C  function  of  temperature ,  in  the  present  version  it  is  input  directly 
C  and  considered  constant. 

“““““““““““““******************* 

SA-PROPS  (19)  /SQART3 

C“““““““““““““““““*““““““*“““*“************* 
C  Calculate  RDOT. 

£****(«*****(•«** ****•***«*•« «*****»»••**•****«*» **t ******************* 

SD-(S/SA)-ONE 
IF  (SD.GT.ZERO)  THEN 

RDOT«AMJ*THIRD*PROPS  (18)  *(I»12  *SD‘  ‘PROPS  (  20  ) 

ELSE 

RDOT-ZERO 
END  IF 

RETURN 

END 
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‘HEADING 

LOAD  BOUNDARY  CONDITIONS  -  20000  N 
‘NODE 
1,0. 0,0.0 
2,0.00635,0.0 
3,0.00635,0.009525 
4,0.0,0.009525 
•ELEMENT ,  TYPE-CAX4H 
1,1, 2, 3, 4 
‘BOUNDARY 
1,1,2 
2,2 

4.1 

‘EQUATION 

2 

3,2,1. ,4, 2,-1. 

‘MATERIAL 

‘USER  MATERIAL, CONSTANTS-20 

.01D0  .0.75D0  ,0.D0  ,  60.0D6  ,1073.14D0  ,32.0D9  .1.59D11  .16.6D-6 

7.3D3  ,0.44D3,6.346D11 .3.7569D4  .5.1125D0  .0.0D0  .3774.D6  .125.1D6 

.06869DO.O.ODO  ,  60.OD6.1.0D0 

‘DEPVAR 
4 

‘AMPLITUDE , VALUE-A  ,TIME-V ,  NAME-TOPLQAD 
0.,  0.,  10.,  -20000,  20.  ,  -20000,  30.  ,  0. 

‘STEP  .NLGECM,  AMP-RAMP ,  INC-1500  .CYCLE-6 ,  SU1MAX 

*VISCO,PTOL-l  .O.CETOL-l . OE-5 

0.0001,30.0, , .06 

‘QjQAD.AMPLITUDE-TOPLQAD 

4,2,1. 

*EL  FILE.FREQ-1 

2.1 

2,1,1 

*EL  PRINT, FREQ-1000 
2,1 
2,1,1 

‘NODE  PRINT, FREQ-1000 
‘PRINT, RES I DUAL-NO, FREQ- 100 
‘END  STEP 
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‘HEADING 

ALIMINTM  GRADIENT  SPECIMEN 

**  NEUTRAL  FILE  GENERATED  ON:  26-MfAR- 86  09:38:09  PATRAN  VERSION:  1.63 

**  SUMMARY  DATA: 

**  52  NODES,  36  ELEMENTS,  0  MATERIALS,  0  PHYSICAL  PROPERTIES 

•NODE,  NSET-CFOOOO 

1  0.0  0.0 

4  0.01270  0.0 

49  0.0  0.04041 

52  0.00374  0.04041 

•NGEN , NSET-LS IDE 

1.49.4 

•NGEN, NSET-RS IDE 

4.52.4 
•NGEN 

1.4.1 

5.8.1 

9.12.1 

13.16.1 

17.20.1 

21.24.1 

25.28.1 

29.32.1 

33.36.1 

37.40.1 

41.44.1 

45.48.1 

49.52.1 

•NSET , NSET-TOP 
49,50,51,52 
•ELEMENT,  TYPE-CAX4H 
1,  1,  2,  6,  5 


•ELGEN 

1,3,1,1,12,4,3 

•BOUNDARY 

1 

1  2 

•OOOE+OO 

2 

2 

. OOOE+OO 

3 

2 

•OOOE+OO 

4 

2 

.OOOE+OO 

5 

1 

.OOOE+OO 

9 

1 

•OOOE+OO 

13 

1 

.OOOE+OO 

17 

1 

. OOOE+OO 

21 

1 

.OOOE+OO 

25 

1 

.OOOE+OO 

29 

1 

.OOOE+OO 

33 

1 

.OOOE+OO 

37 

1 

.OOOE+OO 

41 

1 

.OOOE+OO 

45 

1 

.OOOE+OO 

49 

1 

.OOOE+OO 

•EQUATION 

2 

52.2.1. . 49. 2,-1. 

2 

51.2.1. . 49. 2,-1. 

2 

50.2.1..  49.2, -1. 

•MATERIAL 

•USER  MATERIAL,  OONSTANTS-20 

.01D0,  0.75D0  ,0 .DO  ,  35.0D6  ,  573.16D0  ,25.0D9  ,0.70D11  .22.0D-6 
2.753  , 1 . 00D3  .1.906D7  .2.1086D4  .4.2831D0  .O.ODO  .1115.D6  ,18.9D6 
.07049D0.0.0D0  ,  35.0D6  .l.ODO 
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•DEPVAR 

4 

‘PLOT 

•DRAW 

•STEP  .NLGECM,  AMP-RAMP ,  INC-1 ,  CYCLE-6 ,  SUBMAX 
*  VI SCX) ,  PTOL- 1 . 0 ,  CETOL- 1 .  E  -  4 
0 . 0000065 , 0 . 0000065 , 0 . 0000065 
•BOUNDARY 

TOP  ,2,, -0.00000001 
•NODE  PRINT 
•PRINT, RESIDUAL-NO 
•END  STEP 

•RESTART  ,m  ITE ,  FREQ-300 

•STEP,  NLGECM,  AMP-RAMP,  IN0300,  CYCLE-6,  SUBMAX 
•VI SCO, PTOL- 5.0,  CETOL-1.0E-4 
0.001,  8.0, ,.08 
•BOUNDARY 
TOP, 2, ,-0.0080 
•EL  FILE, FREQ-49 
2  2  11 

2  2  11 

2  2  2  2 

•NODE  FILE, FREQ-49 

2  11112 
•EL  PR  I  NT,  DEPVAR,  FREQ-4  9 
2,1 
2,1,1 

•NODE  PRINT, FREQ-49 
•PRINT, RESIDUAL-NO 

••  PLOTTING  SECTION 

•PLOT .FREQ-49 
•DISPLACED 
1,1.0 
•CONTOUR 

•*  VON  M3SES  STRESS 

9.10 

•*  INTERNAL  VARIABLE 

81.10 
•END  STEP 
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DISTRIBUTION  LIST 


No.  of 

Copies _ To  _ 

Office  of  the  Under  Secretary  of  Defense  for  Research  and  Engineering, 

1  The  Pentagon,  Washington,  DC  20301 

Commander,  U.S.  Army  Materiel  Command,  5001  Eisenhower  Avenue, 

Alexandria,  VA  22333-0001 
1  ATTN:  AMCLD 

Commander,  U.S.  Army  Laboratory  Command,  2800  Powder  Mill  Road,  Adelphi, 

MD  20783-1145 

1  ATTN:  AMSLC-IM-TL 

Commander,  Defense  Technical  Information  Center,  Cameron  Station, 

Building  5,  5010  Duke  Street,  Alexandria,  V A  22304-6145 

2  ATTN:  DTIC-FDAC 

Metals  and  Ceramics  Information  Center,  Battelle  Columbus  Laboratories, 

1  505  King  Avenue,  Columbus,  OH  43201 

Commander,  Army  Research  Office,  P.0.  Box  12211,  Research  Triangle  Park, 

NC  27709-2211 

1  ATTN:  Information  Processing  Office 

Commander,  U.S.  Army  Electronics  Technology  and  Devices  Laboratory, 

Fort  Monmouth,  NJ  07703-5000 
1  ATTN:  SLCET-DT 

Commander  U.S.  Army  Missile  Command,  Redstone  Arsenal,  AL  35898-5247 
1  ATTN:  AMSMI-RD-ST 

1  Technical  Library 

Commander,  U.S.  Army  Armament,  Munitions  and  Chemical  Command,  Dover,  NJ 
07801 

2  ATTN:  SMCAR-TDC 

Commander,  U.S.  Army  Natick  Research,  Development  and  Engineering  Center, 
Natick,  MA  01760 
1  ATTN:  Technical  Library 

Commander,  U.S.  Army  Tank-Automotive  Command,  Warren,  MI  48397-5000 
1  ATTN:  AMSTA-R 

Commander,  U.S.  Army  Engineer  Waterways  Experiment  Station,  P.0.  Box  631, 
Vicksburg,  MS  39180 
1  ATTN:  Research  Center  Library 

Director,  U.S.  Army  Ballistic  Research  Laboratory,  Aberdeen  Proving  Ground, 
MD  21005 

1  ATTN:  SLCBR-DD-T  (STINFO) 


To 


No.  of 
Copies 

Director,  Benet  Weapons  Laboratory,  LCWSL,  USA  AMCCOM,  Watervliet,  NY 
12189 

1  ATTN:  AMSMC-LCB-TL 

Commander,  U.S.  Army  Foreign  Science  and  Technology  Center,  220  7th  Street 
N.E.,  Charlottesville,  V A  22901-5396 
1  ATTN:  AIAST-RA-ST 

Director,  Eustis  Directorate,  U.S.  Army  Air  Mobility  Research  and 
Development  Laboratory,  Fort  Eustis,  VA  23604-5577 
1  ATTN:  SAVDL-E-MOS  (AVSCOM) 

Director,  Langley  Directorate,  U.S.  Army  Air  Mobility  Research  and 
Development  Laboratory,  NASA-Langley  Research  Center,  Hampton,  VA  23665 
1  ATTN:  Aerostructures  Directorate 

Naval  Research  Laboratory,  Washington,  DC  20375 
1  ATTN:  Code  6830 

Office  of  Naval  Research,  800  North  Quincy  Street,  Arlington,  VA 
-22217-5000 

1  ATTN:  Mechanics  Division,  Code  1132-SM 

Naval  Air  Development  Center,  Warminster,  PA  18974-5000 
1  ATTN:  Code  6064 
1  AVCSTD/6043 

U.S.  Navy  David  Taylor  Research  Center,  Bethesda,  MD  20084 
1  ATTN:  Code  172 

U.S.  Air  Force  Office  of  Scientific  Research,  Bolling  Air  Force  Base, 
Washington,  DC  20332 
1  ATTN:  Mechanics  Division 

Commander,  U.S.  Air  Force  Materials  Laboratory,  Wright-Patterson  Air  Force 
Base,  OH  45433 
1  ATTN:  AFWAL/MLLN 

National  Aeronautics  and  Space  Administration,  Marshall  Space  Flight 
Center,  Huntsville,  AL  35812 
1  ATTN:  EHOi,  Dir,  M&P  Lab 

Committee  on  Marine  Structures,  Marine  Board,  National  Research  Council, 

1  2101  Constitution  Avenue,  N.W.,  Washington,  DC  20418 

Director,  U.S.  Army  Materials  Technology  Laboratory,  Watertown,  MA 
02172-0001 
1  ATTN:  SLCMT-TML 

1  SLCMT-PR 

1  SLCMT-IMA-T 
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